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(Unclassified  Foreword) 
FOREWORD 


714*  1.  ?»«*  -  ©•*•♦  tt  nnv««>iB^  the  work  per¬ 

formed  by  Thiokol  Chemical  Corporation,  Reaction  motors  Division, 
Denvllle,  New  Jersey  under  an  extension  to  Air  Force  Contract 
AFG4 (611) -9946.  The  internal  report  number  is  Report  REG 
5801-F  Part  II.  The  work  was  administered  under  the  direction 
of  Mr.  X.  Rimer,  Rocket  Propulsion  Laboratory,  Edwards  Air  Force 
Base  under  the  following  project  and  tawk  numbers: 

Project  Ho.  3058,  Task  Ho.  3058560 
Air  Force  Program  Structure  Ho.  750  0 

The  Research  effort  reported  herein  was  ’.onducted  during 
the  period  15  April  1965  to  15  October  1965  >n  RED  Project  5801. 
The  work  previously  performed  un^er  Part  I  cf  the  program 
covering  the  period  1  May  1964  to  28  February  1965  is  reported 
in  AFRPL-TR-65-105. 

Principal  investigators  of  the  effort  described  herein  are 
Messrs.  B.  R.  Dawson,  T.  F.  Seamans,  and  to.  M.  Vanpee.  Other 
contributors  to  the  program  were  Messrs.  R.  dornslein,  R.  Storms, 
and  f.  Sutton.  Dr.  Vito  Agoste,  Professor  of  Aerospace  Engin¬ 
eering,  Polytechnic  Institute  of  Brooklyn,  served  as  a  consultant 
for  the  establishment  of  the  vaporization  model  and  contributed 
significantly  to  its  development.  The  Project  Leader  was 
Mr.  A.  D.  Corbett  and  the  Program  Manager  was  Mr.  S.  J,  Tunkel. 

This  report  contains  no  classified  Information  extracted 

from  Other  classified  documents. 

This  technical  report  has  been  reviewed  and  is  approved. 

Mary  M,  Racovich 
AF  Contracting  Officer 
Directorate  of  Procurement 
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(Unclassified  Abstract) 
ABSTRACT 


By  theoretical  and  experimental  means,  a  mathematical  model 
of  hyperbolic  ignition  in  reaction  control  systems  has  been 
developed  based  on  physical  kinetics  of  droplet  evaporation  and 
ktzoticr  cf  4 “r icier  *  emotions.  Reasonable  agreement 
between  theoretical  and  experimental  ignition  delays  are 
obtained.  The  dominant  reactions  of  N2O4  with  the  hydrazine 
fuels  at  reduced  pressures  are  found  to  be  thermal,  gas  phase 
reactions  which  are  bimolecular  and  have  low  activation  energies 
and  pre-exponential  factors.  A  pre-ignition  reaction  product 
was  found  which  is  a  clear,  yellow,  viscous  liquid  with  a  very 
low  vapor  pressure.  The  ’’adduct"  has  the  characteristics  of 
a  monopropellant  and  contains  considerable  energy.  Its  relation 
to  pressure  spiking  during  engine  start  transients  is  as  yet 
undetermined.  Of  six  additives  tested,  only  furfuryl  alcohol 
had  a  significant  beneficial  effect  on  the  ignition  character¬ 
istics  of  N204/10IH. 

The  influence  of  thrust  chamber  design  parameters  on  the 
ignition  delay  and  pressure  transients  of  Compound  A  with  N2H4, 
UDHH  and  MHF-5  and  also  gaseous  *2/H2  was  investigated  The 
results  are  reported.  Very  short  ignition  delays  and  smooth 
pressure  transients  were  obtained  with  the  F2/H2  combination 
regardless  of  chamber  configuration.  Also,  short  rise  times 
to  90%  steady-state  chamber  pressure  resulted,  However,  a 
dependency  c-  chamber  L*  only  exists. 
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INTRODUCTION 


Thiokol  Chemical  Corporation,  Reaction  Motors  Division  has 
been  investigating  the  ignition  characteristics  at  reduced 
pressures  of  6  number  of  hypergols  under  Air  Force  Contract 
AF04 (611) -9941? .  The  Final  Report  -  Part  I  covering  the  work 
performed  during  the  initial  ten-month  program  is  contained  in 
Report  AFRPL- Til-65-105,  dated  July  1965.  These  studies  included 
an  investigation  of  hypergolic  ignition  characteristics  in 
(a)  an  unconfined  impinging  stream  apparatus  and  (b)  50-lb. 
thrust  chambers  at  simulated  altitudes  for  the  purpose  of 
determining  attitude  control  engine  design  criteria  to  minimize 
ignition  delays  and  to  eliminate  pressure  spikes  in  the  start 
transients. 

In  brief,  the  results  of  the  initial  ten-month  program 
indicate  that: 

1)  ignition  delays  of  N204/hydrazine-type  fuels  are 
strongly  pressure  dependent, 

2)  impinging  stream  injector  parameters  -  impingement 
angle  and  length,  injection  velocity,  type  of 
manifolding  -  have  a  negligible  effect  on  ignition 
delays, 

3)  for  constant  propellant  flowrate,  design  chamber 
pressure,  i.e.  nozzle  throat  area,  significantly 
affects  ignition  delay  and  rise  time  to  steady 
state  chamber  pressure, 

4)  for  constant  propellant  flowrate  and  constant 
design  chamber  pressure  (throat  area)  thrust 
chamber  geometry  -  L*  and  contraction  ratio  -  has 
little  effect  on  ignition  delay  and  rise  time, 

5)  pressure  spikes  during  start  transients  are  of 
random  occurrence  and  magnitude  and  are  not  affected 
appreciably  by  any  of  the  engine  design  parameters 
studied  (design  chamber  pressure,  L*,  and  contraction 
ratio) , 

and  6)  a  more  fundamental  approach  is  needed  to  understand 
and  then  correct  the  ignition  delay-pressure  spiking 
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problea  that  exists  with  tha  KgOi/hydraslne-tj rpe  fuala 
combinations. 

Tba  reader  la  referred  to  tha  final  roport  of  tha  ten-month 
program  (lef.  1 '  ter  a  data  Had  discussion  of  the  above  findings. 

A  six-months  continuation  program  vaa  under taken  subsequently 
for  tha  purpose  of  construct inf  a  fundamental  Mathematical  modal 
for  tba  ignltloii  of  hypergolio  propellants  at  reduced  presau,  es 
is  reaction  control  systdmr.  The  program  was  performed  in  three 
major  tasks.  Tha  specific  objectives  of  each  task  are  described 
as  follovs: 

Tush  1  -  Measurement  of  Saaction  Antes 

—  —  m...-  -  ■  — — T . . . . 

Ton  objectives  of  this  task  are  to  measure  the  activation 
energy  and  order  of  reaction  for  the  folloTing  propellant 
combinations  and  construct  a  mathematical  model  for  igniti 
of  tisese  propellants  at  reduced  pressures  in  altitude  control 
engines. 


W304/b0-50  s2o4/imh 

f2/h a 

Task  II  -  Evaluation  of  Additives  to  Reduce  activation  Energy 

The  objective  of  thia  task  is  to  evaluate  propellant 
additives  fer  the  reduction  of  activation  energy  requirements 
for  Mg&f/lOffl  and  Jljdi/SC-SO  using  the  experimental  apparatus 
developed  for  Task  I  above. 

Tusk  ill  -  Tkruut  Chamber  Design  Parameter  Study 

muimmamm— — MMMfcm— — amm— — ev  1  m  m—M*— ee— aw—*— 


The  objective  of  this  tack  la  to  evaluate  the  influence  of 
thrufit  chamber  design  parameters  on  the  -gnition  delays  and 
atari  transients  of  Compound  A/K^^  tod  gaseous  Fg/Bj.  The 
specific  parameters*  of  interest  are  characteristic  length, 
chamber  pressure^  contraction  ratio,  and  propellant  lead, 1 *g. 
$5*ts  me*o  made  at  a  thrust  leval  of  '  0  ~i  ha  over  the  rtage  c-f 
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Characteristic  length  (L*) 

Chamber  pressure 
Contraction  ratio  ^ 

Oxidizer  lead 
""“-’•t  (noainal) 

Ambient  pressure 

The  results  of  the  experimental  and 
are  discussed  in  the  following  sections. 

V  ‘ 


G-50  in. 

20-200  psia 
1.5-8 

+2  to  -2  msec 
50  lbs. 

<0.2  psia 
theoretical  program 
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TI 

SUMMARY 


The  objective  of  the  program  reported  herein  van  to  construct 
a  fundamental  mathematical  model  for  the  ignition  of  hypergolic 
pr^pellan  a  at  reduced  pressures.  This  program  is  a  sir-month 
extension  of  an  initial  ten-month  program  to  estaMish  attitude 
control  engine  design  criteria  which  would  minimise  ignition 
delay  and  eliminate  pressure  spikes  in  the  start  transients 
with  byp^r golic  propellants  at  reduced  pressures.  The  results 
of  the*  initial  ten-month  program  are  described  in  detail  in 
flepcrtJ  AFB^l— TR-65-105.  Results  of  the  initial  prog,  am  indicated 
that  fhe  Conventional  injector  and  thrust  chamber  configuration 
parametenm  had  little  influence  on  ignition  and  tnat  t  more 
fundamental  investigation  of  the  pre-ignition  reaction  rates 
and  chemical  kinetics  was  required  in  order  to  develot  a 
suitable  model  for  ignition  of  hyper go lie  propellant?  at  reduced 
pressures . 

This  report  is  the  final  report  of  the  six-nonth  extension 
program  which  was  performed  in  three  separate  tasks  in  order 
to  accomplish  the  objectives.  The  purpose  of  Task  I  was  to 
measure  the  pre-ignition  reaction  rates  for  h204/^nDI,H»  W2Q4/MMH, 
U2Q4/R2H4  and  ?2/®2  *nd  incorporate  these  data  Into  a  fundamental 
mathematical  model  for  the  ignition  ~t  h, ,  rrgolic  propellants 
at  reducod  pressures  in  reaction  control  systems.  The  purpose 
of  Task  D.  was  to  evaluate  selected  additives  for  tb«  1 “Auction 
of  activation  energy  requirements  for  and  HjOj/SO-SQ. 

The  purpose  of  Task  III  was  to  evaluate  the  influence  of  thrust 
eh amber  design  parameter  on  the  ignition  delay  and  pressure 
transients  of  Compound  A/T^Hr  and  gaseous  F2/H2.  The  results 
of  these  ta^fcs  are  summarized  below. 

The  dsvfi'lc' went  of  a  mathematical  model  of  hyperbolic 
ignition  in  reaction  control  systems  followed  two  pates.  The 
d#t*?nainatiou  of  the  pressure  history  in  a”  engine  due  to 
propellant  vaporization  and  the  determination  of  the  chemical 
kinetics  of  ignition  reactions  of  the  various  hypergolic 
combinations.  It  is  necessary  to  consider  bo  ..ft  parts  of  th? 
problem  since  reaction  ti-.es  for  the  hydrazine- type  fuels  with 
nitrog-'-:  t*  reside  at  the  pressure  encountered  oaring  engine 
start -up  are  ca*pars.fele  to  th^  times  required  for  the  pressure 
ia  the  thrust  -chamber  to  reach  levels  at  which  ignition  can 
occur .  ‘.fhe  hyporgolte  ignition  model  results  free  z  merger  of 
the  two  paths. 
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A  computer  program  was  developed  for  predicting  chamber 
pressurization  due  to  propellant  vaporization  talcing  into  account 
the  rate  of  propellant  injection,  the  kinetics  of  propellant 
droplet  evaporation  and  vapor  condensation,  vapor-drop  heat 
transfer,  drop  residence  time,  thrust  chamber  geometry,  and 
propellant  vapor  efflux  through  the  nozzle.  A  three  drop-size 
spray  model  is  used  to  represent  mathematically  the  actual 
propellant  spray  vhich  forms  in  the  thrust  chamber.  The 
principal  quantities,  whose  time-dependent  values  are  calculated, 
are: 

(a',  mass  evaporation  from  each  spray  drop 

(b)  radius  of  each  drop 

(c)  temperature  of  each  drop 

(d)  fraction  of  each  drop  that  is  frozen 

(e)  vapor  mass  flow  through  the  nozzle 

(f)  vapor  mass  condensing  on  the  chamber  walls 

(g)  gas  temperature 
y''' (h)  gas  pressure. 

Using  carfcoa  tetrachloride  as  the  working  fluid  since  its 
properties  are  well  known,  it  was  found  that  the  vapor-drop 
system  is  non-ad iaba tic.  Despite  the  short  times  involved  in 
chamber  pressurization  due  to  propellant  vaporization,  heat 
addition  to  the  vapor-drop  system  occurs  when  drops  strike 
the  chamber  walls.  Using  a  crude  model  to  account  for  the 
heat  addition  to  the  vapor-drop  system,  agreement  to  within 
10%  was  obtained  between  calculated  and  experimental  curves 
of  chamber  pressurization. 

The  computer  program  permitted  an  evaluation  to  be  made 
of  two  potential  problem  areas:  the  effect  on  chamber  pres¬ 
surization  rates  of  accommodation  coefficient  and  of  spray 
drop  sizes.  Using  two  propellants  individually  (one,  CC14; 
the  other?  hypothetical  propellant  differing  only  in  its 
accommodation  coefficient)  it  was  found  that  the  effect  of 
accommodation  coefficient  oa  chamber  pressurization  is  much 
less  than  linear.  The  direct  effect  of  a  low  accommodation 
coefficient-decreased  evaporation  rate  at  a  given  drop 
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temperature  -  is  offset  to  a  great  extent  by  the  resulting 
slower  cooling  rate  of  the  drop.  Thus,  the  vapor  pressure  of 
the  drop  is  higher  than  it  otherwise  would  be  and  this  results 
in  relatively  sore  vaporization.  A  50%  decrease  in  accommo¬ 
dation  coefficient  caused  only  a  10  to  15%  decrease  in  chamber 
pressure  at  any  time.  Since  accommodation  coefficients  are  in 
general  unknown  and  are  extremely  difficult  to  measure,  it  is 
fortunate  their  effect  on  chamber  pressurization  rates  is  small. 

The  effect  of  spray  drop  sizes  on  chamber  pressurization 
was  also  fouud  to  be  small.  Two  computer  runs  were  made  differ¬ 
ing  only  in  the  initial  values  of  radii  of  the  three  drop-size 
sprays.  The  mass-median  drop  in  one  case  was  75  microns,  in  the 
other  50  microns.  The  resulting  pressure  curves  differed  by 
only  7%.  Although  the  smaller  drops  mean  a  greater  total  drop 
surface  area  since  there  are  more  drops  for  a  given  flowrate, 
these  smaller  drops  evaporatively  cool  at  a  faster  rate  which 
leads  to  less  vaporization  due  to  lower  vapor  pressure.  Due 
to  the  uncertainties  in  actual  drop  sizes  of  propellant  sprays, 
this  finding,  too,  is  fortunate. 

Very  short  exposure  photographs  were  taken  of  propellant 
streams  injected  into  a  low  pressure  environment.  The  photo¬ 
graphs  show  not  only  that  ambient  pressure  strongly  affects 
the  characteristics  of  the  sprays  but  also  that  the  breakup 
of  the  injected  stream  into  drops  occurs  quickly  at  low 
ambient  pressures,  up  to  1  1/2  msec  depending  on  the  propellant 
and  injection  velocity. 

For  the  more  volatile  propellants  such  as  N204  and  Compound 
A,  the  calculated  and  experimental  chamber  pressurization  curves 
indicate  that  a  significant  amount  of  time  is  required  for  the 
flow,  once  initiated,  to  build-up  to  the  full  flowrate. 
Appreciable  "flashing"  of  the  propellant  within  the  injector 
volume  slows  the  arrival  of  all-liquid  flow  at  the  injector  face. 
The  reduced  mass  flow  causes  a  slower  pressure  rise  in  the 
chamber  than  would  otherwise  be  the  case.  The  computer  program 
at  present  does  not  include  expressions  to  account  for  transient 
flow  upon  propellant  valve  opening.  With  less  volatile  propel¬ 
lants,  i.e.,  the  hydrazine-type  fuels,  full  flowrates  are 
reached  very  quickly  if  fast-acting  valves  are  used  and  injector 
volume  is  kept  to  a  minimum. 

The  kinetic  factors  required  in  the  hypergolic  ignition 
model  to  describe  the  ignition  reactions  were  obtained  through 
an  application  of  the  theory  of  thermal  explosions  to  experi¬ 
mentally  measure  ignition  pressure  limits.  The  overall  order 
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I  of  reaction  was  determined  from  the  variation  of  ignition 

I  pressure  limits  with  reactor  size,  the  activation  energy  from 

the  effect  of  initial  temperature,  and  the  partial  reaction 
orders  from  composition  effects.  The  pre-exponential  factor 
was  calculated  in  two  ways:  (1)  from  the  theory  of  thermal 
explosions  (steady-state  approximation  equation)  using  the 
experimentally  determined  values  of  overall  reaction  order, 
partial  reaction  orders  and  activation  energy;  and  (2)  from 
a  measured  axial  temperature  profile  in  the  premixed  reactant 
stream  at  low  pressures.  The  good  agreement  between  the  two 
determinations  of  the  pre-exponent ial  factor  substantiates 
the  applicability  of  the  theory  of  thermal  explosions  to  the 
ignition  reactions  of  the  hydrazine- type  fuels  with  nitrogen 
tetroxide.  The  reactions  are  found  to  be  bimolecular  with* 
low  activation  energies  (7.2  kcal/mole  for  UDMH  and  5.2  heal/ 
mole  for  MMH)  and  low  pre-exponential  factors. 

During  the  experimental  determinations  of  the  kinetic 
factors,  the  formation  of  a  pre-ignition  reaction  product  was 
observed  for  each  of  the  four  hydrazine-fuel  combinations 
tested.  The  adducts  formed  at  pressures  .below  minimum 
ignition  pressures,  condensed  on  the  walls  of  the  reactor,  and 
collected  as  a  liquid  at  the  base  of  the  experimental  apparatus. 
In  general,  the  adducts  are  clear,  yellow,  viscous  liquids 
with  very  low  vapor  pressures  and  are  stable  at  room  temperatures 
i  and  pressures.  Preliminary  analyses  of  the  adduct  formed  from 

i  .  premixed  MMH  and  1K>2  vapors  at  5  mm  Hg  indicate  the  adduct  is 

a  simple  additive  product  which  has  the  characteristics  of  a 
monopropellant  and  which  contains  considerable  energy.  Its 
role  in  relation  to  pressure  spiking  during  start  transients 
in  reaction  control  systems  is  as  yet  undetermined. 

The  theoretically  derived  Ignition  delay  equation  (bimolec¬ 
ular  ignition  reaction)  is: 

r*»-  ;&(■* 


where  P«r  “  Pf  +  Pox  and  £  "  r-—  (see  nomenclature  for 

*  meaning  of  symbols) 

The  reactant  partial  pressures,  pf  and  Pox»  are  obtained 
individually  by  the  analysis  for  chamber  pressurization  due 
to  propellant  vaporization. 
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The  above  equation  assumes  that  gas  phase  ignition  reactions 
-  are  the  dominant  reactions.  Hie  justification  for  the  assuaption 
coaes  froa  both  experiaental  and  theoretical  results.  Experi- 
■ental  ignition  delays  of  several  hydrazine  fuels  with  N204  in 
an  iapicging  streaa  apparatus  indicate  that  fuel  vaporization 
is  the  controlling  process  at  one  ataosphere  aabient  pressure, 
thereby  iaplying  gas  phase  reactions.  Also,  theoretical  cal¬ 
culations  show  that  the  collision  frequency  between  unlike 
gas  molecules  is  4000  tiaes  greater  than  the  frequency  of  gas- 
drop  surface  collisions  at  coaparable  conditions. 

Pressure-ignition  delay  curves,  calculated  by  the  above 
equation  for  N3O4/UDMH  and  N204/MMH,  show  the  effects  of  both 
teaperature  ana  vapor  phase  composition  on  ignition  delays. 
Ignition  delays  are  decreased  by  about  22%  per  10°C  increase 
in  teaperature.  The  optiaua  vapor  phase  coaposition  for 
ignition  is  found  to  be  the  equimolecular  mixture  which,  for 
'0DMH,  MMH  and  50-50  with  N204,  is  aarkedly  more  fuel  rich 
than  typical  aixture  ratios  for  optiaua  perforaance. 

A  comparison  between  experiaental  and  calculated  engine 
ignition  delays,  est lasting  for  the  calculations  the  vapor 
composition  and  teaperature,  gives  good  agreement.  To  predict 
the  vapor  coaposition  and  teaperature,  it  is  necessary  to 
Include  in  the  pressurization  analysis  expressions  for  the 
transient  N2O4  flow  upon  propellant  valve  opening  and  for  the 
heat  transfer  between  thrust  chaaber  walls  and  the  evaporating 
vapor  -  drop  systea.  Then,  further  verification  of  the  model 
is  needed. 

Five  fuel  additives  and  one  oxidizer  additive  were  tested 
for  their  effects  on  ignition  characteristics  of 
The  additives  were  tested  in  the  apparatus  used  for  the 
cheaical  kinetics  study.  This  apparatus  is  well  suited  to 
additive  screening.  The  five  fuel  additives  tested  were: 

furfuryl  alcohol 
phenyl  ether 
aethyl  butynol 
ethyl  ether 

and  benzene . 

The  additive  concentrations  were  10%  by  weight  of  the  fuel. 

Of  the  five  additives  only  furfuryl  alcohol  caused  a 
significant  effect.  The  minimum  ignition  pressure  of  NOs/^fMH 
was  reduced  by  25%  with  the  additive. 
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CO  Compound  R,  FC(NF2)3,  was  tested  as  an  oxidizer 
additive  in  amounts  of  2  and  3%  of  N02.  This  additive  not 
only  raised  minimum  ignition  pressures  but  also  made  them 
very  erratic.  It  was  necessary  to  thoroughly  wash  the 
apparatus  in  order  to  obtain  again  the  usual,  very  repeatable 
ignition  pressure  limits. 

(U)  Under  Task  III,  thrust  chamber  tests  were  made  with 
a  variety  of  thrust  chamber  conf igurations  (L*,  contraction 
ratio,  design  chamber  pressure)  to  determine  the  effect  of 
chamber  design  parameters  on  the  ignition  characteristics  of 
Compound  A  with  several  hydrazine-type  fuels  and  gaseous 
F2/H2.  The  design  chamber  pressures  selected  for  the  tests 
were  20,  75  and  200  psia,  each  for  a  nominal  thrust  level  of 
50  lbs.  Thrust  chambers  having  L*'b  of  5,  10,  30  and  50  in. 
and  nominal  contraction  ratios  of  1.5,  3.5  and  8  were 
fabricated  for  each  design  chamber  pressure  (except  for  several 
combinations  of  parameters  which  resulted  in  impractical 
configurations)  and  tested  at  a  nominal  ambient  pressure  of 
10  mm  Hg.  Transparent  thrust  chamber  cylinders  were  used  for 
visual  and  photographic  observation  of  ignition  characteristics. 

A  simple,  single-element  doublet  injector  having  minimum  mani¬ 
fold  volumes  was  used  for  most  of  the  Compound  A/hydrazine-type 
tests.  A  concentric  tube  injector  consisting  of  a  single  axial 
orifice  for  the  oxidizer  and  a  concentric  annulus  for  the  fuel 
was  used  for  the  F2/H2  tests. 

_  t 

(C)  The  thrust  chamber  tests  with  Compound  A  as  the  oxidizer 
were  made  primarily  with  N2H4  as  the  fuel.  In  addition  to  the 
investigation  of  configuration  effects,  tests  were  made  to 
determine  the  influence  of  mixture  ratio  and  propellant  leads 
on  ignition  characteristics.  Additional  tests  also  were  made 
with  UDMH  and  M9F-5  as  the  fuel  for  comparison  purposes.  MHF-5 
is  a  classified  mixed  hydrazine  fuel  consisting  of  26%  hydrazine, 
55%  monomethyl  hydrazine  and  19%  hydrazine  nitrate.  In  general 
it  was  found  that  ignition  delays  with  these  propellants  were 
very  short  and,  contrary  to  the  results  of  similar  thrust 
chamber  tests  with  ^(^/hydrazine-type  propellants  (Ref.  1), 
were  independent  of  pressure  in  the  chamber  prior  to  ignition. 
There  was  no  correlation  between  chamber  pressurization  due  to 
propellant  vaporization  (necessary  for  ignition  of  ^2°4^2^4 
propellants)  and  ignition  delay  in  these  tests. 

(C)  Ignition  delays  and  times  to  90%  of  steady  state 
chamber  pressure  were  shortest  with  Compound  A/UDMH.  Ignition 
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<C> 

occurred  approximately  7-9  msec  after  application  of  the 
electrical  signal  to  the  propellant  valves.  This  corresponds 
to  the  time  at  which  both  propellants  first  enter  the  chasber 
as  a  vapor-liquid  Mixture.  Although  Compound  A  and  UIWH 
vapors  initially  enter  the  chasber  at  3  nsec,  the  concentrations 
are  insufficient  for  ignition,  except  for  several  instances 
of  r~ak,  discontinuous,  preignition  reactions  detected  by  faint 
light  output  in  sose  tests  with  UDtfH  (and  other  fuels) . 

(C)  Ignition  vith  both  and  MHF-5  occurred  10-14  nsec 

after  the  start  signal  which  corresponds  quite  closely  vith 
the  tines  (10-12  nsec)  at  which  both  the  oxidizer  and  the  fuels 
first  enter  in  the  liquid  fttate.  MHF-5  is  first  detected 
entering  as  a  vapor-liquid  mixture  at  about  7  msec  after  start. 

In  the  ease  of  M2H4,  liquid  injection  in  the  liquid  state  is  the 
first  detectable  indication  of  fuel  entry.  M2®4  enters  as  r 
discrete,  cohesive  stream  vith  no  evidence  of  vaporization 
even  at  the  low  ambient  pressures.  In  contrast,  injected 
streams  of  Compound  A,  UDHR  and  MHF-5  spread  out  and  vaporize 
rapidly  upon  injection  even  after  steady-etate.  liquid  flowrates 
have  been  established. 

(C)  Because  of  the  lack  of  dependency  upon  pressure, 
ignition  delay  is  not  significantly  influenced  by  thrust  chamber 
configuration  parameters  (design  chasber  pressure,  L*, 
contraction  ratio,  etc.).  Although  configuration  parameters 
also  do  not  strongly  influence  the  time  to  re*ch  90%  of  steady- 
state  chamber  pressure,  there,  is  a  tre^d  toward  longer  times 
with  larger  volume  configurations  as  might  be  expected.  Mixture 
ratio  tests  indicate  slightly  shorter  ignition  delays  (1-2  usee) 
with  fuel -rich  mixture  ratios  than  oxidizer-rich  mixture  ratios. 
Propellant  lead  tests  with  Compound  A/N2H4  indicate  that  ignition 
occurs  rapidly  upon  Injection  of  (liquid)  hydrazine  regardless 
of  the  state  (vapor,  mixture,  or  liquid)  of  the  Compound  A. 

(C)  Photographic  and  pressure  instrumentation  indicate 
that  ignition  occurs  at  the  injector  with  theee  propellants 
sad  rapidly  spreads  to  the  nozzle  end.  Because  of  the  short 
ignition  delays,  the  start  transients  were  smooth  without 
evidence  of  the  random,  short  duration,  large  amplitude  spikeB 
characteristic  of  the  h^^/hydrarine-type  propellants. 

(0)  Thrust  chamber  tests  with  gaseous  also  resulted 

in  very  short  ignition  delays  and  smooth  pressure  transients. 
Ignition  generally  occurred  within  1  msec  of  injection  of  both 
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prop* 11 »ats  und  was  not  significantly  influenced  by  thruat 
ektsbflr  design  parmtcrs.  la  with  Conpound  A/hydraxin*-typ* 
propellants,  ignition  in  not  dopondont  upon  chnnbor  pressor i- 
sntlon  prior  to  ignition.  Sxcspt  for  propollnnt  load  toato  in 
which  on*  valve  was  deliberately  delayed,  ignition  occurred 
before  any  detectable  rise  is  chamber  pressure  due  to  propellant 
entry  occurred.  Ignition  occurred  at  the  injector  end  of  the 
ebaabe?  and  for  configurations  hawing  long  chanber  lengths, 
delays  up  to  about  0.8  nsec  were  measured  between  start  of 
pressure  rise  at  the  injector  and  noaal*  end*  corresponding 
to  the  epread  of  ignition  through  the  chamber. 

Use  to  90%  of  steady-stats  chamber  pressure 

were  generally  leas  than  7.5  asec,  averaging  about  5.5  nsec 
for  50  in.  L*  configurations  and  about  4.5  aeec  for  the  5  and 
10  In.  L*  configurations.  Contraction  ratio  or  design  chamber 
pressure  had  little  effect  on  the  rise  time. 

There  das  no  indication  of  pressure  spikes  with 
in  any  of  the  tests. 
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III 

TASK  I  -MEASUREMENT  OF  REACTION  BATE 


A.  Developaent  of  Mathematical  Model  of  Hypergolic  Ignition 

A  Mathematical  Model  of  hypergolic  ignition  in  attitude 
control  engines  in  space  is  desirable  for  several  reasons.  A 
Model  that  takes  into  account  the  Many  processes  involved  can 
provide  not  only  design  inforMation  for  ainiaizing  ignition 
delays,  but  also  it  ean  give  valuable  Insight  into  the  pressure 
spiking  problea  by  defining  quantitatively  the  conditions  from 
which  pressure  spikes  result,  i.e.  the  Mass  of  the  propellants 
in  the  chaaber  at  ignition,  the  fraction  of  the  propellants  in 
the  condensed  phase,  etc. 

In  the  following  sections,  a  aatheaatical  aodel  is  described 
as  well  as  the  determinations  of  soae  of  the  quantities  required 
by  the  aodel. 

The  general  approach  followed  in  developing  the  Model  is 
indicated  in  Figure  1.  Since  hypergolic  ignition  of  ^O^-type 
fuels  is  sensitive  to  pressure  (Ref.  1),  the  pressure  build-up 
in  the  chaaber  due  to  propellant  vaporization  aunt  be  deterained. 
Also  required  is  an  expression  of  the  dependency  of  hypergolic 
ignition  delay  on  pressure,  taking  into  account  aixture  ratio 
and  the  reactivity  of  the  specific  propellant  coabination.  The 
Intersection  of  the  two  curves  in  Figure  1  is  taken  to  be  the 
ignition  delay  tiae  in  the  engine  at  space  conditions.  Ignition 
delay  as  used  here  is  that  tiae  which  elapses  between  initial 
entry  of  the  propellants  into  the  thrust  chaaber  and  ignition, 
i.e.  emission  of  visible  light  accoapanied  by  an  increase  in 
pressure.  A  total  ignition  delay  tiae  for  a  given  attitude 
control  systea  is  obtained  by  adding  to  the  above  ignition  delay 
period  the  tiae  elapsed  froa  valve  signal  to  initial  propellant 
entry.  Soae  couaents  regarding  the  latter  as  well  as  the 
assumptions  inherent  in  the  ignition  aodel  indicated  in  Figure  1 
will  be  given  at  the  appropriate  points  in  the  following  sections. 

The  determination  of  the  chaaber  pressurization  curve  due 
to  propellant  vaporization  will  be  discussed  first.  Subsequently, 
the  cheaical  kinetics  aspect  of  the  aatheaatical  model  will  be 
covered  in  detail. 
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B.  Thrust  Chamber  Pressurization  Due  to  Propellant  Vaporization 

The  problem  under  consideration  in  this  section  can  be 
described  in  general  terms  as  follows.  A  propellant  stream, 
issuing  from  an  injector  orifice  into  the  low  pressure  environ¬ 
ment  of  the  thrust  chamber,  forms  a  spray  consisting  of  a  great 
many  drops  of  various  sizes.  The  drops  undergo  evaporation, 
decrease  in  size,  and  cool.  Evaporative  cooling  can  bring  a 
drop  *  s  temperature  to  its  freezing  point.  Subsequent  vaporiza¬ 
tion  occurs  at  constant  temperature  while  the  fraction  of  the 
drop  that  is  frozeu  increases  until  the  drop  becomes  completely 
frozen.  Further  evaporation  can  occur  from  the  frozen  drop 
(sublimation),  the  drop  temperature  resuming  its  decline. 

The  evaporated  gas  raises  the  pressure  in  the  chamber.  The 
pressure  at  any  time  is  influenced  by  the  chamber  geometry  and 
by  the  gas  temperature,  the  latter  in  turn  is  related  to  the 
drop  temperatures.  Some  vapors  pass  through  the  nozzle  and, 
under  certain  conditions,  vapors  can  condense  on  the  thrust 
chamber  walls. 

It  is  necessary,  therefore,  to  determine  the  following 
quantities,  all  as  a  function  of  time: 

a)  the  evaporation  from  each  drop  in  the  system 

b)  the  radius  of  each  drop 

c)  the  temperature  of  each  drop 

d)  the  fraction  of  each  drop  that  is  frozen 

e)  mass  flow  through  the  nozzle 

f)  mass  of  vapors  that  condense  on  the  chamber  walls 

g)  gas  temperature,  and  finally 

h)  gas  pressure . 

The  starting  point  of  the  analysis  is  the  work  of  Agosta  and 
Kraus  reported  in  Ref.  2.  The  present  work  was  performed  in 
close  cooperation  with  Dr.  Agosta  who  served  as  consultant.  For 
the  present  program,  their  analysis  was  modified  somewhat  and 
additional  features  were  added  to  it.  Their  basic  equations 
are  reviewed  first,  and  then  the  features  that  were  added  are 
discussed.  Finally,  comparisons  between  theoretical  and  experi¬ 
mental  results  are  given. 
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The  analysis  takes  into  account  the  rate  of  propellant 
injection,  propellant  evaporation  and  condensation  (assuming 
a  Knudsen-Langmuir  kinetic  model) ,  thrust  chamber  geometry, 
and  propellant  vapor  efflux  through  the  nozzle.  In  essence, 
the  pressure  in  the  chamber  at  any  time  is  obtained  from  a 
mass  balance  on  the  system. 

Steady  state  calculations  are  made  in  successive  time 
intervals  which  are  taken  sufficiently  short  that  all  properties 
remain  essentially  constant  during  each  time  interval.  The 
solution  proceeds  from  known  initial  conditions  to  calculated 
new  values  of  properties  at  the  end  of  the  first  time  interval. 
The  new  values  are  taken  as  initial  conditions  for  the  second 
time  interval,  etc.  to  completion. 

1.  Theory 

a.  Droplet  Evaporation 

The  heart  of  the  problem  is  droplet  evaporation. 
Based  on  the  kinetic  theory  of  gases  and  the  perfect  gas  law, 
the  amount  of  evaporation  occurring  in  a  finite  time  interval 
is  given  by  ^ 

<W  (*  fa  -  °3  fa)  -  -  Vwsr: 

o 

The  assumptions  involved  are: 

1.  The  number  of  molecules  leaving  the  liquid 
surface  per  unit  time  (when  the  ambient 
pressure  is  below  the  liquid  vapor  pressure) 
is  the  same  as  the  number  which  impinge  on 
the  same  surface  when  the  liquid  is  in 
equilibrium  with  its  vapor. 

2.  The  gas  evaporates  at  the  same  temperature 
as  the  liquid  drop  from  which  it  is  evolved. 

3.  The  molecular  weight  of  the  vapor  is  the 
same  as  the  molecular  weight  of  the  liquid, 
i.e.  a  non-dissociating  propellant. 

The  propellant  spray  is  represented  mathematically 
by  a  three  drop-size  model.  The  three  drop  sizes  are  obtained 
by  applying  a  logarithmiconormal  distribution  (with  a  specific 
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geometric  standard  deviation  depending  on  the  type  of  injector) 
to  a  given  mass  median  drop  size  and  then  selecting  three  radii 
such  that  30%  of  the  weight  flow  is  in  drops  of  the  smallest 
radius,  40%  in  the  intermediate  size  drops,  and  30%  in  the 
largest  drops  (Ref.  3).  Thus  there  are  three  classes  of  drops, 

1  <  i  <  3. 

The  continuous  processes  involved  in  the  vapor 
pressurization  of  the  thrust  chamber  are  treated  mathematically 
as  steady  state  processes  in  successive,  very  short  time  intervals. 
Consequently,  in  each  time  interval,  a  new  ensemble  of  drops 
enters  the  thrust  chambe.  These  undergo  evaporation  during 
the  time  interval  as  do  t:*<  drop  ensembles  which  entered  the 
chamber  in  previous  time  Intervals.  Therefore,  in  the  first 
time  interval  (N  -  1) ,  there  is  one  drop  ensemble  (J  »  1)  which 
consists  of  three  classes  of  drops  (1  <  i  <  3) .  In  the  second 
time  interval  (N  -  2)  there  is  not  only  the  first  ensemble 
(j  -  1)  which  continues  to  evaporate,  but  also  a  second  ensemble 
(J  -  2)  which  undergoes  evaporation  for  the  first  time.  This 
second  ensemble  also  consists  of  three  classes  of  drops. 

Therefore,  there  are  six  classes  of  drops  to  consider  in  the 
second  time  interval,  nine  in  the  third  interval,  twelve  in 
the  fourth  interval,  etc.  Each  class  of  drops  has  a  unique 
radius,  temperature  and  fraction  of  the  individual  drops  that  is 
frozen . 


Each  class  of  drops  must  be  identified  and  this 
is  done  by  i,  j  and  N  numbers  which  indicate  the  initial  radius 
of  the  drops  when  the  class  first  entered  the  thrust  chamber, 
the  time  interval  in  which  the  class  entered  the  thrust  chamber 
(or ,  in  other  words,  the  number  of  intervals  the  class  has  been 
in  the  system)  and  the  time  interval  in  question.  Thus  the 
possible  values  of  1,  j  and  N  are: 

i  -  1,  2  or  3 

j  -  1,  2,  3 - N 

N  -  1,  2,  3,  4, - . 

If  N  *  500,  say,  there  are  500  drop  ensembles  (j-numbers)  and 
each  ensemble  consists  of  3  classes  of  drops  (i-numbers) .  There¬ 
fore  there  are  1500  separate  classes  of  drops  to  be  accounted 
for  in  the  500th  time  interval,  each  class  having  a  unique  radius, 
temperature  and  quality  X  (see  below) . 
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Trosi  equation  (1) ,  the  mass  evaporated  from  a 
single  d*op  la  oa»  time  interval  is  given  by: 

£  (5u)«-  *m  a, « H  y^  (2> 

b.  Total  Evaporation  During  One  Time  Interval 

To  determine  the  tota mass  evaporated  from  all 
of  the  drops  during  one  tine  interval,  the  total  surface  area 
of  all  of  the  drops  in  the  system  must  be  known.  The  total 
surface  area  depends,  of  course,  on  the  number  and  radius  of 
the  dropa>  in  each  class. 


For  a  mass-median  drop  size,  rE,  of  73  microns 
and  a  geometric  standard  deviation,  cr  ,  of  2.3  as  given  in 
Ref.  3  for  aa  impinging  stream  injector,  the  initial  radii  of 
tie  three  drop-size  spray  model  are  determined  to  oe: 


30%  of  injectea  mass 
10“6  ft 

40%  of  injected  mass 
10“S  ft 

30%  of  injected  mass 
io~e  ft 

The  number  of  drops  in  each  of  the 
from 

where 

m  - 

73755 TTrlS'  J{  ~ 


in  drops  having  rj  -  83  5 
in  drops  having  r2  *  250  x 
in  drops  having  r3  *  720  x 
three  classes;  is  obtained 
WX  “  3  Wj 

W2  *  .4  Wj  (3) 

W3  -  .3  Wj 


and  Wj  *  w  At,  the  mass  of  propellant  injected  during  one  time 
interval. 
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Tbs  .Mwrtwr  of  drop*  tm  osefe  alas*  is  therefore : 

U%  -  13.51  x  l®10  (f3/  for  r*  -  33  s  10“6  ft 

%  -  „#1  x  10i©  for  rg  -  #80  x  10~*  ft  (4) 

s3  -  „®10  X  IS10  (wy  />!>  for  rg  -  72©  x  10"e  ft 

¥Sw  total  ’ £SH8tt»r  of  drops  injected  in  on*  iiao  interval  is 

^  »  13.14  *  #  <fj/  ^x>.  (5) 

T&«  total  s«rfa#s  sroa  of  the  drops  in  tbs  itfe  class  lsfor« 
evaporation  bogias  Is 

Ay®  4^r  rfNi«  4tTff  Nj  (ni/Nj)  (0a) 


Tbs  t^t&l  srerfa c«  area  of  all  the  drops  bsfors  evaporation  is 

A  -• 


<*tr  Mj 


tc 


♦  f. 


<©b) 


Tb»  total  evaporation  occurring  fro*  all  of  the 
drops  "in  tbs  sycts*  in  sm  ti*e  interval  is  obtained  by  expressing 
ag  of  equation  (1)  by  equations  (6b),  (5)  and  (4),  giving 


or  "  2.0®  x  1©1© 
C3  -  6.443  x  I©*® 
C2  -  .0139  x  l©10 


0 
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Pv(T4ij)N»l  is  the  vapor  pressure  corresponding  to  the  tem¬ 
perature  of  the  il£.  size  drop  of  the  j£&  ensemble  at  the  end 
of  the  N»1  time  interval.  The  density  of  the  liquid,  is 

assumed  constant. 


c.  Vapor  Condensation  on  Chamber  Walls 


The  aass  of  vapors  condensing  on  the  chaaber 
vails  is  given  by 


* 


Py 


(Ty)  «  At 


at 


N 


(8) 


but  with  the  constraint  (KG*.  <  +00.  Condensation  on  the  vail 
is  considered  only  when  Pg(Tg)  >  Py(Tv);  that  is,  the  conden¬ 
sation  tera  cannot  have  a  negative  value  thereby  implying 
evaporation. 


d .  Vapor  Flow  Through  Nozzle 

The  mass  efflux  through  the  nozzle  during  one  time 
interval,  for  constant  c  and  k  for  the  gas,  is  given  by 


e.  Vapor  Mass  in  Chamber 

The  xass  of  vapor  in  the  chamt  r  at  the  end  of 
the  JUdEl  time  i^ter'  al  is  obtained  from  the  perfect  gas  law 


*•  Ss&Jssm***** 

The  temperature  of  the  gas  in  the  chamber  at  the 
end  of  the  Kill  time  interval  is  obtained  by  taking  a  mass  weighted 
average  of  the  temperature,  which  includes  the  temperature  of 
the  gas  generated  du~  to  evaporation  during  the  time  interval 
plus  the  temperature  of  the  gas  left  from  the  previous  time 
interval  minus  the  temperature  of  the  gas  which  condenses  on 
the  chamber  walls  and  that  which  passes  through  the  nozzle. 
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for  n  o  ingle  specie  systen  aid  a  constant  vapor  host  capacity, 

.  jL  fa  fcteW  H-.  Jh j  w-T  (ftv«)«  ~  ("""»»)  *1  (ll) 

*"j«  ^  &U  *  (&w)m  " 

* •  Qan  Pregsnre 

A  nans  balance  of  the  ay a ten  yields  the  gas 
pressure  in  the  ohanber  at  the  eng  of  the  Kill  tine  interval. 

(H,  *u)n  ■  W«  '  =  K)  »- K)«.,  <i« 


The  terns  on  the  left  hand  side  are  giyen  by  agnations  (7) ,  (8) , 
and  (9)  chile  the  right  hand  side  is  equation  (10)  evaluated  at 
the  end  of  the  V  and  3b*  1  tine  Intervals. 

The  gas  pressure  in  the  chanber  at  the  end  of  the 
B&A  tine  Interval,  <pg>i»  im  *iy»®  by; 


«**<*..  />irR(T8jM„ 
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At  the  cz'J  of  every  tiae  interval,  new  drop  radii 
and  drop  teaperaturea  aust  be  calculated  to  account  for  the 
effect  of  evaporation  during  the  tiae  interval  on  the  drop's 
else  and  teaperature.  A  aass  balance  on  the  drop  together  vith 
equation  (2)  yields  the  new  drop  radius: 


1.  Drop  Teaperature  and  "Qualities' 


The  new  drop  teaperature  is  obtained  froa  an 
energy  balance  on  the  drop. 


Clearly  the  teaperature  within  the  drop  is  assumed 
to  be  uniform.  The  heat  of  vaporization,  A*,  the  b«»t  capacity 
of  the  liquid,  cpl,  and  the  liquid  density,  />i,  are  assuaed 
constant.  Values  of  these  quantities  corresponding  to  the  aid- 
teaperature  of  the  range  encountered  are  used. 

At  soae  point,  an  evaporatlvely  cooling  drop  can 
arrive  at  its  freezing  point.  Further  evaporation  of  th*  drop 
occurs  at  constant  teaperature  during  which  tine  the  fraction  of 
the  drop  that  is  frozen  increases.  Thus,  equation  (15)  holds 
until  the  freezing  point  is  reached.  Once  T<j.  .  equals  Typ , 
than  1*^  4  resains  constant  as  the  drop  freezes*?  The  fraction  of 
th*  dropJthat  freenes  in  one  tiae  internal  due  to  evaporation 
is  given  by  an  energy  balance  on  the  drof. 
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Th©  solid  fraction  of  the  drop,  that  is  its  "quality ",  increases 
in  successive  intervals  until  the  drop  is  frozen  solid.  The 
quality,  X*j,  of  a  drop  is  given  by 


Xu -I  * 

H.Nre 


* fit  Kij 


r 

N 


(17) 


with  the  constraint  that  0  <  <  1. 

The  summation  over  successive  tine  intervals  for  the  drop  in 
question  begins  when  the  tenperature  of  the  drop,  T^ja,  reaches 
Tpp.  The  summation  ends  when  the  drop  is  completely  frozen, 
i.e.  Xij  -  1.  Gij  in  equation  (17)  is  given  by  equation  (2). 

Once  the  drop  is  frozen  solid,  further  evaporation 
(sublimation)  causes  the  temperature  of  the  drop  to  decrease. 
Equation  (15),  modified  to  reflect  the  solid  state  of  the  drop¬ 
let,  becomes  applicable  again. 

The  above  equations  were  programmed  for  a  Control 
Data  G-20  computer.  The  time  Intervals,  Atg,  used  in  the 
calculations  ranged  from  1  to  40  microseconds  depending  on  the 
particular  propellant  and  motor  geometry  under  consideration. 

In  general,  the  faster  the  pressure  rise  in  the  thrust  chamber, 
the  smaller  the  time  interval  required.  The  actual  criterion 
for  the  length  of  the  time  interval  is  that  the  change  in 
properties  during  each  time  interval,  i.e.  drop  radius,  drop 
temperature,  etc.,  be  kept  small  enough  so  that  the  properties 
can  be  considered  constant  during  the  tine  interval.  New  values 
are  calculated  at  the  end  of  each  time  interval  and  these  are 
used  for  the  next  interval,  etc.  to  completion  of  the  calculation. 

The  equations  are  used  to  calculate  in  each 
time  interval  the  gas  pressure,  the  gas  temperature,  and  a  new 
radius  and  tenperature  for  each  drop  in  the  system.  Since 
the  injected  propellant  stream  is  represented  by  three  drop 
sizes,  the  number  of  classes  of  drops  increases  by  three  in 
each  successive  time  interval.  Each  class  of  drops  has  a 
unique  radius,  temperature  and  drop  "quality"  X.  In  the  lOOdUl 
time  Interval,  say,  there  are  300  drop  radii,  300  drop  tempera¬ 
tures,  and  300  drop  "qualities"  to  be  accounted  for.  The 
changes  in  the  values  of  each  of  these  must  be  calculated  to 
obtain  the  gas  pressure  at  the  end  of  th?  101®!  time  interval. 
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2.  Result a  and  Discussion 


The  pr seeding  equations  were  solved  to  determine  the 
pressure  history  in  a  thrust  chamber  due  to  vaporisation  of 
carbon  tetrachloride.  CC14  was  selected  because  the  many 
physical  properties,  especially  accommodation  coefficient, 
required  in  the  computer  program  are  known  for  this  substance. 
Corresponding  properties  for  modern  propellants  such  as  N0O4 
and  Compound  A  are  often  unknown  at  present  or  of  obviously 
uncertain  precision.  Thus,  the  model  for  chamber  pressurization 
due  to  propellant  vaporization  was  developed  and  checked 
experimentally  using  CCI4  as  the  working  fluid.  The  vapor 
pressure  of  CCI4  is  very  similar  to  that  of  the  50-50  blend  of 
N2H4-UDMH  for  the  range  of  temperatures  encountered. 

a.  Gas  Temperature  History 

The  initial  calculations  showed  that  the  vapor 
temperature  in  the  thrust  chamber  dec  eases  initially  as  the 
chamber  pressure  increases  but  after  few  milliseconds  (de¬ 
pending  on  motor  geometry  and  the  propellant  properties)  the 
vapor  temperature  reaches  a  minimum  and  thereafter  gradually 
rises  until  steady-state  conditions  are  reached  (Figure  2) . 

The  initial  cooling  is  the  result  of  evaporation 
from  drops  which  have  already  undergone  some  evaporative  cooling. 
The  minimum  in  the  gas  temperature  history  occurs  because  the 
evaporation  from  the  coldest  drops  is  more  strongly  diminished 
by  the  rising  chamber  pressure  than  the  evaporation  from  t..e 
warmer  (newer)  drops.  A  point  is  reached  at  which  the  net 
production  of  cold  vapors  equals  the  net  production  of  warm 
vapors  and  no  the  gas  temperature  remains  unchanged. 

The  gas  temperature  subsequently  increases  gradually 
as  a  consequence  of  the  rising  gas  pressure  which  first  diminishes, 
than  steps,  and  finally  causes  "negative  evaporation"  v>f 
colder  drops.  The  temperature  inversion  is  of  interest  because 
it  occurs  generally  during  typical  ignition  delay  times.  The 
cooler  gas  temperatures  would  tend  to  lengthen  ignition  delay* 
due  to  the  effect  of  temperature  on  reaction  rates. 

»,  Droplet  Condensation 

Condensation  of  the  vapors  was  found  to  occur  on 
drops  which  have  undergone  substantial  evaporative  cooling.  The 
condensation  results  when  the  gas  pressure  in  the  chamber  exceeds 
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the  vapor  pressure  of  the  cold  drop,  see  equation  (1)  or  (2). 
An  illustration  of  condensation  on  a  cold  drop  is  given  in 
Table  I. 


TABLE  1 

EXAMPLE  OF  CONDENSATION  ON  A  COLD  DROP  (CCla  #R) 


Time 

Drop  Radius 

Drop  Temp. 

Drop  V.P . 

Chamber  Pressure 

msec 

x  lO-o  ft 

oft 

mmHg 

mmHg 

2.80 

674 . 907 

450.00* 

8.291 

6.804 

3.20 

674 . 550 

450.00* 

8.291 

7.680 

3.60 

674.495 

450.02 

8.298 

8.484 

674 . 674 

450.44 

8.416 

9.224 

4.40 

675.051 

451.31 

8.671 

9.913 

♦Drop  is  partially  frozen. 


The  drop  under  consideration  in  Table  I  had  a 
radius  of  720  x  10~6  ft  when  it  entered  the  system  which  was  at 
time  zero.  It  is  one  of  the  i  -  3  class  of  drops,  the  largest 
drops  of  the  three  drop-size  spray  model.  After  2.80  msec,  the 
drop  radius  has  decreased  to  the  value  shown  and  its  temperature 
has  reached  the  freezing  point.  The  drop  is  not  frozen  solid  as 
yet  but  is  only  partially  frozen.  At  this  time  the  chamber 
pressure  is  still  below  the  vapor  pressure  of  the  drop.  During 
the  next  0.40  msec,  the  drop  continues  to  undergo  evaporation 
(radius  decreases  further)  since  the  gas  pressure  (column  5)  is 
still  less  than  the  drop's  vapor  pressure  (column  4).  During 
this  time  the  fraction  of  the  drop  that  is  frozen  increases, 
the  drop  temperature  remaining  at  the  freezing  point  value. 

At  3.60  msec,  the  chamber  pressure  has  increased 
to  a  value  greater  than  the  drop's  vapor  pressure  thereby 
causing  condensation  onto  the  drop.  The  drop's  temperature  has 
increased  to  450.02°F  indicating  that  the  condensation  energy 
absorbed  by  th©  drop  not  only  melted  the  frozen  fraction  of  the 
drop  but  also  raised  slightly  the  temperature  of  the  all-liquid 
drop.  Although  the  condensation  increased  the  mass  of  the  drop, 
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the  Nt  eff  et  ow  the  0.40  asec  period  la  a  decreased  radius 
due  to  the  evaporation  which  occurred  during  the  first  part  of 
the  period,  l.e.  until  the  gas  pressure  equalled  the  drop’s 
vapor  pressure .  During  the  next  0.40  ease  period,  that  la  at 
4.00  tsaec,  condensation  on  the  drop  continued  (collier.  5  vs. 
colons!  4).  The  temperature  of  the  liquid  drop  increased  as  did 
the  sis#  of  the  drop,  far tear  condensation  occurred  during  the 
last  iaterval  listed  im  the  table  as  the  gae  preeeure  continued 
to  Ini  greater  than  the  vapor  pressure  of  the  drop . 

c.  yapor-Qrofr  Knur  ary  Transfer 

The  /'ossibilitj'  of  heat  transfer  between  the 
rpor  and  the  drops  la  not  accounted  for  by  the  equations 
given  previously.  Ir.  the  case  of  CCl*,  the  teapersture  differ- 
soce  etweeo  thj  wars  vapors  and  the  coldebt  drops  can  be  as 
sr-i.  as  75°F.  With  such  a  driving  force,  one  would  expect  the 
drops  to  be  waraed  and  the  *apar*  cooled  to  soae  extent . 

The  energy  reaching  the  surface  of  s  drop  during 
en#  ti*^  interval  is.  according  tc  ref.  3: 


m  m  H  [r H  At  N  (^4  t  i)H^ 


(18) 


where  2 


a  i  * 


^  * 


4Tf  h  4t:N  (ru)f:^ 
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Hr  ** 


and  Re 


_  ^  -t  :> 


**-? 


excluding 
drag  effects. 


(21) 


(22) 


Gij  la  equation  (IS)  is  give.;  *^y  <?q*.ation  (2) 
An  energy  balance  on  the  d~ op  is  then 
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fro*  which  the  new  drop  temperature  at  the  end  of  the  ti*e 
interval  in 
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and  |»  ,  z  and  z,  >•  Pr,  &e  and  Q;ij  sr*  given  br  aquations 
(25) ,  (18),  (20),  (£2?  and  (2)  respectively. 


Twa  total  energy  reaching  the  surface  of  all 
drop*  in  oyatess  in  on#  tine  interval  1#; 
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where  Rj ,  the  total  number  of  drop*  in  oash  is  given 

by  equation  (4)  &cd  b  and  1  tre  by  equations  <20}  and  (10) 

respectively.  Tfce  eneryv  Q$  coaes  from  the  gas  in  taa  chamber 
and  therefore  equation  (11).  the  gap  temperature  equation, 
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The  effect  of  including  vapor-droplet  beet 
transfer  on  chamber  pressure  and  temperature  is  alight  %m  in¬ 
dicated  in  the  first  two  lines  of  Table  II. 


As  expected,  the  gas  temperature  is  lower  wh*,a 
heat  transfer  t©  the  cold  drops  i*  included.  Tfc®  change  is  small 
however.  The  higher  gas  pressure  results  fro®  the  drop®  -*hich 
are  slightly  warmer  due  to  the  heat  transfer  and  so  the  drops 
vaporise  slightly  more  s&s®  in  succeeding  time  intervals, 
increased  mass  cf  vapor  more  than  offsets  the  gas  temperature 
effect  on  chtther  pressure  *  -  ■*%  the  net  change  in  the  latter 
is  small. 

:•! ,  Drop  hesidence  Tisse 

T»*e  equations  discussed  to  this  point  put  no 
limit  on  the  lifetise  of  the  drops  «n  the  system.  2b  reality, 
the  lifetime  of  the  drops  is  finite,  the  drops  pass  through 
the  no*»I«  after  some  sort  of  residence  tis#.  Therefore,  a 
residence  time  limitation  was  included  in  the  program-  The 
residence  tine  used  is  defined  simply  as  the  average  time  re¬ 
quired  for  a  drop  to  travel  fro®  the ' injector  to  a  chamber  wail 
travei  *.®g  at  a  coaster*  velocity  -equal  to  the  iajsction  velocity . 
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Despite  the  short  residence  tines  to  which  the  model  leads, 

1-1/2  to  5-1/2  nsec  typically  depending  on  notor  geometry  and 
injection  velocity,  the  effect  on  chamber  pressure  is  very 
small.  A  comparison  of  the  second  and  third  lines  of  Table  II 
shows  the  extent  of  the  effect.  An  arbitrary  reduction  of  the 
residence  time  from  5  to  3.36  msec  (fourth  line  of  Table  II) 
shows  again  the  smallresft  sf  the  residence  time  effect. 

Drop  temperature  histories  show  that  once  a  drop 
has  been  in  the  system  for  one  residence  time  period,  its 
temperature  is  such  that  its  vapor  pressure  is  not  very  different 
fruM  tne  gas  pressure  prevailing  at  that  time.  Thus,  in  accor¬ 
dance  with  equation  (1)  or  (2),  the  amount  of  evaporation  (or 
condensation)  from  that  point  onward  is  relatively  small. 

Table  II  shows  that  the  residence  time  limitation 
causes  a  slightly  higher  chamber  pressure  and  a  slightly  lower 
gas  temperature.  These  effects  result  because  some  of  the 
vapor  condensation  onto  the  cold  drops  is  stopped.  Thus,  some¬ 
what  more  gas  at  low  temperatures  is  in  the  chamber  than  would 
otherwise  be  the  case. 

e.  Effect  of  Drop  Sizes 

The  initial  radii  of  injected  drops  for  the 
three  drop-size  spray  model  are  derived  from  a  given  mass- 
median  drop  size  as  outlined  previously.  The  mass-median  drop 
side  used  initially  in  the  calculations  was  75  microns,  based 
on  the  work  of  Priem  and  Heldman  (Ref.  3)  who  cite  investiga¬ 
tions  of  the  spray  characteristics  of  liquid  streams  emerging 
into  an  environment  in  which  the  pressure  was  substantially 
above  the  vapor  pressure  of  the  liquid.  At  ambient  pressures 
well  below  the  vapor  pressure  of  the  liquid,  one  would  expect 
additional  dispersive  forces  to  be  operative  resulting  in 
smaller  drops  on  the  average. 

A  number  of  short  duration  (approx.  10  micro¬ 
second)  exposures,  with  a  magnification  of  2X,  was  taken  of 
streams  of  various  propellants  issuing  from  an  injector  orifice 
into  a  low  pressure  environment.  Exposures  were  made  at  various 
pressures  from  one  atmosphere  to  1  mi  Hg,  using  a  General  Radio 
Company  High  Speed  Stroboscope,  Type  No.  1533- A,  in  the  single 
pulse  node.  A  Bausch  fc  Lonb  7-1/8"  coated  F/6.8  Proto  Vila  lens 
permitted  magnification  of  the  object.  The  exposures  were 
recorded  os'1  Polaroid ,  Type  57,  300  speed,  4x5  film.  Hie 
streams  are  back-lighted,  a  diffuser  being  used  between  the 
stream  and  the  light  source.  The  direction  of  stream  flow  was 
vertically  upward. 
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A  strong  offset  of  ambient  pressure  on  tho 
sprays  resulting  fro*  CCI4  liquid  Injection  Is  shown  in  the 
pictures  of  Figure  3.  At  high  ambient  pressures,  100  ran  Hg 
and  above,  th*  relatively  low  velocity  stream  is  coherent 
for  the  two-inch  stream  length  viewed.  At  95  ran  Hg,  which 
is  the  vapor  pressure  of  CCI4  at  71°F,  a  vapor  bubble  in  the 
stress  appears  to  have  “exploded".  At  somewhat  lower  ambient 
pressures,  the  "bubble  explosions"  are  more  numerous  and  liquid 
ligaments  appear.  At  still  lower  pressures,  the  ligaments 
are  smaller  as  are  the  drops  on  the  average.  At  10  and  1  mm  Hg, 
the  ligaments  are  short  lived  and  many  of  the  droplets  approach 
the  limit  of  resolution  in  the  pictures  which  is  about  50  microns 
in  diameter. 


Similar  pictures  for  UDMH  sprays  at  various 
ambient  pressures  are  given  in  Figure  4.  The  liquid  injection 
velocity  of  the  streams  in  Figure  4  for  UDMH  is  nearly  the 
same  as  in  Figure  3  for  CCI4. 

Obvious  similarities  exist  between  the  UDMH 
streams  and  the  CCI4  streams.  In  both  cases  the  streams  are 
coherent  at  the  high  ambient  pressures  and  "exploding  bubbles” 
first  appear  when  the  ambient  pressure  is  only  slightly  below 
the"  vapor  pressure  of  the  liquid  (vapor  pressure  of  UDMH  at 
"l0?  is  130  mm  Hg) .  .it  lower  pressures,  slight  differences  are 
acted  buy  these  differences  are  in  degree  only.  UDMH  appears 
to  form  liquid  ligaments  wore  readily  than  CCI4;  and,  on  the 
average,  the  UDMH  droplets  at  the  lowest  ambient  pressures  are 
somewhat  smaller  than  the  corresponding  CCI4  droplets.  This 
is  due,  at  least  partially,  to  the  greater  evaporation  rate  of 
UDMH  which  results  fro®  its  somewhat  greater  volatility. 

The  effect  on  spray  characteristics  of  higher 
flowrates  of  UDMI?  through  the  aarae  injector  orifice  is  shown 
in  Figure  5  for  an  ambient  pressure  of  140  mm  Hg  (slightly 
above  tbs  vapor  pressure  of  UDMH)  and  in  Figure  6  for  an  ambient 
pressure  o#  i  mm  Hg.  At  the  higher  ambient  pressure  the  higher 
flowrates,  which  give  a  higher  injection  velocity,  lead  to 
greater  break-up  of  the  stream.  At  the  low  ambient  pressure, 
however,  any  differences  between  the  sprays  from  the  various 
flowrates  are  much  less  obvious.  The  degr«e  of  streaking  in 
the  pictures  by  tW  droplets  increases  with  injection  velocity 
but  the  drop  sixes  (width  of  streak)  appear  roughly  similar 
in  all  four  pictures.  Tho  distance  to  essentially  complete 
stream  break-up  is  affected  only  slightly  by  the  injection 
velocity,  or  mass  flowrate.  Assuming  the  drops  travel  a*'  a 
coustant  velocity  equal  to  the  injection  velocity ,  break-up 
or.  the  streams  is  complete  in  less  than  a  millisecond. 
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for  comparative  purposes,  similar  spray  pictures 
of  V2O4  are  reproduced  from  Kef.  1  iu  figures  7  and  8.  The  mass 
flowrate  is  the  seme  in  all  six  pictures.  The  different 
injection  velocities  in  the  two  figures  are  obtained  by  different 
injector  orifice  diameter 3.  Again,  the  strong  effect  of  ambient 
pressure  is  obvious.  Since  flowrate  is  constant  in  Figures  7 
and  8,  the  effect  of  injection  velocity  alone  on  stream  break¬ 
up  tine  Can  be  discerned.  In  the  high  velocity  case,  stream 
break-up  is  complete  in  1/3  msec  at  an  ambient  pressure  of 
1$6  mm  Hg.  In  the  low  velocity  case,  the  stream  break-up  time 
is  about  1.5  msec  at  40  an  Eg  ambient  pressure. 

A  complete,  quantitative  drop  size  distribution 
cannot  be  determined  from  the  photographs;  nevertheless  semi- 
quantitative  Information  can  be  gleaned  by  counting  the  re¬ 
solved  drops  and  measuring  their  diameters  at  successive  sec¬ 
tions  of  the  spray. 

Analysis  of  the  pictures  of  the  low  pressure 
CCI4  sprays  indicates  that  the  actual  mass-median  drop  size 
is  less  than  the  75  micron  radius  initially  assumed  for  the 
sprays  in  the  computer  program.  How  much  less  could  not  be 
determined.  Nevertheless,  the  effect  of  drop  size  on  pressuri¬ 
zation  of  a  thrust  chamber  was  determined  by  a  calculation  using 
smaller  drop  sizes.  The  mass-median  drop  size  was  arbitrarily 
reduced  from  the  original  75  micron  radius  to  a  50  micron 
radius.  New  rauii  for  the  three  drop-size  spray  model  used  in 
the  computer  program  were  computed,  again  using  the  logarith- 
miccnormal  distribution,  the  same  ge  metric  standard  deviation, 
and  again  selecting  three  radii  such  that  30%  of  the  weight 
flow  exists  in  drops  of  the  smallest  radius,  10%  of  the  weight 
flow  in  drop®  of  the  Intermediate  *  idiua  and  the  remaining  30% 
in  drops  of  the  largest  r*.Jiu@.  The  resulting  radii  are  given 
in  Table  III  together  with  the  radii  derived  from  the  75  micron 
mass-median  drop  size  used  initially. 

Computer  calculations  were  made  using  each  of 
the  sets  of  radii  in  Table  III  to  represent  the  spray.  The 
calculated  pressure  curves  for  each  case  are  given  in  Figure  f 
As  expected,  the  spray  consisting  of  the  smaller  drops  yields 
the  higher  pressures  due  to  the  increased  Surface  area  re¬ 
sulting  from  the  greater  number  of  drops.  A  compensating 
effect  on  chamber  pressure,  however,  is  that  due  to  t.  ?  higher 
mass  fraction  initially  evaporated  fro*,  the  smaller  drops, 
their  temperature  is  lower  is  succeeding  time  intervals.  The 
lower  temperature  causes  a  relatively  lesser  amount  oi  evapora 
t ion  in  succeeding  interval®.  &&  a  result,  the  surface  area 
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For  comparative  purposes,  similar  spray  pictures 
of  HgO^  are  reproduced  from  Eef.  1  in  Figures  7  and  8.  The  mass 
flowrate  is  the  same  in  all  sis  pictures.  The  different 
injection  velocities  in  the  two  figures  are  obtained  by  different 
injector  orifice  diameters.  Again,  the  strong  effect  of  ambient 
pressure  is  obvious.  Since  flowrate  is  constant  in  Figures  7 
and  8,  the  effect  of  injection  velocity  alone  on  stream  break¬ 
up  time  can  be  discerned.  In  the  high  velocity  case,  stream 
break-up  is  complete  in  1/3  seec  at  an  ambient  pressure  of 
100  mm  Hg.  In  the  low  velocity  case,  the  stream  break-up  time 
is  about  1.5  msec  at  40  tarn  Hg  ambient  pressure. 

A  complete,  quantitative  drop  size  distribution 
cannot  be  determined  from  the  photographs;  nevertheless  semi- 
quantitative  information  can  be  gleaned  by  counting  the  re¬ 
solved  drops  and  measuring  their  diameters  at  successive  sec¬ 
tions  of  the  spray. 

Analysis . of  the  pictures  of  the  low  pressure 
CCI4  sprays  Indicates  that  the  actual  mass-median  drop  size 
is  less  than  the  75  micron  radius  Initially  assumed  for  the 
sprays  in  the  computer  program.  How  much  less  could  not  be 
determined.  Nevertheless,  the  effect  of  drop  size  on  pressuri¬ 
zation  of  a  thrust  chamber  was  determined  by  a  calculation  using 
smaller  drop  sizes.  The  mass-median  drop  size  was  arbitrarily 
reduced  from  the  original  75  micron  radius  to  a  50  micron 
radius.  New  radii  for  the  three  drop-size  spray  model  used  in 
the  computer  program  were  computed,  again  using  the  logarlth- 
miconormal  distribution,  the  same  geometric  standard  deviation, 
and  again  selecting  three  radii  such  that  30%  of  the  weight 
flow  exists  in  drops  of  the  smallest  radius,  40%  of  the  weight 
flow  in  drops  of  the  intermediate  radius  and  the  remaining  30% 
in  drops  of  the  largest  radius.  The  resuittr-^  radii  are  given 
in  Table  III  together  with  radii  derived  from  the  75  micron 
.js  a.crJj.ftn  drop  size  used  initially. 

Computer  calculations  were  made  using  each  of 
the  sets  of  radii  in  Table  III  to  represent  the  spray.  The 
calculated  pressure  curves  for  each  case  are  given  in  Figure  9. 

As  expected,  the  spray  consisting  of  the  smaller  drops  yields 
the  higher  pressures  due  to  the  increased  surface  area  re¬ 
sulting  from  the  greater  number  of  drops.  A  compensating 
effect  on  chamber  pressure,  however,  is  that  due  to  the  higher 
mass  fraction  initially  evaporated  from  the  smaller  drops, 
their  temperature  is  lower  in  succeeding  time  intervals.  The 
lower  temperature  causes  a  relatively  lesser  amount  of  evapora¬ 
tion  in  succeeding  intervals.  As  a  result,  the  surface  area 
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Velocity, 
ft /sec 


Flowrate 

lb/sec 


0.0055 


Figure  5.  Fully  Developed  UDSEH  Flow  at  Various  Liquid 

Injection  Velocities.  Ambient  Pressure  «  140  mm  Eg, 
Orifice  -  0.021  in.  dia .  z  0.042  in.  long, 

Ambient  Temperature,  Magnification  2X. 
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Figure  6.  Fully  Developed  UDMH  Flow  at  Various  Liquid 
Injection  Velocities.  Ambient  Pressure  -  1 
Orifice  -  0.021  in.  dia.  x  0.042  in.  long, 
Ambient  Temperature,  Magnification  2X. 


Flowrate 

lte/asc 


0.0055 


0.0083 


0.0110 


0.0138 


■«  Hg, 
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(b)  150  sun  Hg 


Figur©  7.  Fully  Developed  H2O4  Flow  at  Various 

Ambient  Pressures.  Liouid  Injection  Velocity  - 
110  ft/sec,  Flowrate  -  0.030  Ib/sec,  Orifice  - 
0.024  in.  dla.  x  0.048  in.  long,  Ambient  Tempera¬ 
ture,  Magnification  1.8X. 
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(c)  40  mm  Kg 

Figure  8 „  Fu^y  Developed  N204  Flow  *t  Various  Ambient 

Pressures.  Liquid  Injection  Velocity  -  11  ft/sec 
Flowrate  -  0.030  Ib/sec,  Orifice  -  0.073  in.  dia. 
x  0.146  in.  long,  Ambient  Temperature, 

Maomif ication  1.8X. 


Pres*sv 


A  FRF  I 


x  n  c  r  €■  &  s  fr 

somewhat 


due  to  the  s^lhr  drop  sizes  i fs  part  tally  offset  by 
lower  drop  temperatures.  The  net  effect  of  the 


smaller  drop  sizes  on  chamber  pressurization  is  rather  small . 


TABLE  III 


DROP  SIZES  USED  IN  THE  COMPUTER  PROGRAM 
TO  ZEPRESEVT  TjJE  INJECTED  PROFELLANT  UP  RAY 


Mass-Median 

Radii 

Drop  Size 

i  »  1 

i  «  2 

* 

- 

microns 

micro  feet 

75 

83 

250 

720 

50 

58.3 

171 

384 

Figure  9  shows  that  the  absolute  pressure 
difference  between  the  two  curves  increases  with  pressure. 
However,  the  percentage  pressure  difrerence  is  found  to  be 
nearly  constant  bevond  about  3  msec.  Th#.  constant  «-<•«*»»  - 
tage  increase  in  pressure  is  about  7%  for  the  30%  reduction 
'n  mass-median  drop  size.  Due  tc  the  uncertainties  in  actual 
drop  sizes  for  propellant  streams  issuing  into  a, low  pressure 
environment,  it  is  comforting  to  find  that  the  influence  of 
spray  drop  sizes  on  chamber  pressurization  is  quite  small. 

f .  Effect  of  Accommodation  Coef f icient 

Oversimplified,  the  accommodation  coefficient 
may  be  considered  as  a  measure  of  the  ease  with  which  a  molecule 
of  a  substance  evaporates  from  or  condenses  on  a  liquid  surface 
of  that  specie.  The  accommodation  coefficient  is  unknown 
generally  for  modern  hvpergolic  propellants  as  well  as  for 
©any  common  liquids.  Furthermore,  it  is  a  property  extremely 
difficult  to  measure. 

To  determine  the  importance  of  accommodation 
coefficient  on  chamber  pressurization,  calculations  were  made 
with  CCI4  as  propellant  but  with  two  values  for  its  accom¬ 
modation  coefficient,  c*  .  The  recognized  value  of  unity  for 
CCI4  was  used  in  one  calculation  and  then  it  was  arbitrarily 
reduced  to  one-half  for  a  second  calculation.  All  other 
inputs  were  kept  the  same  in  the  two  computer  runs.  The 
resulting  chamber  pressure  curves  are  given  in  Figure  10. 
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Figure  10.  Effect  of  Accomodation  Coefficient 
on  Pressurization  Rate  -  CCl* 
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4k  in  the  case  of  the  drop  sizes,  the  effect 
of  accommodation  coefficient  on  chamber  pressurization  is 
quite  email.  4gain,  the  absolute  pressure  difference  increases 
with  pressur  ;  but,  here,  the  percentage  pressure  difference 
gradually  decreases  as  the  chamber  pressure  increases.  At  2 
msec,  the  pressure  corresponding  to  the  lower  accommodation 
coefficient  (  be  -  0.5)  is  almost  17%  below  the  pressure  for 
the  <x  -  1.0  case.  The  difference  dirainisnes  to  about  11% 
at  4.00  msec,  and  to  less  than  10%  <*t  6,00  msec  for  the  50% 
reduction  in  <x  . 


Initially,  the  amount  of  vaporization  of  a 
given  size  drop  when  «  -  1.0  is  twice  that  for  »  1/2. 

This  factor  is  not  maintained,  however,  since  the  drop  with  the 
lo~er  e<  is  cooled  to  a  lesser  extent  due  to  the  lesser  amount 
of  evaporation.  Thus,  in  the  second  time  hherval,  the  drop 
with  ©;  -  1/2  is  warmer  than  the  drop  with  <*  -  1.0  and 
therefore  (by  equation  1  or  2)  the  mass  evapor«ted  in  the  aecczd 
time-interval  by  the  drop  with  «  -  1/2  is  more  thau  one-half 
the  mass  e  aporated  by  the  C*  «  1.0  drop.  Consequently ,  the 
direct  effect  of  a  low  accommodation  coefficient  is  partially 
offset  by  th-*  consequent  slower  drop-cooling  rate.  The  net  effect 
on  chamber  pressurization  Ip  much  lees  tnan  linear.  Again,  due 
to  >.he  considerable  uncertainties  in  octuai  accommodation 
coefficients,  it  is  fortunate  tha*  its  effect  un  chamber  pres¬ 
surization  is  rather  small. 

g.  Thrust  Chamber  Wall  Effects 

An  experimental  chamber  pressurization  curve  was 
measured  using  high  response  Kistier  piezoelectric  pressure 
transducers.  Tne  experimental  conditions  are  as  follows  (CCl^ 

Test  *16-19}  : 

prope 1 1 - n  t 

Design  P^,  peia 

A* ,  aq ,  in. 

L*,  inch  10.3 

Contraction  ratio  9.0 

4 set  doublet 


200 
0. 137 


Injector 
Flow  passage 


Outside  holes 
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Orifice  diameter  (4),  inch 
Flowrate ,  lbm/sec 
Flow  control 

< 

Liquid  injection  velocity, 
ft/sec 

Injector  volume,  cu.  in. 
Initial  system  temperature ,°R 
Initial  ambient  pressure,  rnwHg 


0.025 

0.0816 

Cavitating  venturi 

60.7 

0.0077 

530 

1 


Carbon  tetrachloride  was  used  to  check  experi- 
mentally  the  theoretical  chamber  pressurization  model  because 
of  its  known  accommodation  coefficient  and  other  properties 
required  in  the  computer  program.  The  CCI4  mass  flowrate  is 
approximately  one-half  the  total  mass  flow  of  NoO^/UDMH  which 
generates  50  lbs.  thrust  in  the  same  engine.  Fast-acting 
electric  venturi  valves,  which  attach  directly  to  the  injector, 
were  used  (see  Section  VA2) . 

The  experimental  chamber  pressure  curve  for  CC14 
is  given  in  Figure  11  together  with  several  calculated  pressure4 
curves  based  on  the  same  engine  and  CCI4  flow.  The  times  for 
the  experimental  curve  are  measured  from  the  first  indication 
on  the  oscilloscope  record  of  a  pressure  rise  iu  the  thrust 
chamber.  The  time  from  valve  signal  to  propellant  entry  is 
obviously  not  included. 

The  drop  sizes  used  in  the  calculations  are 
those  given  in  Table  III  for  a  mass-median  drop  size  of  50 
microns.  The  residence  time  limitation  for  the  lifetime  of 
the  drops  in  the  system  is  included  in  the  computer  runs  as  is 
heat  transfer  between  the  vapor  and  the  drops.  Also,  the 
appropriate  accommodation  coefficient  (unity)  is  used.  The 
calculated  pressure  curve  under  these  conditions  is  given  in 
Figure  11a;  it  is  labeled  "no  wall  heating".  For  the  first 
1-1/2  msec,  agreement  between  the  calculated  and  experimental 
curves  is  considered  satisfactory;  beyond  that  time  agreement 
is  clearly  inadequate. 
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Flgr  »  11a.  Experimental  and  Theoretical  Chamber 

Pressure  Histories  -  CCI4 
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Figure  lib.  Bxperleental  and  Theoretical  Clumber 

Pressure  Histories  -  CCI4 
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Ac.n'ji ■ 'ing  to  the  definition  of  residence  time 
given  previously  (ta.  average  time  for  a  drop  traveling  at  the 
injection  velocity  to  reach  the  chamber  wall),  the  residence 
time  for  the  case  under  consideration  is  1.4  msec.  Drops  that 
have  been  in  the  system  for  this  length  of  time  have  cooled 
substantially;  in  the  present  case,  after  1.4  msec  of  flow,  the 
temperature  difference  between  the  first  drops  that  entered  the 
system  and  the  chamber  wall  ranges  from  50  to  70°F  for  the  three 
classes  of  drops.  This  driving  force  leads  to  heating  of  the 
drops  when  they  impinge  on  the  chamber  walls,  the  extent  of  the 
heating  depending  on  several  factors.  Two  cases  are  considered. 

In  one  case,  the  drops  are  brought  immediately 
to  the  chamber  wall  temperature  once  they  have  been  in  the 
system  for  one  residence  time  period.  The  drops  then  continue 
to  undergo  evaporation  for  a  second  residence  time  period,  at 
the  end  of  which  t he  drops  are  made  inactive.  This  model  cf 
heat  addition  to  the  drop-gas  system  crudely  simulates  the 
following.  A  drop  undergoes  evaporation  during  the  time  required 
for  it  to  travel  to  the  chamber  wall.  When  the  drop  strikes 
the  chamber  vail,  the  drop  is  warmed  to  the  wall  temperature  and 
then  bounces  off  the  wall  back  into  the  thrust  chamber  where  it 
undergo^*-  further  evaporation.  Evaporation  continues  for  another 
residence  time  period,  at  which  time  the  drop  passes  through 
the  nozzle.  The  chamber  wall  is  considered  an  infinite  heat 
source  and  an  infinite  heat  transfer  rate  from  the  wall  to  the 
drop  is  implied. 


The  chamber  pressurization  curve  calculated  for 
these  conditions  is  given  in  Figure  11a  and  is  labeled  "wall 
heating  once".  This  curve  coincides  initially  with  the  "no 
wall  heating"  curve.  Once  the  first  drops  to  enter  the  system 
reach  the  chamber  walls  (1.4  msec),  the  two  pressure  curves 
begin  to  diverge  because  of  the  greater  evaporation  from  the 
"heated" drops .  The  ensuing  difference  between  the  two  calcu¬ 
lated  curves  is  a  measure  of  the  amount  of  heat  added  to  the 
vapor-drop  system  by  the  chamber  walls.  Agreement  between  the 
experimental  and  calculated  curves  is  Improved  substantially. 

The  second  case  of  heat  addition  considered 
simulates  drops  sticking  to  the  chamber  walls  rather  than 
bouncing  off  the  walls  after  being  heated.  In  this  case,  the 
drops  are  again  brought  to  the  wall  temperature  once  they 
strike  the  wall.  But  here  the  drop  temperatures  are  maintained 
at  the  wall  temperature  while  the  drops  undergo  continued 
vaporization  until  their  radii  decrease  to  one-tenth  their 
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initial  values  (99.9%  of  the  mass  of  the  original  drop  evapor¬ 
ates).  However,  since  the  vaporisation  is  kineticaily  con¬ 
trolled!  ,  this  case  is  not  the  same  as  assusing  eoaplete 
vaporisation,  as  will  be  shown  in  the  following  section.  In 
the  cslcnlation,  the  shape  of  the  drop  is  aaintained  spherical, 
that  is,  although  the  drop  sticks  to  the  wall,  the  shape  of 
the  drop  is  not  altered.  1 

The  chamber  pressurization  curve  calculated 
for  the  condition  of  "constant  wall  heating"  is  also  given  in 
Figure  11a  and  is  labeled  as  such.  Clearly  the  experimental 
curve  is  surrounded  by  the  two  calculated  cases  of  heat  addition 
to  the  system.  One  concludes,  therefore,  that  despite  the 
rather  short  tines  involved  during  chamber  pressurization  due 
to  propellant  vaporisation,  the  vapor-drop  system  ic  nca- 
adiabatlo.  Additional  experimental  evidence  verifies  quali¬ 
tatively  this  conclusion.  Vith  a  Plexiglas  chamber  and 
stainless  steel  nozzle  precooled  with  dry  ice,  the  pressure 
in  the  chamber  rose  during  the  first  residence  tine  period 
only.  \ 

\ 

A  calculation  wa &  aide  in  whxch  the  diops  were 
hi ought  to  the  wall  temperature  twice,  once  after  the  first 
residence  tixe  period  and  again  after  a  second  residence  time 
period.  Evaporation  continued  for  a  tfci.'d  residence  time 
period  at  the  end  of  which  the  drops  **are  made  inactive  The 
calculated  curve  is  given  in  Figure  lln  together  with  the 
experimental  curve  given  in  Figure  17  a.  The  good  a*rr«eaent 
between  the  curves  (10%)  is  a  measure  of  the  amount  of  heat 
additioh  to  the  vapor-drop  system  and  its  time  dependency. 
However,  this  mode'  of  heat  addition  is  of  course  unrealistic 
since  some  oplatterlag  of  the  drops  would  result  wben  they 
strike  a  Fall,  some  of  the  liquid  was*  sticking  to  the  wall 
and  some  bouncing  off  with  very  little  heating. 

h.  Summary  of  Talueo-Steady  State  Pressure  anl 

Temperature 

A  comparison  of  steady  state  values  calculated 
for  CCI4  according  to  various  models  of  chamber  pressux ization 
is  given  la  Table  IT.  An  experimental  value  of  the  steady 
state  chamber  pressure  resulting  from  CCI4  vaporization  in  the 
engine  under  consideration  is  included  also. 

Best  agreement  with  the  experimental  value  of 
chamber  pressure  is  obtained  with  the  "wall  heating  twice" 
model  giver  above. 


TABLE  IV 


Approximate  average  for  all  active  drops  ia 
Calculated  tram  Measured  P_  assuaiag  Tg  -  5: 
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The  lowest  calculated  value  of  steady  state 
chamber  pressure  results  from  a  model  in  which  the  temperature 
of  th®  system  (both  vapor  and  condensed  phases)  is  assumed  to 
be  that  temperature  which  gives  a  vapor  pressure  equal  to  the 
pis  pressure  in  the  chamber.  That  is,  phase  equilibrium  pre¬ 
vails  in  an  adiabatic  system.  Equations  3  through  6  of  Refer¬ 
ence  1  were  used. 


If  one  assumes  complete  vaporization  of  the 
propellant,  in  effect  a  gaseous  feed,  very  high  steady  state 
chamber  pressures  result,  the  values  depending  somewhat  on 
the  assumed  temperature.  This  model  is  very  different  from 
t*e  computer  calculated  case  for  constant  wall  heating  of 
the  drops,  once  they  strike  the  chamber  wall.  The  latter 
model  is  kinetically  controlled  which  results,  in  the  present 
case,  in  only  15%  of  the  propellant  flow  being  vaporized, 
as  given  in  the  table. 

In  summary,  the  model  which  best  predicts  thrust 
chamber  pressurization  due  to  propellant  vaporization  is  the 
modal  discussed  herein  wnich  is  based  on  physical  kinetics 
and  a  non-adiabatic  system.  For  two  prlnicpal  reasons, 
however,  the  model  as  it  presently  exists  is  not  complete. 

The  first  reason  is  that  the  heat  addition  aspect  needs  to 
be  made  general.  Secondly,  a  comparison  between  experimental 
and  calculated  pressure  curves  for  considerably  rore  volatile 
prop-liants  such  as  Nj>04  and  Compound  A  shows  that  a  finite 
time  is  required  for  the  flow,  once  initiated,  to  build  up  to 
the  full  flowrate.  Substantial  "flashing"  of  these  propellants 
occurs  within  the  injector  volume  upon  propellant  valve  opening. 
Therefore,  only  vapors  emerge  from  the  injector  initially,  then 
a  vapor-liquid  mixture,  and  finally  all-liquid  flow  is  achieved. 
During  the  vapor  and  vapor-liquid  periods  of  flow,  the  mass 
flowrate  is  substantially  below  the  full,  nominal  flowrate. 
Consequently  the  pressure  in  t  chamber  rises  more  slowly  than 
would  be  thb  case  without  the  vapor  and  two-phase  flow  condition. 

C.  Reaction  Kinetics 


The  prediction  of  ignition  oelays  in  hypergolic  systems 
necessitates  an  accurate  knowledge  of  the  rate  constants  of 
the  reaction  leading  to  ignition.  A  first  step  toward  this 
goaX  is  the  determination  of  the  overall  reaction  rate  as  a 
function  of  the  concentrations  of  reactants  and  temperature, 

A  second  step  is  to  asssss  the  exact  mechanism  of  reaction. 
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This  part  of  the  report  describes  the  progress  which  has 
been  accomplished  toward  the  solution  of  these  two  problems. 

It  is  organized  as  follows.  First  a  background  of  thermal 
explosion  theories  with  regard  to  their  use  for  the  deter¬ 
mination  of  reaction  kinetics  is  presented.  These  considera¬ 
tions  are  followed  by  a  description  of  the  experimental 
apparatus  and  procedure.  A  following  section  analyzes  the 
experimental  data  and  then  the  experimental  results  are 
discussed,  including  our  present  knowledge  a  product  that 
is  formed  by  preignition  reactions  of  MgO^NgK^-type  fuels. 

1.  Theory  of  Thermal  Explosions 

Thermal  explosions  may  be  expected  to  develop  in  » 
reacting  system  whenever  the  heat  liberated  from  exothermic 
reactions  exceeds  the  rate  of  heat  dissipation  by  conduction 
or  by  convection.  Because  of  the  exponential  dependence  of 
the  reaction  rate  on  temperature,  the  rate  increases  rapidly 
as  the  temperature  rises  until  an  explosion  results. 

The  quantitative  mathematical  txeatment  of  thermal 
explosions  has  been  developed  mainly  by  Semenoff  (Ref.  4), 

Todes  (Ref.  5),  Frank-Xamenetskii  (Ref .  6),  Rice  (Ref.  7), 

Chambre  (Ref.  8),  and  Thomas  (Ref.  9,  10).  The  progress  in 
the  theory  has  consisted  mainly  of  finding  methods  of  inte¬ 
gration  of  the  differential  heat  equation  which  governs  the 
process.  A  good  review  of  the  use  of  explosion  limit 
phenomena'  for  elucidation  of  reaction  mechanism  is  given  by 
W.  Roth  and  D.  Scheer  (Ref.  11)  in  the  Advances  in  Chemistry 
Series  of  the  A.C.S. 

Two  questions  are  generally  considered  in  the  theory. 
The  first  is  to  find  the  so  called  critical  conditions  for 
ignition.  Th.'s  aspect  of  the  theory  does  not  consider  the  time 
variable.  All  one  is  concerned  with  is  the  minimum  values  of 
the  parameters,  (pressure,  concentrations,  temperature,  charac¬ 
teristic  longth)  above  which  an  explosion  develops.  This 
problem  is  relatively  simple  and  has  been  solved  analytically 
for  a  number  of  configurations  (slab,  cylinder,  and  sphere). 

An  will  be  seen,  the  critical  conditions  for  ignition  (or 
explosion  limits)  say  serve  to  determine  the  global  kinetic 
factors  of  the  reaction. 

The  second  question  concerns  the  time  Interval  which 
elapses  prior  to  the  development  of  a  thermal  explosion.  This 
problem  is  more  complex  and  analytical  solutions  are  known  only 
for  the  case  of  a  reaction  of  aero  order.  A  numerical  solution 
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has  been  given  by  Rice  (Ref.  7)  in  the  case  of  a  first  order 
reaction.  No  attempt  to  treat  the  case  of  reactions  of  higher 
order  has  ever  been  made. 

In  the  following,  we. will  be  concerned  only  with 
the  explosion  limits.  We  will  discuss  their  use  for  the 
establishment  of  reaction  kinetics. 

a.  Explosion  Limits 

Let  us  consider  a  reactive  mixture  contained 
in  a  closed  vessel  of  constant  wall  temperature  T0.  In  this 
case  the  transfer  of  heat  is  only  by  conduction  as  convection 
currents  are  assumed  to  be  negligible.  The  appropriate  equation 
for  thermal  balance  between  the  heat  generated  by  the  chemical 
reaction  and  heat  conducted  away  is  then: 

AV2T  -  -QW  (29) 


where 

«  thermal  conductivity  of  mixture  (cal  sec"’1cm“1<>K“1) 
T  -  gas  temperature  (°K) 

Q  -  heat  of  reaction  (cal  mole”1) 

W  -  reaction  velocity  (mole  sec-lcnT^) 

V2  •  Laplacien  operator  (cm”2) 


We  have  now  to  find  a  solution  of  (29)  under  the 
given  boundary  conditions.  As  long  as  such  a  solution  exists, 
the  system  is  stable.  The  value  of  the  parameters  at  which  such 
a  solution  becomes  impossible  is  then  taken  as  a  condition 
of  infladdiation.  For  a  wide  class  of  reactions,  one  may  write 
the  ArrhefciuS  reaction  rate  expression  as 

W(ca,  cb,  T)  -  A  ca  ®  cb  11  exp  (-  )  (30) 

H 


and  A 

is  the  sb-called  frequency  factor.  As  a  first  approximation 


where  c  and  cb  are  concentrations  of  reactants  A  and  B. 
respectively.  S  is  the  energy  of  activction  (cal  mole-1) 
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we  will  consider  A  as  constant41.  The  symbols  a  and  n  are  the 
partial  order*  relative  to  A  and  B, respectively .  By  definition 
the  total  order  of  reaction, N, is  N  -  a  +  n. 

Substituting  expression  (30)  for  W  in  equation 
(22)  balance  equation  becomes: 

V2T  -  A  ca  m  cb  °  exp  (  ^  )  (31) 

Analytical  solutions  of  equation  (31)  can  be  found  if  the  two 
following  approximations  are  made. 

1)  The  exponent  in  the  term  is  expanded 

as  first  proposed  by  Todes  in  a  power  series  (T-T0)/T0  (where  Tc 
is  the  wall  temperature)  and  powers  higher  than  one  are  neglected. 
The  following  expression  results: 

-JL  -  -JL 

RT  RTC 

It  is  worthwhile  to  note,  sin'**  the  expansion 
is  valid  only  for  values  of  T/T0  close  to  unity,  that  Todes 
approxiaatlon  implies  RTQ/£  «  1. 

2)  The  concentrations  of  reactants  are  assumed 
to  vary  only  slightly  durin*  the  induction  period  so  that  the 
effective  concentration  of  A  and  B  at  the  ignition  point  may  be 
approximated  with  the  initial  concentrations  clo  and  c^.  , 
respectively. 


Under  the  above  assumptions  it  has  been  shown 
first  by  Prank-Kamenetskil  that  the  criterion  for  thermal  igni¬ 
tion  Is: 


s 


ai  r2  cao*  c^H  A  exp  (-  - — ) 
A  RTd2  RT0 


(32) 


Strictly  speaking  A  is  only  constant  in  the  case  of  a 
reaction  of  the  first  order.  When  the  order  of  reaction 
is  superior  to  one,  A  is  a  weak  function  of  temperature. 
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pressure,  P, 
B, 


Expressing  cftC  and  cbo  in  terms  of  total 
and  initial  mole  fractions  and  Xbo  ot  A  and 


c  * 

S.O 


'ae 


Rbo 

RTft 


RT 


Ound 


P 

RT0 


where 

-  initial  moles  of  A  in  the  vessel 
ao 

V  volume  of  vessel 

p#o  «  partial  pressure  of  A 

Pbo  *  partial  pressure  of  B 

Substituting  these  values  of  c*0  and  cbo  in  ^32)  *ivcs: 

6  *  f  fv  r‘  x~-  >C  (fcj *  (**:)  <S3’ 


and 

6 

1.  2 


In  this  expression  r  is  a  characteristic  length 
S  is  a  non-dimensional  constant.  The  numerical  value  of 
depends  on  the  geometry  of  the  vessel  and  has  the  value  of 
and  3.32  in  the  case  of  a  slab,  cylinder  and  sphere, 


respectively. 


b.  Determination  of  Reaction  Kinetics  from  Explo¬ 
sion  Limits 

Equation  (33)  is  the  basic  equation  used  for 
the  determination  of  reaction  kinetics  from  explosion  limits 
data.  In  the  following  sections  three  categories  of  explosion 
limit  data,  depending  on  the  parameter  which  is  maintained 
constant,  are  considered. 


i.  Overall  Order  of  Reaction,  W 


If  the  temperature  and  composition  of 
the  mixture  are  maintained  constant  and  the  vessel  diameter 
is  varied,  equation  (33)  can  be  written  in  the  form 
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•*  * 


P  f 


R  5  a  exp 


Ms! 


E  aAX 


40 


cor\/CaM"t 


or 


In 


i  *  — 


m  r 


c  vrwt  an  i 


(34) 


Thus,  «  plot  of  In  Pvb  In  r  yields  a  straight  line  fros  whose 
slope  one  can  determine  H,  the  overall  order  of  reaction. 


ii.  Partial  Orders  of  Reaction,  n  and  ■ 

If  the  temperature  and  reactor  size  are 
maintained  constant  and  the  composition  is  varied,  equation  (33) 
can  be  written  as 


with  conltat*! 


conjTant 


~~  He-*  -Nil  r  , 

R  5  A  exp 
£  Ok  f\rt- 


It  can  easily  be  shown  that  the  pressure  limit  has  a  minimum 
for  a  value  of  X*0  equal  to  5,  allowing  a  determination  of  the 
partial  orders  a  and  n  if  tbBtotal  order  of  reaction,  M,  is 
known. 


iii.  Activation  Energy,  E 

If  the  composition  and  the  reactor  size  are 
maintained  constant  and  the  intial  temperature,  T0,  is  varied, 
equation  (33)  can  be  written  as 


_E _  J_  g  \  R _ 

NRT.  N  £  d  A  !•*  x~  X;. 
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Thus,  a  plot  of  In 


(p/ Te'+',N) 


vs.  y 


giyes  with  a  good 
E 


approximation  a  straight  line  with  slope  equal  to  .  This 
permits  a  determination  of  the  activation  energy,  fry  if  the  over¬ 
all  order  of  reaction,  N,  is  known. 


To  summarize,  the  overall  order  of 
reaction  N  can  be  obtained  from  the  effect  of  vessel  size  on 
ignition  pressure  limits.  Also,  the  composition  effect  allows 
a  determination  of  the  partial  orders  relative  to  each  reactant. 
Finally,  the  effect  of  temperature  allows  a  determination  of 
the  activation  energy. 


2.  Experimental  Apparatus 

The  considerable  reactivity  of  N204  with  the  hydra- 
zine-type  fuels  requires  that  an  experimental  technique  be  used 
that  slows  the  ignition  reaction  to  a  point  which  permits  valid 
measurements to  be  made.  Several  techniques  are  possible:  low 
temperatures,  considerable  dilution  by  an  inert  substance,  or 
low  pressures.  The  latter  was  chosen  for  several  reasons. 

Since  the  reactions  of  Interest  are  gas  phase  reactions,  low 
temperatures  are  ruled  out.  Dilution  by  an  inert  substance  is 
undesirable  in  general  for  kinetic  investigations  due  to 
possible  "third-body"  effects,  or  inf luence,  on  reaction  mechan¬ 
isms.  The  low  pressure  technique  was  adopted  so  that  the 
ctesired  gas-phase  reactions  of  the  undiluted  reactants  could  be 
studied  at  temperatures  within  the  range  encountered  during  the 
start-up  of  operational  reaction  control  systems. 

The  apparatus  used  for  the  determination  of  the 
kinetic  factors  required  in  the  mathematical  model  is  shown 
schematically  in  Figure  12.  It  consists  of  vaporizing  and 
flow  metering  sections,  a  flow  reactor  and  a  low  pressure 
bousing  with  observation  windows  and  instrumentation. 

a .  Vaporization  of  Propellants 

i.  Fuels:  Since  the  substances  studied  are 
liquids  at  room  temperature ,  whereas  vapor  is  desired,  it  is 
necessary  to  provide  a  means  for  vaporizing  the  propellants. 

All  four  hydrazine-type  fuels  investigated,  H2H4.  MMH,  UDMH, 
and  50-50,  are  hazardous  and  could  give  rise  to  explosive 
decomposition  when  vapors  are  in  contact  with  metal.  Therefore, 
an  all  glass  vaporization  system  was  used.  The  liquid  was 
vaporized  in  a  glass  coil  heated  at  90°C,  and  was  allowed  to 
expand  in  a  1  liter  flask  maintained  at  the  same  temperature. 
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The  prv»#urv  in  the  veporixor  whs  monitored  by  a  mercury 
M9999ter  as  sfc'vv'n  la  Figure  12,  The  vapors  were  then  admitted 
te  an  sll-glss** heated  line  connected  to  the  flow  reactor. 

i.i .  >304 :  Gaseous  MoOa  was  taken  directly 

from  a  cylinder  ikjMrwnrin  a  constant  temperature  bath  at 
JC°C.  The  vapor  pressure  at  thi*  ten  per at u re  is  a  little 
above  1  atmosphere .  The  heat  source  was  sufficient  to  saint a in 
steady  floe  conditions  up  to  0,2  gr  /sec.  This  flow  was 
adequate  for  the  purpose  of  the  experiments. 

111  •  Flowraticg  Systems 

i.  Fuel:  The  hydrasines  wmre  flow  rated  in 
the  liquid  phase  with  calibrated  rotameters  prior  to  entering 
the  vaporiser  described  above.  Sach  of  the  fuels  was  calibrated 
individually.  It  took,  in  general,  a  few  sinutes  before  steady 
flow  conditions  were  reached  in  the  whole  flow  system.  Stabilized 
flows  were  established  in  each  test  before  data  were  taken. 

it.  BaQf:  Kitrogen  tetroxide  was  flow  rated 
in  toe  gaseous  form  at  a pressure  slightly  below  ambient 
tc«G 0  m  Hf.) .  The  rotameter  was  maintained  at  constant  tempera¬ 
ture  slightly  beve  axbient,  30°C.  Under  these  conditions  the 
vppoes  are  «  fixture  of  and  ITOg.  Therefore,  the  flcwrator 

wae  calibrated  in  terms  of  mbs  flow  by  collecting  the  vapors 
in  a  condenser  at  liquid  nitrogen  temperature  and  weighing  the 
Bess  collected  during  a  measured  tins  period. 

c .  The  Flow  Reactor  \ 

TLie  flow  reactcr  shown  schematically  in  Figure  12 
was  a  pyr  >x  tube  open  at  one  end  to  the  ambient  atmosphere  of 
the  low  pressure  housing  described  below.  The  standard  flow  tube 
In  which  most  it  the  experiments  have  been  conducted  had  ia 
diameter  of  4.7  cm  and  a  mixing  length  of  51  cm.  Two  other 
diameters  of  2.2  cm  and  6.6  cs  also  were  used.  The  flow  tube 
was  located  in  an  insulated  housing  and  could  be  heated  elec¬ 
trically  over  ita  entire  length  up  to  temperatures  of  300°C. 

Hefore  entering  the  flow  tube,  the  reactants  were  preheated  to 
the  reactor  temperature  and  were  premixed  in  a  small  section 
whos*  diameter  is  approximately  1/3  that  of  the  reactor. 

d.  Concentric  Tube  Arrangement 

Scat  ignition  experiments  were  perforned  using 
the  concentric  tube  arrangement  shown  in  Figure  13.  The  concen¬ 
tric  tube  arrangement  differed  from  the  flow  reactor  discussed 
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Fuel 


Figure  13.  Concentric  Tubf Reactor 


"»F! 


above  in  that  the  exits  of  the  inner  and  outer  tubes  were  in 
the  same  plane  and  mixing  took  place  in  the  vacuum  tank.  These 
experiments  were  not  intended  for  reaction  kinetic  study  but 
were  performed  in  order  to  compare  the  hypergolicity,  of  the 
N204/hydrnzines  combinations  with  the  F2/H2  system  on  the  same 
scale.  The  oxidizer  was  fed  through  the  inner  tube  (1.2  cm 
in  diameter)  and  emerged  as  a  free  jet  in  a  concentric  flow 
of  fuel  (4.7  cm  in  diameter). 

e.  Low  Pressure  Housing 

The  low  pressure  housing  was  made  of  stainless 
steel  for  corrosion  resistance.  It  is  a  cylinder  three  feet 
high  and  12  inches  in  diameter.  The  reaction  tube  and  insulated 
housing  are  mounted  in  a  stainless  steel  tube  which  is  water 
cooled  to  prevent  heat  transfer  to  the  surrounding  shroud  gas. 
The  reactor  assembly  is  installed  through  the  lower  flange  of 
the  low  pressure  vessel  by  means  of  a  vacuum  feed-through  and 
can  be  moved  up  and  down  automatically.  The  pressure  in  the  low 
pressure  vessel  is  regulated  both  by  a  throttle  valve  in  the 
connecting  line  to  the  high  capacity  pump  and  by  an  appropriate 
bleed  of  shroud  air  or  nitrogen  which  has  the  additional  func¬ 
tion  of  diluting  the  corrosive  ga3es. 
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£.  Temperature  Measurements 

The  temperature  in  the  flow  tube  or  in  the 
exhaust  stream  is  measured  by  means  of  Pt  .  Pt/Rh  thermocouples 
aade  of  wires  of  0.2  bub  in  diameter.  In  one  experiment  the 
thermocouple  bead  was  covered  with  a  coat  of  Si02  to  prevent 
the  possibility  of  catalytic  effects.  As  the  readings  did  not 
differ  from  those  obtained  with  the  bare  thermocouple  it  is 
concluded  that  no  catalytic  surface  heating  of  the  thermo¬ 
couple  took  place  and  that  the  thermocouple,  coated  or  uncoated, 
indicated  the  true  temperature  of  the  gas.  Subsequently,  most 
experiments  were  performed  with  an  uncoated  thermocouple. 

g.  Gas  Sampling 

The  sampling  probe  was  made  of  a  pyrex  tube, 

0.8  cm  in  diameter,  which  had  been  drawn  At  the  extremity  to 
reduce  its  diameter.  The  diameter  of  the  throat  was  0.2  mm. 

As  it  was  not  intended  to  quench  the  reactions,  the  throat 
diameter  was  not  of  critical  importance.  The  gases  were  with¬ 
drawn  into  a  pre-evacuated  vessel  and  were  later  analyzed  mass- 
spectroscopically . 

3.  Experimental  Results 
a.  Ignition  Limits 

.  * 

i.  Flow  Reactor:  The  procedures  for  the  flow 
reactor  tests  are  as  follows.  T5e  apparatus  was  first  evacuated. 
Then  a  bleed  of  nitrogen  of  about  1  cfm  was  admitted  into  the 
pressure  vessel  and  by  adjusting  the  throttling  valve  in  the 
vacuum  line  a  pressure  of  about  0.5  mm  Hg  was  maintained  in 
the  vessel.  Gaseous  fuel  and  oxidizer  were  then  allowed  to 
flow  into  the  reactor.  When  steady  flow  conditions  were  reached 
the  ambient  pressure  in  the  chamber  was  gradually  raised  by 
bleeding  air  into  the  vacuum  line  through  an  auxiliary  valve. 
Ignition  occurred  at  some  minimum  pressure  and  was  observed 
visually.  It  was  somewhat  difficult  to  distinguish  if  ignition 
occurred  in  the  tube  or  in  the  free  jet.  However,  the  occur¬ 
rence  of  ignition  was  unmistakeable  as  it  resulted  in  a  visible 
flame  either  anchored  at  the  open  end  of  the  reaction  tube  or 
flashing  back  and  forth  in  the  tube. 

All  experimental  data  obtained  with  the 
hydrazine-type  fuels  are  reported  in  Tables  V,  VI,  VII,  and 
VIII.  Minimum  ignition  pressures  for  the  vapor-mixtures  of 


* 
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NO 2 /Mifll,  N02/UDMH,  N02/50-50  and  NO2/N2H4  are  given  in  the 
respective  tables.  The  important  parameters  from  a  kinetic 
point  of  view  are,  as  discussed  above,  the  equivalence  ratio 
0,  the  initial  temperature  of  the  reactor,  To,  and  the  reactor 
diameter.  These  arc  listed  respectively  in  Columns  I,  II, 
and  III.  Column  IV  gives  the  length-to-diameter  ratio  (L/D) , 
and  Column  V,  the  mean  flow  velocity  of  the  gas  in  che  reactor. 
The  minimum  ignition  pressures  are  reported  in  Column  VI.  The 
order  in  which  the  data  are  reported  in  each  table  has  been 
organized  so  as  to  show  the  effect  of  the  individual  parameters 
listed  above. 


ii.  Concentric  Flows 


Minimum  ignition  pressures  using  the 
concentric  flow  arrangement  of  Figure  13  were  obtained  for 
NO2/MMH  and  for  F2/H2.  The  procedure  with  NO2/MMH  was  similar 
to  that  used  for  the  determination  of  the  ignition  limit  in  the 
flow  tube  arrangement.  First,  the  desired  flows  were  established 
with  a  low  ambient  pressure  in  the  low  pressure  housing.  Then 
the  pressure  was  gradually  raised  until  ignition  occurred.  The 
result  of  the  ignition  was  a  diffusion  flame  which  burned  at 
the  exit  of  the  inner  tube. 

With  F2/H2  the  experiments  were  performed 
in  the  large  tank  used  in  Task  III  of  the  program.  The  flows 
of  H2  and  F2  were  metered  by  means  of  calibrated  orifices  in 
a  manner  similar  to  that  described  for  the  impinging  stream 
tests.  The  control  valves  were  preset  for  a  running  time  of 
1  second.  Each  experiment  was  conducted  with  a  given  pressure 
in  the  tank.  The  limit  of  ignition  was  found  by  trial  and 
error.  Here  again  the  ignition  resulted  in  a  diffusion  flame 
anchored  at  the  open  end  of  the  inner  tube.  The  minimum 
ignition  pressures  with  the  pertinent  flow  parameters  are  given 
in  Table  IX.  It  is  seen  that  under  similar  flow  conditions 
F2/H2  has  an  ignition  pressure  lower  than  that  of  NO2/MMH. 

b.  Temperature  Profiles 

A  typical  longitudinal  profile  taken  along  the 
center  line  in  the  4.7  cm  diameter  tube  with  a  flow  of  3.2  meters/ 
sec  of  stoichiometric  NO2/MMH  at  5.1  mm  Hg  is  shown  in  Figure  14. 
The  figure  also  shows  profiles  taken  without  flow  in  the  reac¬ 
tion  tube,  one  (Curve  I)  immediately  before,  and  another  (Curve 
III)  immediately  after  the  experiment.  The  data  shows  clearly 
that  a  low  temperature  reaction  with  heat  evolution  is  taking 
place  before  ignition  and  corroborates,  therefore,  the  general 
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TABLE  V 


MINIMUM  IGNITION  PRESSURES  OF  N09/MMH  VAPORS 
IN  FLOW  REACTORS  z 
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TABLE  V  (continued) 


X.o 


2.5 


1.0 


h 


203 


298 


2.5 


1.0  298 


2.5 


2.5 


0.5 

0.75 

1.0 

1.5 

2.0 

3.0 


298 


D 

cm 


6.6 


6.6 


L/D 


4.8 


U 

a/sec 


0.45 

0.94 

2.17 

3.82 

4.71 


1.69 

2.20 

3.75 


9.5 


2.2 


2.2 


11.5 


30 


0.46 

1.21 

2.63 

4.05 

4.85 


1.94 
2.64 

3.94 


2.16 

3.25 

6.49 

7.82 

9.41 


4.54 
5.80 

9.54 


3.56 

7.20 


5.29 

7.28 

12.5 


5.25 

6.44 

8.33 

10.03 

10.05 


MOTES 


:  6  ■  r (7/0)  Actual 

TT/Cr  gfbich. 

w  — — 


P 

Ml  Hg 


6.0 

4.35 

3.60 

3.55 

3.6 


2.30 

2.65 

3.10 


5.85 

3.37 

3.10 

3.35 

3.50 


2.0 

2.2 

2.95 


11.3 

11.3 

11.3 

12.5 

13.0 


7.7 

8.9 

11.0 


10.3 

10.2 


6.6 

7.2 

8.4 


14.0 

11.4 

8.6 

7.3 

7.3 

>16.0 


^Equivalence  Ratio 


vol 


pavg. 
n»  Sg 


4.22 


2.68 


3.83 


2.38 


11.9 


9.2 


10.25 


7.4 


?o  ■  Initial  vail  temperature  of  flow  reactor 
D  a  Diameter  of  reactor 
L  ■  Length  of  reactor 
0  ■  Initial  gae  velocity 
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TABLE  VI 

MINIMIIM  IGNITION  PRESSURES  OF  NOo 

REACTORS 


t 

k 

D 

ca 

Ia/D 

0.25 

298 

4.7 

4.76 

H 

1.00 

. 

1.50 

2.00 

3.00 

4.00 

IH 

- 

m 

1.00 

341 

4.7 

4.76 

0 

a/aec 


3.98 

4.23 

6.40 

5.85 


2.98 

4.67 

8.31 


2.39 

4.78 

6.60 

7.50 

11.61 


4.38 

4.20 


3.44 

6.89 


2.91 

4.01 

7.72 


3.33 

4.59 


3.44 

4.59 


3.68 

4.95 


5.09 


5.13 

6.44 

8.64 

10.95 
13. 69 
17.21 
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VAPdB  IN  FLOW 


P 

■«  Hb 

m 

5.1 

4.8 

5.4 

5.9 

5.3 

4.4 

4.6 

4.4 

4.5 

3.0 

3.0 

2.9 

3.1 

3.5 

3.1 

2.45 

2.5 

2.5 

El 

2.3 

2.5 

2.6 

2.67 

2.3 

2.5 

2.4 

2.5 

2.5 

2.5 

2.6 

2.9 

2.75 

3.1 

1.6 

1.7 

1.9 

2.0 

2.0 

2.7 

1.98 
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TAELS  VI  (continued) 


0 

D 

ca 

L/D 

T  ] 

m/sec 

P 

«*  Hg 

Pavg. 
an  Hg 

1.00 

403 

4.7 

4.78 

19.4 

1.0 

j 

23.6 

1.1 

1.2 

24.9 

1.3 

1.00 

293 

4.7 

10.8 

2.17 

3.3 

3.54 

2.7 

4.78 

3.0 

2.9 

5.62 

3.4 

7.71 

3.1 

3. 00 

4.18 

1.6 

6.27 

1.6 

1.73 

10.03 

2.0 

4.00 

4.78 

1.6 

6.37 

1.8 

1.76 

8.25 

1.9 

1.00 

•  298 

2.2 

11.5 

6.33 

6.9 

9.77 

6.7 

12.7 

6.9 

6.9 

16.3 

6.7 

18.9 

7.2 

4.0 

13.2 

5.3 

5.6 

17.8 

5.9 

_ i 
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TABLE  VII 

MINIMUM  IGNITION  PRESSURES  OF  KOo/50-50  VAPORS  IN  FLOW 

REACTORS 


0 

£ 

D 

CB 

0 

n/sec 

P 

■■  Be 

El 

1.0 

268 

•4.7 

10.8 

2.30 

3.6 

2.77 

4.0 

3.68 

4.3 

3.93 

6.34 

3.33 

7.16 

3.85 

8.01 

4.3 

3.26 

. 

3.82 

2.0 

5.46 

2.1 

6.24 

2.45 

2.56 

6.20 

3.7 

■nnn 

2.2 

11.3 

7.8 

■If 

8.0 

8.03 

12.2 

8.3 

2.20 

10.3 

5.1 

11.3 

6.2 

6.33 

■ 

13.6 

7.7 

j  TABLE  VIII 

j  MINIMUM  IGNITION  PRESSURES  OF  NO2/N0H4  VAPORS  IN  FLOW 

REACTORS 
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TABLE  IX 

MINIMUM  IGNITION  PRESSURE  IN  TREE  JETS 


Type 

Oxidizer 
Flo* 
cm3 /sec 

STP 

Fuel 

Flo* 
cm3 /See 

STP 

Jet 

Velocity 

cm/sec 

OTP 

Ignition 
Pressure 
mm  Hg 

no2 

in  MMH 

12 

24 

520 

15.3 

F2  in  H2 

J 

50 

200 

7500 

12.0 

BB 

25 

300 

3750 

11.5 

Diameter  inner  tube  -  1.2  cm 
Diameter  outer  tube  -  4.7  cm 

STP  -  Standard  Temperature  and  Pressure 

OTP  -  Operating  Temperature  and  Pressure 


Ambient  Pressure 
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Figure  14.  Axial  Temperature  Profiles, 

N02/MMH 
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picture  of  thermal  ignition  which  is  postulated  in  this  work. 
However,  some  experimental  difficulties  for  the  quantitative 
study  of  the  reaction  kinetics  are  apparent.  Namely,  a  gradual 
increase  of  the  wall  temperature  is  taking  place  during  the 
run  as  can  be  seen  by  the  comparison  of  the  profiles  taken 
without  flow  before  and  after  the  experiment.  Secondly,  during 
the  time  needed  to  take  the  profile  with  the  traversing 
mechanism  (2  to  4  minutes)  a  condensate  of  low  vapor  pressure 
was  observed  on  the  walls  of  the  tube  and  also  on  the  thermo¬ 
couple.  The  importance  of  this  condensate  or  adduct  will  be 
discussed  further  later  in  this  report.  False  readings  of 
temperature  were  a  very  likely  possibility  under  these  conditions. 

In  order  to  eliminate  these  eventual  causes  of 
error,  it  was,  therefore,  decided  to  take  the  temperature  point 
by  point  allowing  only  a  minimum  running  time  so  that  the  wall 
temperature  during  the  run  would  have  changed  only  by  a  negli¬ 
gible  amount.  Such  a  technique  was  used  to  determine  the 
temperature  profile  in  the  tube  Just  prior  to  ignition.  With 
the  thermocouple  located  at  a  fixed  position  in  the  flow  tube, 
the  desired  propellant  flow  rates  were  established  at  a  low 
ambient  pressure  (  ^  1  mm  Hg)  in  the  stainless  steel  vacuum 
tank.  Ambient  pressure  then  was  increased  gradually  until 
ignition  occurred.  The  temperature  rise  was  recorded  as  a 
function  of  time  and  the  maximum  temperature  reached  Just  prior 
to  ignition  was  noted. 

The  results  of  these  experiments  are  reported 
in  Figures  15  and  16  which  give  the  temperature  profiles  obtained 
with  stoichiometric  NQ2/MMH  and  NO2/UDMH,  respectively.  It 
should  be  noted  that  the  temperature  rises  immediately  from  the 
mixing  point  and  that  it  reaches  a  maximum.  Such  behavior  in¬ 
dicates  that  the  reaction  starts  immediately  without  any  in¬ 
duction  period.  This  point  is  of  importance  with  regard  to 
the  mechanism  of  reaction  as  it  tends  to  indicate  a  pure  thermal 
mechanism,  without  chain  branching. 

c.  Analysis  of  Reaction  Products 

The  analyses  reported  herein  refer  to  the  products 
of  the  low  temperature  reaction  occurring  prior  to  ignition. 

No  analyses  of  the  end  products  of  the  explosive  reaction  were 
performed. 


i.  Gas  Phase 

The  probe  was  located  at  the  center  lj|ne 
and  exit  plane  of  the  4.7  cm  tube  of  L/D  -  10.  Only  the 
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Distance  from  Mixing  Point, cm 


lagure  15.  Axial  Temperature  Profile  in  Plow 
Tube,  NCL/MMH 


Mixture:  Stoich 

Pressure:  4.5.nw  Hg 
Tube  Bis;  4.7  cat 
Initial  Linear 
Floe  Velocity:  3.13  */sec- 
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stoichiometric  NO2/MMH  was  investigated  at  pre-ignition  pres¬ 
sures  (  ^  5  mm  flg) .  A  shroud  of  helium  gas  was  used  in  order 
to  eliminate  any  doubt  concerning  the  origin  of  No  in  case 
this  gas  would  have  been  present  in  the  reaction  products. 

The  result  of  the  mass  spectrometric  analysis  are  given  in 
Table  X.  It  is  seen  that  only  2 %  of  CO2  appears  in  the  pre¬ 
ignition  reaction  products.  No  nitrogen  was  found  and  no 
trace  of  MMH  could  be  detected.  It  must  be  noted  that  the 
composition  reported  cannot  by  any  means  represent  the  exact 
composition  of  the  exhaust  gases  as  the  reaction  was  not 
quenched,  and  in  reality,  continued  in  the  sampling  container. 
The  analysis  is  revelant,  however,  in  the  sense  that  no  nitrogen 
or  hydrogen  was  formed  and  that  only  CO2  appeared  as  an  end 
product  but  in  rather  small  amounts. 

ii.  Condensate 


Upon  mixing,  the  vapors  of  the  hydrazines 
with  NO2  give  a  large  amount  of  smoke.  This  has  been  reported 
by  others  previously.  In  our  apparatus  the  liquid  droplets 
condensed  further  on  the  walls  and  collected  finally  as  a 
viscous  yellow  liquid  at  the  base  of  the  flow  tube  (placed 
vertically) .  The  amount  of  condensate  was  the  greatest  in  the 
case  of  NO2/MMH  and  amounted  to  more  than  20%  by  weight  of 
the  initial  reactants  in  some  experiments. 

Results  of  the  elemental  analyses  of  the 
condensate  are  given  in  Table  X  for  MMH  and  UDMB  as  fuels. 

The  elemental  analyses  agree  reasonably  well  with  additive 
compounds  of  formula  MMHNO2  and  UDMH • NO2 . 

The  infrared  absorption  spectra,  Figures  17 
and  18,  show  the  presence  of  nitrate  and  amine  groups.  No  free 
MMH  or  UDMH  was  detected. 

d.  Thermocheaical  Data  for  Condensate 


Some  preliminary  thermochemical  data  for  the 
MMH  adduct  was  also  obtained.  The  heat  of  combustion  with  oxygen, 
as  determined  from  limited  bomb  calorimetric  tests,  is  3090 
cal/gr.  This  gives  a  heat  of  formation  of  the  adduct  of 
A  Hf  •  -0.8  kc&l/gr  based  on  the  elemental  analysis. 
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TABLE  X 


ANALYSIS  OF  REACTION  PRODUCTS  OF  PRE- IGNITION 
REACTION,  N02AtMH 


GASEOUS  PRODUCT  OF  STOICHIOMETRIC 

NO  o /MMH 

Al*  5  mm  tig 

N2°4 

98% 

co2 

2% 

Helium 

Trace 

MMH 

None 

ELEMENTAL  ANALYSIS  OF  CONDENSATE 

NOg/MHH 

Carbon 

14.6% 

Hydrogen 

8.4% 

Nitrogen 

46.4% 

Oxygen  (by  difference) 

30.9% 

no2/udmh 

Carbon 

23.0% 

Hydrogen 

9.1% 

Nitrogen 

26 . 8% 

Oxygen  (by  difference) 

41.1% 
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4.  Discussion  of  Reaction  Kinetics  Results 


a .  Ignition  Limits 

i .  Effect  of  Flow  Velocity 

Each  set  of  data  consists  of  pressure 
limits  for  a  number  of  flow  velocities,  all  other  parameters 
being  maintained  constant.  It  will  be  seen  that  the  minimum 
Ignition  pressure  has  a  tendency  to  increase  with  flow  velocity 
within  the  limit  of  the  flow  rates  investigated.  This  is  a 
normal  effect  as  it  is  known  (?.ef.  1)  that  there  is  an  upper 
flow  velocity  limit  called  blow  off  limit  above  which  no  stable 
flame  can  be  maintained  in  the  stream.  The  important  point 
is  that  in  approaching  this  limit  the  ignition  pressure  does 
not  vary  drastically  with  the  stream  velocity.  This  is 
especially  true  for  the  stoichiometric  mixtures.  However,  for 
the  fuel  rich  mixture  the  velocity  effect  is  somewhat  greater. 

ii.  Effect  of  L/D 

Disregarding  as  a  first  approximation  the 
effect  of  velocity,  mean  values  of  pressure  limits  over  the 
velocity  ranges  investigated  may  be  computed.  These  values  are 
listed  in  Column  VII  of  Tables  V  to  VIII.  This  averaging  pro- 
cedure  is  used  to  determine  the  effect  of  L/D  on  the  minimum 
ignition  pressures.  In  Table  XI,  the  mean  values  of  the  pressure 
limits  which  are  tabulated  as  functions  of  L/D  and  the  equiva¬ 
lence  ratio  are  compared  where  possible.  Three  sets  of  data 
exist  for  NO2/MMH  and  one  for  NOo/UDMH.  It  will  be  seen  first 
that  no  appreciable  effect  of  L/D  can  be  detected  for  the 
stoichiometric  mixture  (0  -  1).  For  rich  mixtures,  however, 
the  minimum  ignition  pressure  increases  when  L/D  decreases. 

iii.  Effect  of  Reactor  Size 


The  fact  that  L/D  and  velocity  effects  are 
negligible  for  stoichiometric  mixtures  makes  it  now  possible  to 
determine  the  true  effect  of  diameter  on  the  explosion  limits. 
This  is  important  since  it  permits  the  total  order  of  reaction, 

N,  to  be  deduced  from  these  data.  Since  the  L/D  has  only  a 
slight  effect  for  0-1,  it  is  possible  to  average  again  the 
pressure  limit  values  of  the  two  L/D's  for  each  reactor  size. 

The  resulting  averages  are  plotted  in  Figure  19  as  In  P  vs.  In  r. 
In  the  case  of  KO2/HMH, for  which  three  data  points  are  available, 
a  straight  line  is  obtained  with  a  slope  indicating  an  order  of 
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TABLE  XI 


EFFECT  OF  REACTGS  L/D  OH  MIHIMOM  IGRITION  PRESSURES, 

RO2/MMH  and  NO2/UDMH 


0.6  cm  Tube,  MKH 


\l/D 

4.8 

9.5 

1 

4.22 

3.83 

avg.= 

4.02 

2.5 

2.68 

2.38 

4.7  cm  Tube,  KMH 


\VD 

0X 

4.76 

10.8 

a 

5.33 

5.87 

m 

B 

fB 

3.57 

4.7  cm  Tube,  UDMH 


\l/d 

0  \ 

4.76 

10.8 

1 

m 

2.9 

BPS! 

KBl 

3 

N 

* 

• 

_ a 

1.73 

- 

4 

2.4 

1.76 

2.3  cm  Tube,  KMH 


\l/d 

0  \ 

11.5 

30 

1 

11.9 

10.28 

KTg.- 

11.0a 

2.5 

9.2 

7.4 

Rote:  NiniBUM  ignition 

pressure  units 
are  mm  Hg 
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reaction  close  to  2.  For  NOg/UDMH  and  N02/50-50,  only  two 
data  points  exist  but  the  lines  joining  the  two  points 
indicate  again  an  order  of  reaction  close  to  two. 

iv .  Effect  of  Equivalence  Ratio 

.  .  For  the  NO2/MMH  mixtures  two  sets  of  data 

are  available,  one  taken  in  the  4.7  cm  tube  with  an  L/D  of  4.76 
and  the  other  taken  in  the  2.2  cm  tube  with  an  L/'D  of  30.  Using 
the  averaging  procedure  described  above,  the  minimum  ignition 
pressures  are  plotted  in  Figure  20  against  the  equivalence  ratio. 
In  both  cases  the  minimum  is  found  on  the  fuel  rich  side.  It 
is  somewhat  difficult  to  define  exactly  this  minimum  because 
the  curves  are  quite  flat  but  it  can  be  estimated  at  about 
0  -  2.5,  i.e. ,  the  equimolecular  mixture  NO2  +  CH3N2H3.  In 
the  case  of  NO2/UDMH  as  indicated  in  Figure  21,  if  seems  again 
that  the  minimum  is  close  to  the  equimolecular  composition, 
which  is  0  -  4.  Thus,  for  these  two  combinations  the  data 
indicate  partial  orders  of  reaction  close  to  1  relative  to 
fuel  and  oxidizer. 

No  complete  set  of  data  are  available  for 
N02/50-50  and  }i02/^2^4’  The  on^y  statement  which  can  be  made 
regarding  the*e  combinations  at  the  present  time  is  that  the 
equimolecular  mixtures  ignite  somewhat  easier  than  the  other 
mixture  ratios  tested. 

v.  Effect  of  Temperature 

The  effect  of  initial  temperature  was 
investigated  for  NO2/MMH  and  NOo/UDMH.  The  measurements  were 
made  using  the  4.7  cm  tube  of  L/D  -  4.76.  Averaging  again  the 
minimum  ignition  pressures  over  the  velocity  ranges  invest!-  „ 
gated,  the  results  are  shown  in  Figure  22  in  a  plot  of  In  P/T_ 
against  1/Tp.  The  three  experimental  points  for  NO2/UDMH  fall 
on  a  straight  line  whose  slope  corresponds  to  an  energy  of  acti¬ 
vation  of  E  -  7.2  kcal/mole.  Two  points  only  were  taken  for 
NO2/MMH.  The  slope  of  the  line  indicates  an  energy  of  activation 
of  5.2  kcal/mole.  The  temperature  effect  for  N02/50-50  and 
NO2/N2H4  was  not  investigated. 

vi .  Comparison  of  Reactivities  of  Hypergolic 

Combinations 


This  comparison  is  best  given  by  the 
pressure  limits  obtained  for  the  stoichiometric  mixtures  in 
the  4.7  cm  reactor  of  L/D  -  10.  The  comparison  is  shown  in 
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Tablos  £11  which  includes  kls©  the  concentric  flow  testa.  Tba 
following  order  of  reactivity  is  indicated:  J^Og/TOMB  >  MO«/ 
50-80  >H02/»BI  >  NC>2/lf2H4 . 

b.  Beteraaiaation  of  the  Beat  Release  Factor  (AQ 

So  far  the  analyses  of  the  experimental  data 
on  explosion  lisits  has  permitted  us  to  establish  the  overall 
kinetics  of  the  reaction.  The  total  order  of  reaction  for  the 
two  system  completely  investigated  (NOg/MMH  and  ifOg/HDUH) 
has  been  found  to  be  close  to  2  and  their  partial  orders 
relative  to  fuel  nnd  oxidizer,  close  to  1.  It  is  now  possible 
to  write  the  rate  express4 on  for  reaction  rate  in  the  form: 


p 

W  -  A  cjcg  exp  v-  g-) 

whore  c^  and  r ^  represent  the  concentrations  of  fuel  and 
oxidizer,  respectively.  Substituting  this  expression  into 
equation  (32) ,  the  Frank -Kaaenetskil  criterion  for  thermal 
ignition  becomes: 

8  *  f  tv  ^  *•  **  *  *xp  (■  ity  <3s> 

In  this  expression  all  quantities  are  known  from  the  explosion 
Halts  data  with  the  exception  of  A  and  Q.  It  is  thus  possible 
to  calculate  the  product  A Q.  This  calculation  will  be  performed 
in  detail  for  the  case  of  flOg/BMK .  The  explosion  limit  is  taken 
as  5.5  ks  Hg  which  is  the  mean  value  obtained  ior  the  stoichio¬ 
metric  mixture  at  Tc  -  228°£  in  the  tube  of  4.7  cm  diameter. 
Taking  for  B  the  value  of  5.2  kcal/mole  and  taking  for  Tv  a 
value  of  4  x  10-3  cal  see-5,  est-1  d»gree~5,  the  following  data 
are  known: 


G  -  2 

T0  -  298°K 

r  -  2.35  cm 

Cj  -  0,84  x  10-7  asole/cc 

eg  “  2.1  x  10"7  aole/cc 
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?i  -  4  x  10“5  cal/sec-ca-°C 

E  «  5.2  x  10^  cal /mole 

Substituting  these  values  into  equation  (35) ,  the  product  AQ 
becomes : 


AQ-  1.7  x  lO*4  cal  cc  _ 

rnole^  sec 

The  product  AQ  may  also  be  evaluated  using  a 
quite  different  method.  Considering  the  temperature  profile 
at  the  center  line,  as  shown  in  Figure  15,  it  is  noted  that 
the  first  part  of  the  curve  is  quite  linear.  In  this  region 
it  might  be  assumed  as  a  first  approximation  that  heat  losses 
by  conduction  to  the  wall  are  negligible.  This  is  certainly 
true  in  the  central  region  of  the  tube  before  the  boundary 
layer  is  fully  developed.  Taking  for  the  kinematic  viscosity 
the  value  of  V  -  13.3  CE^sec-1  and  for  the  initial  linear 
velocity  the  value  of  UQ  -  245  cm/sec,  it  is  found  that  the 
boundary  layer  merges  in  the  center  at  a  distance  of  16.3  cm 
from  the  origin.  For  the  central,  initial  portion  of  the  tube 
the  heat  equation  can  be  written  in  the  form 


/><>  U  —  -  WQ  -  AQc,c*  exp  (36) 

with  />  -  density,  cp  -  specific  heat  and  U  -  velocity  at 
the  center  line 

Applying  this  equation  at  x  -  10  cm  from  the  origin  and  taking 
into  account  the  pressure  and  composition  of  the  mixture,  it 
is  found  that: 

2 . 4°C/cm 

U  -  375  ca/sec 

T  -  347°K 

c^  -  0.7©  x  10”7  mole/cm3 

c2  •  1.99  x  10“7  mole/cm3 
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Cp  -  0.24  cal/gr-°C 

/°  -  1.4  x  10“®  gr/cm^ 

These  values,  substituted  into  equation  (36),  give 


AQ  =•  3.4  x  1014  cal 

mole*5  sec. 

Considering  the  uncertainties  in  the  values  of  cp  and  A,  the 
agreement  between  the  two  methods  must  be  considered  as  very 
satisfactory.  It  should  be  noted  that  the  value  from  the 
temperature  profile  method  is  preferred  as  it  is  more  direct 
and  is  independent  of  the  approximations  involved  in  the  theory 
of  thermal  explosions. 

c.  Pre-ignition  Gas  Phase  Reaction  Mechanism 

In  defining  a  scheme  of  reaction  for  the  processes 
leading  to  hypergolic  ignition  in  the  gaseous  mixtures  of  the 
hydrazines  with  nitrogen  dioxide  the  following  facts  have  to  be 
accounted  for: 

(a)  A  low  temperature  reaction  with  heat  generation  takes 
place  in  the  gaseous  mixture  prior  to  ignition. 

(b)  The  energy  of  activation  of  this  reaction  is  low  and  of 
the  order  of  5  to  10  kcal  per  mole. 

(c)  The  total  order  of  reaction,  N,  is  close  to  2  and  the 
partial  orders  relative  to  fuel  and  oxidizer  are  close  to 
unity. 

(d)  Hie  reaction  starts  without  any  apparent  induction  period. 

(e)  Chemical  additives  tested  have  in  general  no  other 
Influence  than  the  small  influence  expected  from  a 
diluent. 

(f)  A  liquid  condensate  of  low  vapor  pressure  and  showing 
a  1/1,  F/0  molecular  composition  is  formed  prior  to 
ignition. 

The  fact  that  no  induction  period  is  observed  and  that  chemical 
additives  have  in  general  no  influence  on  the  rate  of  reaction 
leads  first  to  the  conclusion  that  the  low  temperature  reactions 
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are  not  of  the  chain  type.  Further,  the  order  of  reaction 
and  the  adduct  formation  point  toward  a  molecular  association 
reaction  between  fuel  and  oxidizer.  It  is  difficult  at  present 
to  decide  whether  this  association  reaction  is  a  necessary 
step  or  only  a  side  reaction  in  the  mechanism  leading  to 
ignition.  However,  we  will  venture  to  propose  the  following 
physical  picture  as  a  working  hypothesis. 

When  fuel  and  oxidizer  come  into  contact,  they 
first  form  an  adduct  which  further  dimerizes  to  form  the  con¬ 
densate  of  low  vapor  pressure.  These  steps  are  accompanied  by 
a  significant  heat  evolution.  The  heat  liberated  is  only  a 
fraction  (10%)  of  the  total  heat  liberated  by  the  explosive 
reaction  but  it  is  sufficient  to  bring  the  mixture  to  a 
temperature  level  where  decomposition  and  explosive  chain 
reactions  are  initiated.  It  is  well  known  that  nitrogen  tetroxidc 
forms  addition  compounds  with  most  organic  compounds  that  are 
able  to  donate  electrons.  Electron  donors  usually  form  a  2:1 
compound  which  may  be  represented  as 

DONOR 

\ 

N  —  N 

\ 

DON  OR 

It  is  well  possible  that  the  condensate  isolated  from  the 
NO2/MMH  reaction  could  be  a  compound  of  this  type  as  the 
elemental  analysis  gave  C,  H  and  N  in  the  proportions  of  14.6, 

8.1  and  46.4  which  is  in  reasonable  agreement  with  the 
theoretical  values  of  14.5,  8.25,  and  47.3  for  a  compound  of 
the  formula  (CH3N2H3*H02)2* 

Some  interesting  conclusions  concerning  the 
mechanisms  of  reaction  might  further  be  drawn  by  considering 
the  results  of  the  thermochemical  data.  The  data  were  obtained 
from  only  a  limited  number  of  tests  and  certainly  require 
verification.  Nevertheless,  we  will  use  the  existing  preliminary 
information.  From  the  determined  heat  of  formation 

AHf0  -  -0.8  Kcal/gram 

and  the  known  heats  of  formation 

AHf0  -  23.2  Kcal/mole  for  MMH(g) 

and  AHfo  -  7.9  Kcal/mole  for  N02(g) 
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we  find  for  the  reaction,  assuming  the  adduct  is 

MMH(g)  +  N02(g)  -  1/2  (CH3N2H3NO2) 2  (1)  +  AH  -  105  kcal 

Assuming  that  this  reaction  is  the  dominant  exothermic  process 
in  the  low  temperature  regime  prior  to  ignition,  the  value  of  Q 
appearing  in  the  heat  equation  (36)  and  (35)  might  be  identified 
with  -  AH  -  105  kcal..  Thus  taking  3.4  x  10*4  cai  cc  mole-2sec~^ 
for  the  value  of  AQ,  we  find  for  the  frequency  factor  of  the 
bimolecular  association  reaction 

A  -  0.3  x  1010  cc  mole“lsec“* 

Such  a  value  is  substantially  lower  than  the  collision  fre¬ 
quency,  Z,  which  for  moderate  size  molecules  is  of  the  order 
of  Z  -  5  x  1Q14  cc  mole~)-sec~l .  Thus,  only  one  efficient 
collision  occurs  for  every  50,000  collisions  between  the 
reactant  molecules  and  the  steric  factor  f  -  A/Z  has  a  value  of 
6  x  10“6.  In  most  bimolecular  reactions  of  free  radicals  with 
molecules,  0.1  <  f  <  1.  However,  certain  additional  reactions, 
for  instance  at  double  bonds,  have  low  f  values  (10-3  -  10-5) . 
Furthermore,  it  has  been  found  that  molecular  association 
reactions  have  in  general  low  A  values.  For  instance,  Kistiakowsky 
and  Stouffer  (Ref.  12)  found  for  the  A-factors  of  the  association 
of  hydrogen  bromide  and  hydrogen  chloride  with  2-methyl  propane 
values  of  1010*2  and  10H  mole'^-cc-sec-l  which  are  fairly  close 
to  the  value  obtained  for  the  reaction  between  N02  and  MMH  from 
the  preliminary  thermochemical  data  on  the  adduct. 

5.  Summary  and  Conclusions 

The  analysis  of  the  experimental  data  on  explosion 
limits  has  permitted  us  to  establish  the  overall  kinetics  of 
the  low  temperature  reactions  leading  to  ignition  of  the  two 
systems  N02/MMH  and  NO2/UDMH.  The  total  order  of  reaction  for 
these  two  systems  is  close  to  2  and  the  partial  orders  relative 
to  fuel  and  oxidizer  are  close  to  unity.  The  energies  of 
activation  are  5.2  and  7.2  kcal/mole  for  NO2/MMH  and  N02/UDMH, 
respectively.  Partial  data  for  the  systems  NC>2/50-50  and 
NO2/N2EU  have  been  obtained.  The  total  order  of  reaction  for 
N02/50-50  has  been  found  to  be  close  to  2.  The  minimum  ignition 
pressure  for  NO2/H2H4  has  been  determined  in  the  case  of  the 
stoichiometric  mixture,  which  is  also  the  equimolecular  mixture, 
and  indicates  a  reactivity  less  than  those  of  the  other  hydra¬ 
zines. 
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The  formation  of  an  adduct  of  low  vapor  pressure 
and  of  near  1/1  molecular  composition  relative  to  fuel  and 
oxidizer  has  been  observed  and  analyzed  in  the  case  of  NO2/MMH 
and  N02/UIMH.  A  condensate  was  also  observed  in  the  case  of 
NOg/CO-SO  but  it  was  not  analyzed.  Both  its  appearance  and 
quantity  formed  are  similar  to  the  NOg/MMH  adduct.  A 
considerably  lesser  amount  of  the  NO2/UDMH  adduct  collected  in 
the  liquid  state  at  the  base  of  the  burner  assembly. 

A  synoptic  representation  of  these  conclusions  is 
given  in  Table  XIII. 

D.  Statement  and  Results  of  Mathematical  Model  of  Hypergolic 

Tgnit ion 

In  this  section,  the  results  of  Sections  B  and  C  on 
chemical  kinetics  and  chamber  pressurization  due  to  propellant 
vaporization  are  combined  into  a  mathematical  model  of  hyper¬ 
golic  ignition  in  reaction  control  systems  at  space  conditions. 
The  model  and  the  assumptions  involved  are  discussed  first. 
Subsequently,  calculated  ignition  delays  are  given  and  compared 
to  experimentally  measured  values.  Finally,  overpressures 
which  occur  during  engine  start-up  transients  are  considered. 

1.  Hypergolic  Ignition  Model 

a.  Dominant  Ignition  Reactions 

It  is  assumed  in  the  mathematical  model  that  the 
dominant  ignition  reactions  are  gas-phase  reactions.  Some 
experimental  justification  for  this  assumption  comes  from  the 
work  of  the  preceding  ten-month  program  reported  in  Ref.  1. 

Figure  31  of  that  reference  is  a  log-log  plot  of  average  ignition 
delays  vs.  pressure  for  a  number  of  hypergols  tested  in  an 
unconflned  impinging  stream  apparatus  (i.e.  Impinging  stream 
injectors  without  thrust  chambers) .  Curves  for  the  N2O4/ 
hydrazine-type  fuels  combinations  are  reproduced  in  Figure  23 
of  the  present  report. 

Essentially,  the  curves  as  drawn  simply  connect 
the  averaged  data  points  at  each  pressure  by  straight  lines. 
However,  in  the  pressure  range  from  150  to  60  mm  Hg,  the  three 
data  points  for  N2O4/UDMH,  N2O4/MMH  and  N2O4/MHF-5  (a  mixed 
hydrazine  fuel)  fall  on  single  straight  lines.  This  linearity 
is  found  despite  the  fact  that  the  greatest  data  scatter  was 
observed  in  this  pressure  range  (see  Ref.  1).  1*2^4  /  50- 50 

shows  a  nearly  linear  relationship  in  the  log- log  plot  in  this 
pressure  range. 
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The  fairly  high  degree  of  linearity  in  thv  190 
to  00  am  Hg  presvu re  range  suggests  that  a  common  process  is 
rate  controlling  in  this  range.  Since  the  delays  at  one  at- 
no  sphere  are  longer  in  every  case  than  those  to  be  expected  from 
an  extrapolation  of  the  straight-line  low  pressure  curves,  a 
different  controlling  process  is  indicated  for  the  ignitions 
at  high  pressure. 

The  results  of  the  nearly  150  tests  aade  with 
these  hypergolic  combinations  using  various  iaplnging  stress 
injector  configurations  (injection  velocities,  ispingenent 
angles,  inpingenent  lengths  and  type  of  nanlfoldlng)  shoe  that 
nixing  processes  are  not  controlling  at  any  pressure  investi¬ 
gated.  Neither  the  controlling  process  at  one  atnosphere, 
nor  the  one  at  lover  pressures  is  propellant  nixing. 

Figure  24  is  a  log-log  plot  of  average  ignition 
delays  at  one  atnosphere  vs.  relative  volatilities  of  the  five 
hypergolic  conbinations.  Relative  volatility  as  used  here  is 
sinply  the  ratio  of  the  vapor  pressures  of  oxidizer  to  fuel.  The 
figure  shows  that  the  greater  the  discrepancy  betveen  fuel  and 
oxidizer  volatilities  the  longer  the  ignition  delay  at  one 
atmosphere.  Two  observations  are  suggested  by  the  staple 
correlation  of  Figure  24.  Since  the  oxidizer  is  the  sane  for 
the  five  conbinations,  the  abscissa  of  the  figure  is  effectively 
the  reciprocal  of  the  vapor  pressure  of  the  fuel.  The  Ignition 
delays  of  these  hypergols,  then,  are  seen  to  depend  on  the  vapor 
pi^osure  of  the  fuel.  Tho  figure  indicates,  therefore,  that  fuel 
vaporization  is  controlling  at  one  atnosphere.  Then,  since 
fuel  vaporization  is  controlling,  the  dominant  ignition  reactions, 
even  at  one  atnosphere,  must  be  gas-phase  reactions. 

At  reduced  ambient  pressures,  no  correlation 
exists  betveen  ignition  delay  and  relative  volatility.  Con¬ 
sequently,  the  controlling  process  at  low  pressures  is  not 
propellant  vaporisation,  a  result  compatible  with  the  curves 
of  Fleurs  23.  Since,  as  stated  above,  mixing  processes  are  not 
controlling  at  low  pressures,  one  concludes  that  gas-phase 
chemical  kinetics  is.  The  success  of  the  mathematical  model, 
which  assumes  that  gas-phase  chemical  kinetics  is  controlling, 
supports  the  conclusion. 

Further  support  comes  fros  a  comparison  betveen 
the  number  of  gas-drop  collisions  and  the  auaber  of  gas-gas 
collisions.  Due  to  tbs  .  ..go  differences  betveen  the  volatilities 
of  Ng04  and  the  hydraxiaa  fuels,  ve  will  consider  the  possibility 
of  reactions  resulting  from  collisions  of  gaseous  molecules 


.93- 


Avar ago  I  fill  t  ion  Dalny  at  1  At»o«pUr* 


^•lative  Volatility  at  looa  Twaperatur*, 
V.P.  0*/V.P.fu#i 


Flyura  %4 .  Iaflotfoca  of  Propallact  Volat lilt  lea 
oe  Ignition  Delay  at  0e«  Atnonpbore 


/iF*PL-T»~«5-aS7 


fitl  dnps  of  Mi  (gas-drop  eollisi  mi),  for  «x»apl«.  Ths 
freewency  of  vhooo  collisions  will  bo  compared  to  the  fre¬ 
quency  of  gaseous  ROy-itstona  MB  collision*  (gas- gas 
collisions) . 

Making  css  of  computer  run  CC14  #E,  fo*  which 
w  -  0.0916  lb^/sec,  M  “  <  *  10“-  hoc,  Tq  2  540®F, 

?c  *  1.015  x  10-2  ft-*,  sad  rB  -  75  si  cron*  (see  Table  III 
for  initial  radii  of  the  three  drop  size  spray  model) ,  so 
first  ctlcolaU  t  he  total  crop  let  surface  area  per  unit  chamber 
toIsm.  Be  select  a  tine  at.  0.8  0  nsec  for  which  P»  -  1.74  an 
T.  -  513°F  -  28£°K  and  the  approximate  average  radii  of  all 
the  drop*  ia  the  system  in  each  class  are 

?!  -  77  x  10“®  ft  (range  of  actual  ri's:  83xl0~®  to 

7Lx10-6  ft) 

?.  -  238  x  )0~®  £t  (range  of  actual  ry'a:  250x10“®  to 

333x10-6  ft) 

and  r3  -  705  x  10“®  ft  (range  of  actual  r3's:  720x10“®  to 

687x10-®  ft) 

8y  equations  (3) ,  (4)  and  (6c)  ant  noting  that  20  ensenbles 
have  undergone  evaporation  during  tae  0.800  nsec  period  (tine/ 

A  tB  *  20),  we  obtain  tb*j  total  eurface  area  of  the  drops 
in  each  claee.  Dividing  the  'emitting  areas  ly  the  ch»nber 
volume  fives  the  eurface  are ae  per  unit  volume  ior  each  of  the 
thrm?  aggregate  drep  classes: 

Total  Aj/Tj,  -  1.28  x  10“2  caVcc  for  Fj  -  77  x  i0-®  ft 

Total  %/Tc  -  0.91  t  1Q~2  cw2/cc  for  rj  -  236  x  10“®  ft 

and  Total  A3/Y,.  -  0.26  x  10“2  c*2/cc  for  ?3  -  70S  x  10“®  ft 

The  total  drop  eurface  area  par  cc  ie  therefore  3.15  x  10“2  cm2/ 
cc. 


The  concentration  of  gas  no locales  at  0.800  nsec 
when  S^g  -  1.74  ae>  Sg  and  Tg  «  285-oS  (above)  is 

molecular  concentration  -  c  -  SL?2 .  -  f tfo  - 

T  ITg 

5.20  x  10l®  gas  »c lecoles/cc  (37) 

where  -  Ivogadre  number  -  6.034  x  1022  aolecules/mole 
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By  the  klietltf  theory  «f  frees ,  the  auaber  ef  collisions  p*r 
unit  tiwa  &s*we*2  gas  «ole«ules  and  a  surface  of  unit  are*  is 

^gas-surface  -  c  u  *»  c  /&Tg  \  (33) 

*  */ 

where  c  is  molecular  concentration  sad  u  is  Mean 
'  Molecular  velocity, 

««*  • 

Multiplying  04.  (38)  by  the  total  drop  surface  area  par  cc  as 
doternined  above,  and  taking  iato  account  the  fact  that  the 
Taper  Molecules  are  CCi.4,  we  obtain,  for  7  -  285°E: 

o  > 

%a««*drop  -  9.2  x  1018  gas-drop  collisions/cc-sec  (39) 

SquwUon  (39)  is  tbs  awaber  of  OCI4  gas-drop  collisions  per 
unit  tine  per  unit  voluse  at  0.800  »sec  of  the  computer 
calculation. 

m  f 

ffe  assu  aow  that  the  vapor  Molecules  are 
I0o  and  the  drops  are  IhH.  To  account  for  the  difference  in 
Molecular  weights  of  CCI4  and  KO*j  'v  c  «& .  38),  eq.  (39)  is 
Multiplied  by  (ICC^/KjlOg)!/^.  ju  }  jince  tbo  nuwber  of 
drops  ia  each  class  is  inversely  proper ..renal  to  density 
(see  eq.  4),  c-.  (39)  is  nultipled  also  hjr  (  f*  CCI4/  fi  MSB) . 
The  ssssdit  gives  the  frequency  of  collisions  between  the  B02 
Molecules  and  HMH  drops  and  is 

n  ga«-KB8  drops  -  3.1  x  lO*9  gas- drop  (40) 

colliaiona/cc-sec 

1 1  r)”C  ,  /  ■ 

for  JP*  -  1.7e  mm  Hg  and  Tg  -  285°K 

it  <*'-0.1  v 

Ftom  the  kinetic  -heory  of  gases,  the  total 
gar*#a*  collisions  between  unlike  Molecules  per  unit  veluae 
and  per  unit  tine  is  given  by 

(',<1*  00H  .  0  ,  • 

2gas-gas  -  cf  *QX  (rt  +  r0*>2 


(xn 

-■  f  on  \  o> 
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whore  Cf  and  eox  are  molecular  concentrations  of  fuel  and 
oxidizer 

Yf  and  ros  are  molecular  radii 

k  “  Boltzmann  constant  “  1.38  x  10”^®  er g/°K  molecule 
and  cjf  and  mOK  are  molecular  masses. 

For  r^us  «  4  x  10”8  cm ,  rN02  “  3  x  10~®  cm,  and  T  -  285°K,  we 
obtain 

Sgas-gas  -  7.88  x  10-10  Cf  cox  gas-gas  collisions/  (42) 
cc-sec 

As  a  "worst  case",  we  assume  the  molecular  concentrations  of 
gaseous  MMH  and  N<?2  are  in  proportion  to  their  vapor  pressures 
at  the  temperature  of  interest  (285°K)  ,  cnOjj/cMMH  -  21.  Since 
the  total  molecular  concentration  is  5.90  x  1016  (eq.  37),  the 
total  gas-gas  collisions  between  unlike  molecules  is 

Zgas-ga3  -  1.2  x  10^3  gas-gas  collision/cc-sec  (43) 

This  is  the  second  of  two  collision  frequencies  sought.  The 
results  are  summarized  as  follows: 


Type  of 
Collision 

r- 

Gas 
Press, 
mm  Hg 

Gas 

Composition 

Gas 

Temp. 

sir 

Collision 
Frequency ,  Z 
molecules/cc-sec 

Pressure 
Dependency 
of  Z 

(constant  A«j) 

Gas-Drop 

1.74 

All  NO2 

285 

3.1  x  10*9 

P 

Gas-Gas 

1.74 

/3  -  21 

285 

1.2  x  1023 

p2 

The  frequency  of  the  gas-gas  collisions  is  about  4000  times 
greater  than  the  gas-drop  collisions  at  the  low  pressure  of 
the  calculations.  At  higher  pressures  the  discrepancy  would 
be  greater,  as  seen  from  equations  (41)  and  (38)  for  a  constant 
drop  surface  area.  Of  course,  not  every  collision  between 
oxidizer  and  fuel  (whether  gaseous  or  liquid)  results  in 
reaction.  However,  if  the  activation  energies  and  steric 
factors  are  even  roughly  similar  for  the  two  cases,  the  gas 
phase  reactions  would  be  dominant. 
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b.  Ignition  Delay  Tiae 

The  analytical  expressions  yielding  ignition 
delays  in  an  engine  are  based  on  the  concept  indicated  in 
Figure  1.  In  essence,  an  equation  giving  ignition  delay 
times  is  developed,  the  vapor  concentrations  required  in  the 
equation  being  given  in  terms  of  the  partial  pressures  of 
the  reactants  in  the  thrust  chamber.  The  partial  pressures 
of  each  propellant  are  obtained  by  the  analysis  for  chamber 
pressurization  due  to  propellant  vaporization  (Section  III  B)  . 

Ve  consider  a  volume  element  of  vapor-phase 
reactants  in  a  thrust  chamber  and  assume  heat  losses  from  the 
volume  are  small.  The  heat  equation,  for  a  bimolecular 
reaction  and  no  heat  losses,  becomes  (eq.  36): 

yOg  dT  ”  QW  “  AQClc2  e  -Ea/KT  (44) 


Introducing  non-dimensional  variables 

0  -  ~  €  -  At  J  C]Co  ,  oc  -  Is.  , 

Tg  ,  v  1  *  RTg 


equation  (44)  becomes 


A  .  9  yf  cic2 

SpgT& 


(45) 


In  the  case  where  ot  >  17,  Todes  has  shown  that  the  ignition 

delay  is  very  closely  approximated  by: 

e<* 

ign  -  —  (46) 

<x  X 

Although  in  our  case  the  values  of  OC  are  between  9  and  12  due 
to  the  low  activation  energies  of  the  ignition  reactions,  we 
will  use  equation  (46)  in  the  absence  of  a  better  criterion  of 
ignition.  Reintroducing  the  usual  variables  in  equation  (46), 
one  obtains: 
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r 


ign 


RTg2 

EaAQ 


/°gcpg  e  ®a>^Tg 
clc2 


(47) 


Expressing  the  concentrations  and  C2  in  terms  of  reactant 
partial  pressures  pf  and  pox,  and  expressing  the  density  and 
heat  capacity  of  the  r^/'tant  mixture  in  terms  of  the  properties 
of  the  individual  specie«t  we  obtain  finally  for  the  ignition 
delay  time: 

T  ign  -  r2t83  (1  +  i)  (Cpf  +  /3  CpOJ[  )  e  Ea/RTg  (48) 

PgEaAQ 

where  Pg  -  Pf  +  Pox 

^  -  Pox 
Pf 

and  and  CpQX  are  molar  heat  capacities. 

The  reactant  partial  pressures,  pf  and  pQX  are  obtained 
individually  by  the  analysis  in  Section  III  B  for  chamber 
pressurization  due  to  propellant  vaporization. 

The  theoretical  ignition  delay  time  in  a  given  engine 
operating  under  specified  conditions  is  obtained  by  determining 
the  point  of  intersection  of  the  two  curvbs  of  Figure  1.  From 
the  individual  partial  pressure  curves  for  each  reactant  (by  . 
Sect.  Ill  B) ,  the  total  pressure  curve  for  the  chamber  can 
be  constructed  assuming  no  reaction,  as  well  as  a  vapor  phase 
composition  curve  (  /3  values)  and  an  average  gas  temperature 
curve  for  the  vapor  mixture.  Values  of  Pg,  Tg  and  /9  for 
various  vaporization  times  /are  substituted  into  equation  (48) 
until  the  calculated  Tign  agrees  with  the  time  associated 
with  the  particular  values  of  Pg,  Tg  and  /3  used. 

Obviously,  equation  (48)  applies  to  a  system  in 
which  not  only  the  pressure  and  temperature  are  constant  during 
the  ignition  delay  period  but  also  the  composition  of  the 
reactant  mixture.  Although  the  gas  temperature  variation  is  not 
too  great  during  chamber  pressurization  due  to  propellant 
vaporization  (Fig.  2) ,  the  gas  pressure  varies  considerably  as 
it  rises  toward  a  steady-state  vaporization  chamber  pressure. 
Thus,  the  intersection  of  the  two  curves  of  Figure  1  where  the 
dotted  curve  is  given  by  eq.  (48)  would  be  expected  to  give 
ignition  delays  shorter  than  corresponding  experimentally 
determined  Ignition  delays. 
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In  addition  to  the  pressure  variation,  the  vapor 
composition  can  be  expected  to  vary  due  to  the  different 
vaporization  rates  of  the  two  reactants.  In  general,  the 
vapor  composition,  for  usual  liquid  mixture  ratios,  is  quite 
lean  since  the  oxidizer  is  much  more  volatile.  The  sen¬ 
sitivity  of  calculated  ignition  delays  to  vapor  composition, 
temperature  and  pressure  is  given  in  the  following  section. 

2.  Results  of  Hypergolic  Ignition  Model  and  Comparison 
with  Experiment 

a.  N9O4/UDMH 

Several  pressure-ignition  delay  curves,  cal¬ 
culated  by  equation  (48) ,  are  plotted  in  Figure  25  for  N2O4/ 
UDMH.  Three  /3  values  and  two  temperatures  have  been  used 
in  calculating  the  cury.es.  The  lower  temperature  (286°K) 
is  an  estimated  value  for  the  gas  temperature  in  a  typical 
engine  during  the  ignition  time  based  on  the  chamber 
pressurization  calculations.  The  three  values  of  &  used 
with  this  temperature  indicate  the  effect  of  vapor  composition 
on  ignition  delay.  The  most  favorable  vapor  composition  for 
ignition  is  the  equimolecular  mixture  for  which  /%  -  1 

(Curve  3) .  The  stoichiometric  vapor  mixture  (  /3  -4)  causes 
somewhat  longer  delays,  (Curve  2).  In  a  rocket  chamber  at 
space  conditions,  one  might  expect  that  the  vapor-phase 
composition  (for  simultaneous  propellant  entries)  would  be 
roughly  proportional  to  their  vapor  pressures  at  the  tempera¬ 
ture  in  question.  On  this  basis,  a  /9  of  about  6  for  N2O4/ 
UDMH  gives  the  vapor  composition.  The  ignition  delays  for  this 
case  (Curve  1)  are  seen  to  be  about  40%  longer  than  the  delays 
of  the  equimolecular  mixture  at  the  seme  pressure. 

The  effect  of  temperature  is  indicated  by  a 
comparison  of  Curves  2  and  4  which  are  both  for  the  stoichio¬ 
metric  vapor  mixture.  The  ignition  delays  for  the  lower 
temperature  case  are  roughly  50%  longer  than  the  delays  for 
the  warmer  mixture  at  the  same  pressure.  The  temperature  • 
difference  is  12°C. 

The  pressure-ignition  delay  curve  based  on 
estimated  preignition  vapor  phase  composition  and  temperature 
in  an  engine,  (Curve  1  of  Fig.  25)  is  given  in  Figure  26  together 
with  experimental  ignition  delays  reported  in  Ref.  1.  The  x 
experimental  ignition  delay  times  start  when  the  lagging  v 
propellant  first  emerges  as  a  liquid  from  the  injector  face. 


-100- 


Pressure 


Chamber  Pressure  at  Ignition, 


AFRPL-TR-6 5-2  57 


Figure  26.  Coaiparison  of  Experimental  and  Theoretical  Ignition 
Delays  for  Estimated  Preignition  Vapor  Phase 
Coaposition  and  Temperature 
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These  delays  are  plotted  against  the  pressure  in  the  chamber 
when  ignition  occurred.  The  delays  were  measured  with  engines 
of  various  thrust  chamber  configurations  (design  Pc>  L*,  C.R.) 
operating  at  1  mm  Hg  ambient  pressure. 

The  theoretical  curve,  calculated  for  a  fixed 
gas  temperature  and  vapor  phase  composition,  does  predict 
delays  shorter  than  experimental  ones  as  anticipated.  One 
could  use  some  average  value  of  the  pressure  in  the  chamber 
during  the  ignition  delay  time  in  oruV;  to  get  better  agreement 
bet weep  experimental  and  theoretical  results.  However,  this 
would  be  of  little  value  until  the  actual  vapor  phase  composi¬ 
tion  and  temperature  can  be  predicted  by  the  chamber  pressuriza¬ 
tion  analysis  in  Section  III  B.  As  stated  previously,  an 
accounting  of  the  transient  build-up  to  full  flowrate  for  the 
very  volatile  propellants  (N2O4)  must  be  included  in  the  analysis 
in  order  to  properly  predict  chamber  pressurization,  vapor 
composition  and  temperature.  Since  the  less  volatile  hydrazine 
fuels  undergo  no,  or  only  slight,  "flashing''  inside  the  injector 
volume,  the  analysis  predicts  their  chamber  pressurization  rates 
adequately. 

b.  N2O4/HMH 

Theoretical  pressure-ignition  delay  curves  for 
N2O4/MMH  are  given  in  Figure  27.  The  same  two  temperatures 
are  used  here  as  in  the  N2O4/UDMH  case.  Again  the  equinolecular 
mixture  (  -  1)  gives  the  shortest  delays  at  a  given  tempera¬ 

ture  (Curve  3).  The  stoichiometric  mixture  (  /3  -  2.5)  gives 
only  slightly  longer  delays  (Curve  2).  Using  the  same  criterion 
as  before  for  estimating  the  preignition  vapor  phase  composition, 
a  high  /3  results  due  to  the  lower  vapor  pressure  of  MMH.  The 
calculated  curve  for  this  case  indicates  exceptionally  long 
delays.  The  temperature  effect  on  ignition  delays  for  this 
combination  is  seen  by  comparing  Curves  2  and  4. 

A  comparison  of  Figure  27  for  N204/10IH  with 
Figure  25  for  N2O4/UDMH  shows  the  latter  combination  has  sig¬ 
nificantly  shorter  delays. 

3*  Pressure  Spiking  Consideration 

Two  possible  causes  of  pressure  spiking  are  (1) 
explosion  of  the  "pre-ignition  reaction  product"  (adduct)  at 
ignition  and  (2)  very  fast  reaction  of  accumulated  propellants. 
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a .  Pre-ignition  Reaction  Product 

During  the  course  of  the  experimental  determination 
of  the  overall  kinetic  rate  factors  for  ignition,  it  vas  found 
that  all  four  hydrazine-type  fuel/N2(>4  combinations  formed  a 
reaction  product  under  conditions  at  which  ignition  did  not 
occur.  The  "adduct"  condensed  on  the  walls  of  the  reactor  tubes 
and  generally  collected  in  the  annulus  at  the  base  of  the  tubes. 

The  adduct,  which  formed  at  pressures  below  the  minimum 
ignition  pressures  (Tables  V  through  VIII),  is  a  clear  yellow, 
viscous  liquid  when  it  accumulates  at  the  base  of  the  reactor. 

¥i;th  N2H4  as  fuel,  a  negligible  quantity  of  the 
adduct  collected  in  the  reactor  annulus.  Some  material  did 
condense  on  the  reactor  walls  but  the  droplets  "boiled"  vigorously. 
Particularly  with  N2H4,  readily  observable  smoke  was  entrained 
in  the  flowing  gas  stream  at  the  sub-ignition  pressures. 

The  liquid  adducts  that  collected  with  MMH, 

UDMH,  and  50-50  as  fuels  are  stable  at  room  temperatures  and 
pressures  and  have  very  low  vapor  pressures. 

Chemical  analyses,  which  were  performed  mainly 
on  the  MMH  adduct,  show  the  adduct  has  the  characteristics  of 
a  monopropellant  and  contains  considerable  energy.  Since  the 
adduct  forms  at  pressures  which  exist  during  ignition  delay 
periods  in  attitude  control  engines,  it  is  reasonable  to  suspect 
that  the  adduct  is  the  cause  of  the  pressure  spikes  that  occur 
during  engine  start  transients. 

A  qualitative  correlation  between  the  quantity 
of  adduct  formed  and  tendency  toward  spiking  appears  to  exist. 

It  was  found  in  the  case  of  NO2/MMH  that  the  adduct  which 
collected  at  the  base  of  the  reactor  annulus  was  about  20%  of 
the  total  mass  of  reactants  flowed  during  the  measured  time 
period  (9  minutes).  Qualitatively,  N(>2/50-50  formed  a  similar 
amount  of  adduct  but  NOo/XJDMH  formed  only  about  one-tenth  as 
much.  Engine  tests  indicate  that  spiking  tends  to  be  less 
severe  with  UDMH  as  fuel  than  with  either  MMH  or  50-50  as 
fuel.  The  tendency  of  N2O4/N2H4  toward  spiking  was  not  inves¬ 
tigated  in  this  program. 

It  is  interesting  to  note  that  the  50-50  blend 
of  N2H4  and  UDMH  forms  an  appreciable  amount  of  liquid  adduct 
while  neither  fuel  individually  with  NO2  forms  more  than  about 
10%  as  much . 
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It  should  be  mentioned  that  the  experimental 
apparatus  in  which  the  adducts  vere  formed  did  not  allow  rates 
of  formation  of  the  adducts  to  be  determined.  The  adducts 
formed  ire 11  within  100  msec  but  whether  they  can  form  during 
typical  Ignition  delay  times  (up  to  10  msec  or  so)  cannot  be 
stated  at  present. 

b.  Fast  Reaction  of  Accumulated  Propellants 

A  second  possible  cause  of  pressure  spiking  is, 
of  course,  very  fast  reaction  of  accumulated,  unreacted  propel¬ 
lants.  Once  the  analysis  of  chamber  pressurization  due  to 
propellant  vaporization  is  completed  (Section  III  B.2.h),  the 
mathematical  model  of  hypergolic  Ignition  can  give  the  quantities 
of  fuel  and  oxidizer  in  the  chamber  at  ignition,  the  mass  of 
each  in  the  vapor  phase,  etc.  In  essence,  the  hypergolic 
ignition  model  gives  not  only  the  ignition  delay  time  but  also 
the  conditions  in  the  chamber  at  ignition  from  which  pressure 
spikes  result.  These  are  then  the  initial  conditions  for  a 
pressure  spiking  model  (autoignition/detonation  transition 
concept)  which  would  permit  evaluation  of  chamber  geometry  and 
ignition  hydraulics  from  a  pressure  spiking  viewpoint. 
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IV 

TASK  II  -  EVALUATION  OF  ADDITIVES  TO  REDUCE  ACTIVATION  ENERGY 


(U)  Hie  experimental  apparatus  used  in  the  chemical 
kinetics  study,  Figure  12,  is  well  suited  to  the  screening 
of  additives.  Since  the  fuels  are  metered  in  the  liquid  state 
and  are  subsequently  vaporized  drop-by-drop  before  entering 
the  reaction  chamber,  it  is  assumed  that  the  fuel-additive 
vapor  stream  has  the  same  composition  as  the  liquid  fuel-addif ive 
mixture.  Preferential  vaporization  is  eliminated.  Also, 
pressure  ignition  limits  are  a  sensitive,  reliable  measure 
of  ignition  delays.  The  pressure  limits  are  quite  easily 
measured,  and  are  obtained  much  less  expensively  and  more 
accurately  than  engine  ignition  delays.  Screening  tests  of 
many  additives  can  be  accomplished  in  a  relatively  short 
time.  Promising  additives  should  then  be  checked  in  actual 
engine  tests  to  ascertain  the  improvement  actually  realized. 

(U)  Five  fuel  additives  and  one  oxidizer  additive  were 
tested  in  the  low  pressure,  premixed-vapors  apparatus  to 
evaluate  their  effect  on  hypergolic  ignition  of  NO2/MMH  at 
low  pressures.  Generally,  tests  were  made  for  both  the 
stoichiometric  and  equimolecular  mixtures. 

A.  Fuel  Additives 

(U)  The  five  organic  fuel  additives  are  miscible  with 
MMH  and  so  mixtures  of  the  liquids  were  prepared  and  tested 
in  the  same  manner  as  the  neat  fuels  themselves,  i.e.,  liquid 
phase  flow  metering  followed  by  drop-wise  vaporization  in  a 
heated  glass  spiral. 

(U)  The  results  of  the  tes^s^are  given  in  Table  XIV. 

The  minimum  ignition  pressures  For  the  two  "neat"  mixtures 
were  measured  prior  to  the  additive  testing. 

(U)  Hie  data  of  the  table  shows  that  of  the  five  fuel 
additives,  only  furfuryl  alcohol  had  a  significant  beneficial 
effect.  The  minimum  ignition  pressure  of  the  neat  stoichio¬ 
metric  mixture  was  reduced  by  about  25%  and  that  of  the  neat 
equimolecular  mixture  by  about  17%.  The  other  four  additives 
had  no  significant  effect,  one  way  or  the  other. 
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TABLE  XIV 

EFFECT  OF  FUEL  ADDITIVES  ON  IGNITION  OF  N^/MMH  MIXTURES 


Additive 

Minimum  Ignition 

Pressure,  mm  Hg 

Stoichiometric 

0  -  1.0 

Equimolecui&r 

0  *  2.5 

None 

6.0 

4. 4, 4. 3 

Furfuryl  Alcohol 

4 . 4 , 4 . 6 , 4 . 5 

4. 0,3. 5, 3. 3 

Phenylether 

5. 8, 5. 3, 5. 8 

4. 5,4. 5 

Methylbutynol 

6. 0,6.0 

4. 0,4. 2 

Ethylether 

6.0,6. 5,6.0 

4.2 

Benzene 

6. 0,5. 8, 6.1 

4. 6, 4, 3, 4.0 

4.7  cm  burner  diameter,  ^  *  10.8 

Total  propellant  flow  rate  -  33.6  cc/sec  at  STP 

All  additives  are  in  amount  of  10%  by  weight  of  fuel 


B.  Oxidizer  Additives 

(C)  The  very  short  ig-ition  delays  experienced  with 
Compound  A  and  the  hydrazine-type  fuels  in  the  unconfined 
impinging  stream  tests  of  the  preceeding  program  (Ref,  I) 
suggested  that  delays  of  N2O4/MMH  might  be  shortened  if  a 
fluorine-containing  additive  were  used.  Compound  R,  i.e., 
CF(NF2)o,  looks  very  attractive  on  paper  as  it  contains  seven 
atoms  01  fluorine  per  molecule,  has  roughly  the  same  volatility 
as  N2O4,  and  is  miscible  with  N2O4  giving  a  mixture  that  is 
safe  to  handle  at  concentrations  of  the  additive  to  over  20%. 

(U)  Because  »2°4  is  flowmelered  in  the  vapor  state  in 
the  experimental  system,  it  was  necessary  to  meter  the 
additive  flow  separately  to  ensure  a  known  concentration. 

The  two  vapor  streams  were  mixed  at  the  base  of  the  burner, 
well  before  the  mixing  point  with  the  fuel  vapor. 
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(C)  The  tests  with  Compou;;J  R  produced  reetiiu  unlike 
any  others  in  the  system.  The  minimum  ignition  pressure  was 
not  only  increased,  but  the  good  reproducibility  of  the  tests 
observed  consistently  to  this  point  no  longer  was  obtained. 

(C)  After  a  control  run  with  the  neat  N0Q4/MME  combina¬ 
tion,  two  tests  with  Compound  R,  at  concentrations  of  3  wgt.  % 
N2O4  and  2  wgt.  %  NgO^,  were  made  using  the  stoichiometric  . 
mixture.  Minimum  ignition  pressures  for  the  doped  system  were 
higher  than  for  the  neat  system.  A  repeat  test  of  the  neat 
combination,  after  purging  the  additive  line,  gave  minimum 
ignition  pressure  of  35%  higher  than  the  corresponding  pressure 
measured  just  pH  or  to  the  Compound  R  additive  tests.  The 
pressure  was  even  higher  than  those  obtained  with  the  additive. 
Two  subsequent  check  tuns  with  the  neat  combination,  this  time 
the  equimolecular  mixture,  gave  again  minimum  ignition  pressures 
well  above  the  value  obtained  prior  to  the  additive  tests. 

(C)  Only  after  thoroughly  washing  the  burner  tubes  did 
results  agree  with  the  previously  measured  values  for  the 
neat  system.  Reproducibility  became  good  once  more.  Compound 
R  apparently  contaminated  the  system.  It  was  noted  that  the 
liquid  which  slowly  collected  at  the  base  of  the  burner  tubes 
during  the  Compound  R  tests  was  darker  in  color  than  the  usual 
clear,  yellow  liquid  from  the  neat  N204/MMH  mixtures. 

(U)  The  Compound  R  flow  system  was  well  suited  to  testing 
gaseous  additives  to  the  oxidizer.  Thus  a  few  tests  with  air 
and  with  oxygen  were  made.  Minimum  Ignition  pressures  with 
both  of  these  additives  were  higher  *'  .  1  for  the  system. 

Increase  in  the  minimum  Ignition  pressure  was  found  to 
depend  on  the  amount  of  gaseous  additive  added-  These  two 
additives  apparently  act  simply  as  diluents. 
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V 

TASK  III  -  THRUST  CHAMBER  DESIGN  PARAMETER  STUDY 


A.  Teat  Apparatus  and  Instrumentation 
1.  Experimental  Facilities 

(U)  The  thrust  chamber  design  parameter  studies 
were  performed  in  a  large  v*«uu*  tank  so  that  the  effect  on 
ambient  pressure  due  to  propellant  vaporization  and  combustion 
would  be  negligible.  The  vacuus  tank  is  7  ft.  in  diameter 
and  25  ft.  long  providing  a  volume  of  approximately  1000  cu.  ft. 
It  ie  fabricated  of  stainless  steel  for  use  with  corrosive 
propellants  and  products.  The  tank  hv.v  four  5  in.  diameter 
ports  along  each  side  for  instrumentation  and  observation 
purposes.  Per  the  seblieren  pictures,  3/4  in.  thick  plate 
glass  windows  which  had  suitable  optical  qualities  vere  used. 

(U)  The  vacuum  pump  system  consists  of  a  Roots  E  225 
pump  and  a  1000  blower.  An  adsorption  column  or  trap  for 

removing  any  residual  and  HP  from  the  exhaust  products  is 
installed  between  the  vacuum  chamber  and  pumps.  The  column 
contains  successive  sections  of  activated  sodium  bifluoride, 
soda  lime,  and  drlerite.  The  vacuum  pumps  and  adsorption 
column  are  shown  in  Figure  28. 

(U)  The  durations  of  the  thrust  chamber  tests  were 
approximately  80  ?isec  which  was  sdequate  to  achieve  steady-stare 
chamber  pressure  to  allow  determination  of  the  90%  Pc  time. 

At  flow  rates  corresponding  to  a  design  thrust  level  of  50  lb. , 
the  pressure  rise  in  the  tank  after  a  test  at  a  10  mm  Hg  pressure 
level  was  less  than  0.5  «w  Hw  of  the  ?hcrt  run  duration 

and  large  tank  volume.  In  most  cases  the  short  duration  of 
the  tests  also  prevented  damage  to  the  transparent  chambers, 
permitting  them  to  be  used  for  a  series  of  tests  without  the 
necessity  of  entering  the  tank  and  replacing  the  cylinder. 

fU)  The  oxidizer  and  fuel  propellant  systems  for  the 
liquid  Com<>>ound  A/hydrazine-type  propellants  and  gaseous  F2/H2 
propellants  are  shown  in  Figures  29  and  30,  respectively.  Each 
liquid  system  consisted  cf  a  300  cc  stainless  steel  tank, 
safety  valve,  propellant  solenoid  valve,  and  associated  hand 
valves  and  tubing.  Each  set  was  mounted  on  a  separate  plate 
which  coQinf >lnod  provlsloDf  Cor  nil  D6C680iry  6xt6rml  conncc** 
tlons  such  is  pressurization,  vents  and  loading  so  that  a 
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Figure  28.  Vacuun  Pumps  and  Adso^j^lon  Column  for 

Fluorinated  Oxidizer  Propellan c  Combinations 


Figure  29.  Schematic  of  Conpound  A/Hydrazine  Propellant  Systems 


Figure  30.  Schematic  of  Fluorine/Hydrogen  Propellant  Systems 
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propellant  temperature  conditioning  bath  could  be  installed, 
if  desired.  The  tank  sizes  were  increased  to  1000  cc  for 
the  fluorine  and  2250  cc  for  the  hydrogen  to  minimize  pressure 
decrease  during  the  tests  with  gaseous  propellants. 

Special  solenoid  x'alves  were  precured  to  assure  fast 
response  rates  and  compatibility  with  both  of  the  oxidizers  used 
The  Compound  A  tests  were  made  with  1/4  in.  venturi  valves 
manufactured  by  the  Fox  Valve  Development  Co.,  Inc.,  Hanover, 

N.  J.  These  valves  were  fabricated  from  stainless  steel  and 
employed  a  pintle  providing  metal-to-metal  propellant  shutoff 
in  the  throat  of  the  oxidizer  venturi  and  Tef lon-to-metal 
shutoff  in  the  fuel  valve.  As  a  result,  the  propellant  volume 
in  the  valve  downstream  of  the  seat  was  quite  small  and  the 
valve  was  mounted  directly  to  the  back  face  of  the  injector 
which  minimized  the  time  required  for  the  propellant  to  fill 
the  system.  The  normal  valve  opening  time  was  3-4  msec  at 
24  volts  D.C.  Initially,  attempts  were  made  to  use  Teflon 
poppets  in  the  oxidizer  valve  also  to  ensure  zero  leakage. 
However ,  it  was  found  that  Teflon  could  not  withstand  the  high 
flowrate  conditions,  necessitating  the  use  of  aluminum  alloy 
(6061-T6)  poppets. 

The  gaseous  propellant  valves  were  similar  in  design 
but  were  sized  for  1/2  in.  lines  so  that  flowrates  equivalent 
to  50  lb.  thrust  could  be  obtained  with  gaseous  propellants 
without  excessive  pressure  drop.  Orifices  also  were  used  in 
place  of  the  venturi  sections.  Again,  aluminum  alloy  was  used 
fer  the  oxidizer  poppet  to  ensure  suitable  compatibility. 

Teflon  was  used  for  the  fuel  poppet.  Both  the  1/4  in.  and  the 
1/2  in.  sets  of  propellant  valves  were  powered  by  a  separate 
90  VDC  power  supply  to  decrease  valve  response  time  to  about 
3  msec  with  the  1/2  in.  valves  and  somewhat  less  with  the  1/4  in 
valves . 


,  Standard  1/4  in.  Jamesbury  ball  valves  having  Teflon 

seats  were  used  as  safety  and  vent  valves  for  all  propellant 
systems.  The  Teflon  seats  eroded  in  the  F2  vent  valve  about 
half  way  through  the  F2  tests  and  the  Compound  A  safety  valve 
seats  were  partially  eroded  at  the  end  of  the  test  program 
although  tne  latter  was  not  severe  enough  to  leak  and  may 
have  occurred  during  passivation  with  F2.  The  F2  valve  seats 
were  replaced  and  the  valve  functioned  satisfactorily  for  the 
balance  of  the  F2/H2  tests. 
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2.  Experimental  Hardware 

The  50-lb.  thrust,  attitude  control  chambers  consisted 
of  an  injector-valve  assembly,  a  chamber  flange,  a  transparent 
chamber  section  for  photographic  observation,  and  a  stainless 
steel  nozzle.  A  typical  thrust  chamber  assembly,  including  tne 
Fox  valves  used  for  the  Compound  A  tests  and  two  Kistler  pressure 
transducers  installed  in  the  nozzle,  is  shown  in  Figure  31. 

An  exploded  view  of  the  same  thrust  chamber  is  shown  in  Figure 
32.  The  Kistler  water  cooled  adapters  shown  in  the  pictures 
were  used  only  during  the  initial  portion  of  the  original 
10-month  program.  In  the  present  program,  the  Model  601  and 
603  transducers  were  flush-mounted  in  the  chamber.  A  thin 
layer  of  grease  over  the  sensing  surface  delayed  until  after 
shutdown  the  onset  of  Spurious  signals  due  to  thermal  effects 
from  the  corabusion  gases. 

The  nozzle  throat  diameter  and  chamber  length  and 
diameter  were  varied  to  provide  a  range  of  design  parameters 
including  characteristic  lengths  of  5,  10,  30  and  50  in.  with 
nominal  contraction  ratios  of  1.5,  3.5  and  8  at  design  chamber 
pressures  of  20,  75  and  200  psia.  The  transparent  chamber 
sections  were  fabricated  from  commercially  available  sizes  cf 
acrylic  tubing  which  were  simply  cut  to  the  desired  length. 

Since  standard  tube  diameters  were  used,  actual  contraction 
ratios  ranged  from  1.4  to  9. 

A  separate  stainless  steel  nozzle  section  was 
fabricated  for  each  of  the  three  contraction  ratios  at  design 
chamber  pressures  of  20  and  75  psia  and  for  contraction  ratios 
of  3.5  and  8  at  the  200  psia  design  chamber  pressure.  Tests 
were  not  made  at  a  contraction  ratio  of  1.5  at  200  psia  chamber 
pressure  because  resulting  chamber  diameter  was  too  small  to 
be  used  with  the  desired  injector  configurations.  Tests  also 
were  not  made  with  the  1.5  contraction  ratio  nozzles  at  an  L* 
of  50  in.  at  any  of  the  design  chamber  pressures  because  the  L/D 
of  the  chamber  became  impractically  long.  Also,  the  8  contract¬ 
ion  ratio  at  an  L*  of  5  in.  was  not  physically  possible  to  test 
because  the  volume  in  the  nozzle  section  was  more  than  equi¬ 
valent  to  a  5  in.  I.*. 

Each  nozzle  section  included  a  3/4  in.  length  -with 
the  same  internal  diameter  as  the  inside  of  the  respective 
chamber  section  to  permit  mounting  of  the  Kistler  pressure  trans¬ 
ducers.  The  volume  of  this  section  as  well  as  that  of  the  con¬ 
vergent  portion  of  the  nozzle  was  included  in  the  I*  calculations 
for  sizing  the  transparent  sections.  Typical  nozzles  are  shown 
in  Figures  33  and  M.  The  nozzles  shown  in  Figure  33  are  designed 
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5801-13 

Figure  33.  Set  of  Nozzles  for  Same  Design  Chamber  Pres¬ 

sure  (75  psi)  and  Different  Contraction  Ratios  (1.5-8) 


5801-16 

Figure  34.  Set  of  Nozzles  for  Same  Contraction  Ratio 

(3.5/1)  and  Different  Design  Chamber  Pressures  (20-200  psi) 
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for  the  UM  nominal  chaaber  pressure  (75  psia)  but  for 
different  contraction  ratios.  The  nozzles  in  Figure  34  are 
designed  for  tbe  same  contraction  ratio  (3,5) ,  but  for 
different  nominal  chaaber  pressures. 

The  entire  chaaber  asseably  was  assembled  with 
1/4  in.  threaded  rods.  Buna-N  O-rings  were  used  to  seal  the 
flanges  and  tbe  transparent  chaaber  sections  with  satisfactory 
results.  This  type  of  construction  and  asseably  was  used  to 
facilitate  rapid  changes  of  cbaabers  and  nozzles  because  of 
the  large  number  of  configurations  to  be  tested.  In  some 
tests  the  chaaber  flange  at  the  injector  end  was  replaced  with 
a  thicker  flange  to  allow  room  for  installation  of  a  Kistler 
pressure  transducer  at  the  injector  end  of  the  chamber  as  well 
as  the  sozsle  end.  In  these  tests  the  transparent  chamber 
length  was  reduced  to  maintain  the  design  1>*. 

The  twenty-seven  thrust  chamber  configurations  tested 
in  the  Task  III  program  are  shown  in  Table  XV. 

The  two  injectors  used  in  the  liquid  propellant  tea's 
are  shown  in  Figure  35.  The  majority  of  the  tests  were  mad'~ 
with  the  single  element  doublet  injector  having  0.052  in. 
diameter  orifices.  The  impingement  angle  was  60°,  the  impinge¬ 
ment  length  was  0.12  in,,  oxidizer  velocity  was  15  ft/sec  and 
fuel  velocity  80  to  100  ft/sec.  The  single  element  triplet 
Injector  maintained  the  60°  angle  between  the  center  oxidizer 
orifice  and  each  outside  orifice  end  the  impingement  length 
was  also  0.12  in.  with  the  same  propellant  velocities.  In 
both  designs,  the  Fox  valves  were  mounted  directly  on  the  back 
of  the  injector  with  the  outlet  section  of  the  venturi  inside 
the  injector  body  ao  the  throat  was  approximately  in-line  with 
the  back  face  of  the  injector.  This  minimised  the  volume  be¬ 
tween  the  valve  poppet  at  the  venturi  throat  and  the  injector 
orifices. 


The  injector  for  the  gaseous  propellant  tests,  shown 
in  Figure  35,  consisted  of  a  single  axial  orifice  for  the 
fluorine  surrounded  by  a  concentric  annulus  for  the  hydrogen. 
The  orifice  wan  0.20  in.  in  diameter  and  the  annulus,  angled 
at  30®  to  the  a are  so  tbit  the  hydrogen  shee*  impinged  or  the 
fluorine  Stream,  wet  0.020  in.  wide.  Tbe  larger  propellant 
valves  Counted  directly  on  the  rear  fac*  of  the  injector  to 
mlniialm4  volume  downstream  of  the  valve. 
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Oxidizer 

Venturi 

Valve 

Fuel 

Venturi 

Valve 

Bggjgg^jggg 

Both  Orifices  =  0.052  in.  diameter 


(a)  Doublet  Injector  with  Valves  Assembled 


Fuel  Inlet 


Oxidizer  Orifice  ■  0.052  in.  diameter 
Fuel  Orifice  -  0.033  in.  diameter 
(Two  sections  taken  90°  apart) 

(b)  Triplet  Injector 

Figure  35.  Injectors  for  Compound  A/N2H4  -  Type  Propellants 
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tteorlne  Orifice  *  0.209  ia.  diameter 
Bydfogea  Aaaalae  -  0.020  ia.  vide 
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3.  Instrumentation 


In  addition  to  conventional  instrumentation  required 
for  propellant  pressurization  and  system  monitoring,  special 
instrumentation  was  employed  to  determine  ignition  delays  and 
monitor  pressure  transients  during  the  ignition  process.  This 
instrumentation  included  high  speed  schlieren  movies  of  the 
propellant  impingement  zone,  photomultiplier  tube  to  detect 
ignition,  high  speed  direct  photography,  and  high  response 
pressure  transducer  and  recording  equipment. 

a.  Schlieren  System 

The  purpose  of  the  schlieren  system  was  to  investigate 
propellant  stream  characteristics  and  determine  the  time  at 
which  the  propellants  impinged.  Not  only  can  the  initial  time 
of  contact  be  ascertained  by  this  method,  but* information  on 
vaporization  and  mixing  of  the  propellants  could  be  obtained. 
Since  it  was  anticipated  that  initial  injection  of  the  propel¬ 
lants  would  be  in  the  vapor  phase  at  the  ambient  pressures 
below  10  mm  Hg,  direct  photography  alone  was  not  adequate  for 
this  purpose. 

Although  the  schlieren  instrumentation  could  not 
be  used  to  detect  propellant  entry  during  the  actual  thrust 
chamber  tests  because  of  the  poor  optical  properties  of  the 
plastic  chambers,  it  was  used  to  study  the  Injection  charac¬ 
teristics  of  each  injector  with  each  of  the  propellants  that 
were  investigated.  Schlieren  films  were  made  of  each  injector 
without  the  chamber  installed  in  order  to  determine  the  time 
from  valve  signal  to  propellant  entry  as  both  vapor  and  liquid 
to  aid  in  analysis  of  the  thrust  chamber  ignition  tests. 
Schlieren  films  were  also  made  in  which  both  propellants  were 
injected  simultaneously  in  order  to  determine  ignition  delay 
in  these  unconfined  tests  for  comparison  with  the  thrust 
chamber  tests. 

A  two-mirror,  parallel-path  schlieren  system  was 
used  with  a  Fastax  high-speed  camera  to  obtain  suitable  time 
resolutions.  The  essential  components  of  the  system  included 
a  Unertl  Model  BH6  Normal  and  Color  Schlieren  Source,  a  pair 
of  eight-inch  front-surface  parabolic  mirrors  of  64  in.  focal 
length,  knife-edge,  and  the  Fastax  camera  capable  of  up  to 
approximately  16,000  pictures  per  second  using  split-frame 
optics.  Time  resolution  at  ignition  was  approximately  7  frames 
(14  pictures)  per  msec. 
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Wjw#  37.  SokMltic  of  Thrust  Clumber  Test  oetup 
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The  schlieren  light  source  and  one  parabolic 
■irror  ar£  shown  in' Figure  38.  The  second  mirror  and  Fastax 
camera  with  its  controls  are  shown  in  Figure  39. 

b.  Direct  Photography 

A  Kodak  camera  having  a  speed  of  approximately 
3000  frames  per  second  and  a  second  full-frame  Fastax  camera 
were  used  for  direct  photography  in  selected  tests  with  the 
complete  thrust  chamber  assembly.  They  were  used  to  determine 
the  location  of  ignition  in  the  chamber  and  the  nature  of  the 
propagation  of  the  flame.  The  cameras  viewed  the  chamber  through 
tbs  same  tank  windows  as  the  schlieren  and  were  located  about 
six  feet  from  the  combustion  area.  Tank  lights  were  not  on  so 
ignition  location  could  be  detected  from  the  combustion  light 
and  therefore  propellants  could  not  be  seen  before  ignition  in 
these  tests. 


c.  Ignition  Detection 

To  determine  the  time  at  which  ignition  occurred, 
a  flame  detector  consisting  of  an  RCA  IP  28  photomultiplier  tube 
was  used.  The  tube  is  sensitive  to  wavelengths  from  about  2200A 
to  6000A.  No  filter  was  required  since  the  tank  was  not  illumi¬ 
nated  for  the  thrust  chamber  tests.  The  flame  detector  was 
located  in  the  vacuum  chamber  approximately  6  feet  from  the 
injector  with  an  unobstructed  view  from  the  injector  to  the 
bottom  of  the  tank,  *  distance  of  about  3  1/2  feet.  It  could 
sense  ignition  through  the  plastic  chambers  in  the  thrust 
chamber  tests  or  elsewhere  in  the  tank  in  the  unconfined  tests 
regardless  of  origin. 

d.  Pressure  Instrumentation 

Chamber  pressure  was  measured  with  Model  601  and 
603  Kistler  pressure  transducers  and  model  566  Kistier  charge 
amplifiers  The  outputs  were  read  out  on  an  oscilloscope  and 
photographed  for  analysis.  Two  transducers  were  mounted  in  the 
same  plane  at  the  noxzle  end  of  the  thrust  chamber,  as  shown 
in  Figure  31,  for  all  tests.  A  model  601  transducer  having  a 
frequency  of  150,000  cps  and  rise  time  ->f  3  microseconds  was 
calibrated  for  345  am  Bg  pressure  at  3  cm  deflection  to  define 
the  preignition  chamber  pressure  rise  due  to  propellant  entry 
and  vaporisation.  A  Model  603  transducer  with  only  1/5  the 
sensitivity,  but  with  a  natural  frequency  of  200,000  cps  and 
rise  time  of  1  microsecond,  was  calibrated  for  several  chamber 
pressure  ranges,  depending  upi  the  design  chamber  pressure,  to 
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5801-6 

Figure  38.  Schlleren  Light  Source  end  Oscilloscope 

Setup 


5801-4 


Figure  39.  Schlieren  Csaera  Setup 
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record  steady-state  pressure  as  veil  as  any  possible  pressure 
spikes  which  night  be  experienced.  A  second  Model  603  trans¬ 
ducer  was  used  at  the  injector  end  of.  the  chamber  in  some 
tests  to  compare  the  chamber  pressure  at  both  ends  of  the 
chamber.  The  transducers  were  uncooled  because  of  the  short 
duration  tests,  but  the  diaphragms  were  protected  from  over¬ 
heating  with  a  layer  of  silicone  grease  to  prevent  drift  after 
ignition.  Although  the  grease  would  melt  after  several  tests, 
particularly  at  the  highest  chamber  pressure,  and  have  to  be 
replaced,  it  did  not  appear  to  react  with  the  propellaats  nor 
did  it  affect  the  response  of  the  transducers  as  indicated  by 
comparative  tests. 

e.  Recording  Equipment 

In  addition  to  th i  schlieren  and  direct  photo¬ 
graphy  films,  the  primary  means  of  data  acquisition  was  a 
Tektronix  Type  551  Dualr-Beam  Oscilloscope  with  suitable  ampli¬ 
fiers  to  permit  simultaneous  display  of  four  parameters.  The 
output  from  the  photomultiplier  tube  and  the  two  Kistler  trans¬ 
ducers  at  the  nozzle  end  of  the  chamber  were  monitored  on  the 
oscilloscope  and  recorded  by  a  Polaroid  camera.  The  fourth 
channel  was  used  for  the  third  Kistler  transducer  or  for  the 
current  from  the  propellant  valve  circuit,  depending  upon  the 
nature  of  the  test. 

The  oscilloscope  and  camera  are  visible  in 
Figure  38.  Also  shown  is  the  single-flash  strobe  light  attached 
to  the  parabolic  mirror  pedestal.  The  flash,  having  a  duration 
of  about  2  msec,  was  detected  by  the  cameras  and  by  the  photo¬ 
multiplier  tube  and  was  used  as  a  reference  signal  to  correlate 
the  various  instrumentation. 

B.  Experimental  Thrust  Chamber  Program 

The  objective  of  the  Task  III  thrust  chamber  tests  was  to 
evaluate  the  influence  of  thrust  chamber  design  parameters  on 
the  ignition  delay  of  two  advanced  propellant  combinations, 
gaseous  fluorine/hydrogen  and  Compound  A/hydrazine-based  fuels. 
The  specific  design  parameters  of  interest  were  characteristic 
length,  chamber  pressure 9  contraction  ratio,  and  oxidizer  lead. 
The  programmed  range  of  these  parameters  is  shown  below: 

Characteristic  Length  5-50  in. 

Chamber  Pressure  20-200  psia 

Contraction  Ratio  1.5-8 
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Oxidizer  Lead 


+.002  to  -.002  sec 


Thrust 


50  lb 


Ambient  Pressure 


0.2  psia  or  less 


(U)  The  principal  propellant  combinations  investigated 
during  the  program  were  gaseous  fluorine/hydrogen  and  Compound  A/ 
N0H4.  For  comparison,  tests  were  also  made  with  Compound  A/ 

UDMH  and  Compound  A/MHF-5. 

(U)  Since  unconfined  impingement  tests,  which  have  been 
reported  in  AFRPL-TR-65-105 ,  indicated  that  injection  para¬ 
meters  had  no  significant  influence  on  ignition  characteristics, 
simplified  injectors  were  used  so  that  injection  characteristics 
(vaporizing,  mixing)  of  the  various  propellants  could  be  readily 
determined  photographically.  No  attempt  was  made  to  optimize 
the  Injectors  from  a  performance  standpoint.  The  majority  of 
the  teets  were  made  with  a  single  concentric  injector  for  the 
F2/B2  tests  and  a  single-element  doublet  for  the  Compound  A 
tests. 


(U)  All  tests  were  made  with  ambient  pressure  in  the  large 
vacuum  tank  between  7  and  10  mm  Hg  (0.13  to  0.19  psia).  The 
twenty-seven  thrust  chamber  configurations  shown  in  Table  XV 
were  tested  with  the  two  principal  propellant  combinations 
(P2/H2,  Compound  A/N2H4)  to  investigate  the  effects  of 
characteristic  length,  chamber  pressure,  and  contraction  ratio 
with  simultaneous  propellant  injection.  Duplicate  tests  were 
made  with  most  of  the  configurations  to  determine  reproducibility 
of  the  results.  Selected  combinations  of  thrust  chamber  con¬ 
figurations  were  then  used  to  evaluate  the  ignition  charac¬ 
teristics  of  the  other  propellant  cdmbinations,  the  effect  of 
oxidizer  lead  or  lag  and  the  effect  of  mixture  ratio  with  the 
Compound  A  combinations. 

(U)  Results  of  the  experimental  program  are  discussed  in 
the  following  sections. 

C.  Experimental  Results  -  Compound  A/Hydrazine-Based  Fuel 

Tests 

(C)  A  complete  series  of  ignition  tests  was  made  with 
Compound  A/hydrazine  for  the  range  of  thrust  chamber  configura¬ 
tions  which  are  tabulated  in  Table  XV.  The  design  chamber 
pressures  selected  for  these  configurations  were  20,  75,  and 
200  psia  for  a  nominal  thrust  level  of  50  lbs.  Comparative 
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teats  were  unde  with  Goapenms  A/TJDKH  and  Co^ound  A/MKF-5  at 
all  three  chamber  prepares  In  chambers  with  a  30  in.  L*  and 

3.5  contraction  ratio. 

1.  Schlieren  Characterisation  of  Propellants 

(C)  In  order  to  determine  the  ignition  delays  for 
the  war  ion#  pn&f*ellant  combin'*  Hows  and  chanber  configurations , 
the  propel  taint  tmtry  times  were  measured  with  the  schlierea 
system .  Since  the  sebliereu  techniques  are  ineffective  with 
the  plastic  chamber©  instaiJ.ec.,  each  propellant  and  injector 
no  evaluated  wi-neat  the  chamber  attached  and  the  results  are 
a bo  vs  in  table  XVI.  £*cb  propellant  was  injected  individually 
to  determine  the  time  at  which  vapor  and  liquid  injection  first 
appeared  without  being  masked  by  the  second  propellant  or 
ignition.  Several  additional  schlieren  tests  were  made  with 
both  propellants  injected  simultaneously  to  determine  ignition 
delay  and  ignition  location  in  the  unconfined  state.  The 
column  lf.be, led  "resistance"  in  'Table  XVI  is  the  value  of  the 
series  resistor  used  in  the  valve  circuit  to  delay  the  second 
propellant  valve,  In  all  schlieren  tests,  s  dummy  propellant 
reive  was  used  in  place  of  the  first  propellant  valve  to 
duplicate  the  electrical  conditions  in  the  actual  chamber 
tests. 


(C)  The  physical  properties  of  the  propellants 
greatly  effect  the  abate  in  which  the  various  propellants  are 
first  injected,  particularly  at  the  low  ambient  pressure 
involved.  Significant  differea-es  between  the  propellants  were 
detect**,  ia  tne  schlierea  films.  For  example,  hydrazine  does 
not  vsipc-riSv  sufficiently  for  the  schlieren  system  to  detect 
any  gameens  propellant  before  entry  of  the  liquid  stream, 
lydrarinc  first  enters  as  ?.  thin  coherent  stream  even  at  low 
ambient  pres jure  with  no  indication  of  vaporisation  or  break¬ 
up  within  the  2  inch  length  visible  in  the  schlieren  field 
of  view  as  shown  in  Figure  40,  Entry  time  of  the  liquid  was 
very  easy  to  detect ,  but  the  entry  time  varied  over  a  range 
of  several  milliseconds  at  low  ambient  pressure.  Three  tests 
*lth  nominal  flow*  of  0.06  Ih/sec  at  20  mm  Hg  pressure  produced 
entry  times  of  16.7,  S.O  and  1C. 4  asec.  These  tests  were  longer 
than  anticipated  so  additional  tests  were  maoe  at  higher  ambient 
pressures  to  determine  if  pressure  influenced  the  timing.  The 
flowrate  remained  constant  with  the  changing  back  pressure  sinr« 
ii  was  controlled  by  the  pressure  upstream  of  the  ve*»C  uri  which 
remained  constant.  Entry  times  at  the  higher  ambient  pressures 
i-  Table  XVI  were  shorter  and  more  reproducible  probably  due  to 
less  vaporisation  of  the  hydrazine. 
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TABLE  XTI 


SCHLIEREN  DATA  -  UNCOHFINED  TESTS  OF  COMPOUND  A,  N2H4,  UDMH,- 
AND  MHF-3^  USING  DOUBLET  INJECTOR 


Propellant 

1 - 

Ambient 

Pressure 

mm  He 

Resistance 

ohms 

* 

lb/sec 

Valve  Signal  to 

Vapor 

msec 

Vapor  ft 
Liquid 

msec 

Liquid 

■sec 

Compound  A 

n 

0 

.120 

2.7 

6.4 

10.2 

20 

0 

3.2 

7.3 

20 

0 

2.9 

7.3 

10.8 

20 

10 

I  .120 

4.8 

8.2 

20 

20 

4.9 

8.2 

n.8 

20 

30 

■ 

4.4 

8.7 

12.0 

20 

0 

KB 

- - 

e .  8 

■  1  M 

20 

0 

.108 

— 

6.7 

10.2 

H2H4 

750 

0 

6.4 

750 

0 

.060 

6.4 

750 

0 

7.3 

750 

0 

.060 

6.4 

150 

0 

.060 

8.6 

100 

0 

.000 

8.2 

50 

0 

.060 

10.0 

20 

*  0 

.060 

9.0 

20 

0 

.060 

10.4 

20 

0 

iiifHALi'fli 

10.7 

20 

5 

.060 

11.2 

20 

10 

.060 

11.0 

20 

20 

.060 

11.7 

20 

30 

.060 

13.5 

20 

40 

15.9 

20 

0 

.072 

10.1 

20 

0 

.052 

11.2 

OMH 

17 

0 

.060 

mm 

■ 

17 

0 

.060 

S^FWill 

17 

0 

.060 

2.9 

Mill 

KHF-5  t 

11 

0 

.060 

7.1 

8.0 

11 

0 

.060 

9.1 

Era mi 

11 

0 

.060 

7.8 

9.8 

11 

30 

.060 

11.7 

13.9 

11 

40 

.060 

14.2 

16.3 

11 

0 

.072 

8.0 

9.2 

11 

0 

.052 

9.3 

12.3 
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5801-22 

Figure  40.  Schlieren  of  Hydrazine  Injection, 

Ambient  Pressure  -  20  mm  Hg 


5801-21 

Figure  41.  Schlieren  of  Compound  A  Injection, 

Ambient  Pressure  ■  20  mm  Hg 
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Varying  hydrazine  flowrate  to  obtain  other  mixture 
ratios  did  not  change  entry  tine  by  more  than  1/2  msec  at  the 
20  mm  Hg  pressure.  Entry  time  was  shorter  for  low  0/F  and 
longer  for  high  0/F  as  would  be  expected  for  the  varying  flow- 
rate.  Fuel  entry  could  be  delayed  up  to  5  msec  by  varying 
the  resistor  in  the  propellant  valve  circuit. 

In  the  ether  extreme,  Compound  A  vaporized  very 
drastically  at  low  pressure  and  the  precise  time  of  Compound  A 
entry  was  more  difficult  to  define.  At  20  mm  Hg  ambient* 
pressure,  vapor  could  barely  be  detected  about  3  msec  after 
valve  "on"  signal  which  corresponds  to  expected  time  of  valve 
opening.  A  mixture  of  vapor  and  liquid  appeared  at  about  7  msec 
and  a  definite  change  to  a  more  dense  liquid  phase  could  be 
seen  between  10  and  11  msec.  In  contrast  with  the  cohesive 
hydrazine  stream,  however,  the  steady-state  oxidizer  flow 
appears  to  vaporize  immediately  upon  injection  at  low  pressure 
and  the  stream  expands  very  rapidly  as  shown  in  Figure  41. 

The  rapid  spread  of  the  oxidizer  flow  obscures  the  entrance  of 
the  fuel  in  the  ignition  tests. 


The  schlieren  tests  with  UDMH,  which  has  a  higher 
vapor  pressure  than  hydrazine,  appear  to  be  similar  to  the 
Compound  A  tests.  Vapors  first  are  visible  at  about  3  msec 
after  valve  "on"  signal  and  they  appear  to  become  a  mixture 
of  liquid  and  vapor  at  7  msec.  However,  there  was  no  clear- 
cut  transition  to  all-liquid  as  with  the  Compound  A. 

The  initial  state  of  MHF-5  appears  to  be  a  mixture 
of  vapor  and  liquid  based  on  the  density  of  the  schliere  and 
the  time  at  which  it  appears  (about  7  msec)  which  corresponds 
quite  closely  with  the  appearance  of  the  UDMH  mixture  at  the 
same  flowrate.  The  MHF-5  also  vaporized  rapidly  after  injection 
when  steady  state  flow  was  established.  The  steady-state 
appearance  of  both  UDMH  and  MHF-5  is  quite  similar  to  that  of 
Compound  A  shown  in  Figure  41 . 

The  schlieren  tests  with  both  propellants  are  reported 
in  subsequent  sections  describing  the  results  of  the  chamber 
configuration  tests  for  the  various  fuels  with  Compound  A.  In 
general,  ignition  times  of  the  unconfined  tests  were  comparable 
to  those  for  the  chamber  tests. 


i 

r 


i 


r 
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2.  Thrust  Chamber  Tests  -  Compound  A/N2H4 

Table  XVII  summarizes  the  results  of  the  thrust 
chamber  tests  with  Compound  A/N2H4  for  the  various  thrust 
chamber  configuration  parameters  with  the  doublet  Injector. 

These  tests  were  made  at  a  mixture  ratio  of  2.0. 

The  significance  of  the  various  data  columns  in 
Table  XVII  and  the  subsequent  tables  is  as  follows.  The 
ambient  pressure,  Pa,  is  the  pressure  in  the  vacuum  tank. 

It  was  measured  before  the  test,  but  it  did  not  increase  more 
than  0.5  mm  Hg  during  a  test.  Chamber  pressure,  Pc,  is  the 
steady-state  value  and  90%  Pc  is  the  time  from  valve  "on" 
signal  to  90%  of  the  steady-state  chamber  pressure.  These 
values  were  determined  from  the  high-range  pressure  transducer 
which  was  also  used  to  determine  the  shape  of  the  pressure 
transient  after  ignition  and  the  magnitude  of  any  pressure 
peaks  which  might  occur.  The  time  from  electrical  signal  to 
the  propellant  valycs  to  both  the  first  indication  of  a  signal 
on  the  photomultiplier  tube  and  the  vertical  deflection  of  the 
photomultiplier  signal  are  designated  by  "PMpirst"  and  "PMVert", 
respect ivelj'.  The  Pign  or  pressures  at  ignition  corresponding1 
to  the  preceeding  photomultiplier  signals  are  also  included 
in  the  table.  These  pressures  were  measured  by  the  pressure 
transducer  having  the  expanded  low  pressure  range  for  detection 
of  initial  pressure  rise  due  to  propellant  vaporization.  The 
final  column  in  the  table  is  the  time  at  which  the  low-range 
pressure  transducer  rose  nearly  vertically,  indicating  a  rapid 
chamber  pressure  rise  due  to  ignition.  The  correlation  between 
this  indication  of  ignition  and  that  of  the  photomultiplier 
tube  is  discussed  further  later  in  this  section. 

Tables  XVIII  and  XIX  summarize  the  additional  mixture 
ratio  and  propellant  lead  tests  made  with  Compound  A/N2H4  using 
the  doublet  injector.  In  Table  XVIII,  the  results  are  shown 
for  tests  at  mixture  ratios  of  1.5,  2.0,  and  2.5.  These  tests 
were  made  with  chambers  having  a  contraction  ratio  of  3.5  and 
L*  extremes  of  5  and  50  in. ,  at  all  three  chamber  pressures. 

Also  included  are  four  unconfined  tests  which  were  monitored 
with  both  the  oscilloscope  and  the  schlieren.  In  all  the 
mixture  ratio  tests  except  test  963,  the  individual  propellant 
flow  rates  were  adjusted  to  maintain  a  constant  total  flow  rate 
of  0.18  lb/sec.  In  test  963,  the  oxidizer  flow  rate  was  reduced 
to  extend  the  low  mixture  ratio  range  to  1.0.  In  this  test  the 
total  flowrate  was  reduced  from  0.18  lb/sec  to  0.144  lb/sec. 

The  tests  at  a  mixture  ratio  of  2.0  are  repeated  from  Table  XVII 
to  provide  a  direct  comparison  with  the  other  tests. 
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Table  XIX  contains  the  oxidizer  lead/lag  tests. 

These  tests  vere  made  by  delaying  the  lagging  propellant  valve 
by  adjusting  the  appropriate  resistor  in  the  valve  circuit. 

The  lead  time  was  computed  from  the  previously  discussed 
schlieren  tests  (Table  XVI)  of  the  individual  propellants. 

Most  of  these  tests  were  performed  at  the  lowest  design 
chamber  pressure  to  more  readily  detect  differences  in 
ignition  time. 


a .  Ignition  Characteristics 


Although  the  photomultiplier  tube  is  used  to 
detect  ignition  in  the  thrust  chamber  tests,  ignition  also  is 
Indicated  in  each  test  by  a  sharp  increase  In  the  rate  of 
pressure  rise  measured  by  the  Klstler  pressure  transducer 
calibrated  for  the  low  pressure  range  to  detect  vaporization 
and  ignition.  Generally,  both  of  these  indications  of  Ignition 
coincide  although  there  may  be  up  to  a  0.75  msec  lag  in  the 
ignition  time  indicated  by  the  pressure  transducer  installed 
in  the  nozzle  section.  This  lag  is  due  to  the  finite  time 
required  for  propagation  of  the  flame  front  from  the  injector 
to  the  nozzle  and  is;  discussed  in  more  detail  In  the  discussion 
of  the  Compound  A/MHF-5  tests. 

In  some  tests,  however,  there  were  slight 
fluctuations  in  the  photomultiplier  trace  on  the  oscilloscope 
up  to  3-4  msec  in  advance  of  the  vertical  deflection  associated 
with  ignition.  These  fluctuations  were  small  in  amplitude  and 
even  returned  to  zero  deflection  for  one  or  more  milliseconds 
in  some  tests  before  the  final  strong  deflection  occurred. 

Both  the  time  from  valve  signal  to  these  initial  small  fluctua¬ 
tions  (when  they  occurred)  and  the  vertical  deflection  are 
reported  under  the  "PMfirst"  and  "PMyert"  columns,  respectively, 
in  Tables  XVII,  XVIII  and  XIX.  The  pressure  in  the  chamber 
at  each  of  these  times  is  also  reported.  In  the  tests  in  which 
the  small  fluctuations  did  not  occur  the  time  indicated  in  the 
"PMfirst”  column  is  the  same  as  in  the  ”PMvert”  column. 

Schlieren  tests  were  made  without  the  chamber 
and  nozzle  installed  to  determine  ignition  delay  in  an  uncon¬ 
fined  condition  and  to  investigate  these  fluctuations  to 
determine  which  signal  from  the  photomultiplier  tube  should 
be  considered  as  the  indication  of  ignition.  These  tests 
are  reported  in  Table  XVIII  as  Tests  928-931  and  are  summarized 
below: 
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The  time  resolution  of  the  schlieren  system  was 
approximately  0.07  msec  per  picture.  ’  As  indicated  in  the  table 
the  schlieren  indication  of  ignition  and  the  vertical  deflec- 
tloc  of  the  photomultiplier  tube  agreed  very  closely  in  each 
of  the  tests.  In  test  930,  however,  the  photomultiplier 
detected  slight  light  output  at  9.4  msec  and  again  at  9.7  msec 
but  returned  to  zero  before  the  strong  deflection  at  12.0  msec. 
The  schlieren  pictures  indicate  that  liquid  oxidizer  entered 
at  9.2  msec  but  that  there  was  no  evidence  of  ignition  until 
12.1  msec.  Ignition  occurred  at  the  impingement  point  of  the 
two  propellants,  well  within  the  field  of  view  of  the  camera. 
Similar  photomultiplier  fluctuations  were  experienced  in 
Test  931  although  the  tine  between  the  initial  fluctuations 
and  the  vertical  rise  was  much  shorter. 

On  the  basis  of  these  tents  it  appears  that  the 
small  fluctuations  in  the  photomultiplier  trace  are  the  result 
of  very  small  localized  and  discontinuous  reactions  and  should 
not  be  considered  as  a  true  indicat 'an  o*  ignition.  This  is 
also  indicated  by  the  relatively  small  increase  in  chamber 
pressure  during  the  period  in  which  these  fluctuations  wccur 
prior  to  the  vertical  deflection  of  the  photomultiplier  tube 
output . 

Ignition  in  a  chamber  was  also  studied  photo¬ 
graphically  with  high-speed  motion  pictures  at  about  4000 
frames  per  second.  In  testa  947  to  950,  reported  in  Table 
X7III,  the  firnt  visible  sign  of  ignition  was  a  small  flame, 

1  to  2  inches  long,  at  the  injector  end  of  the  chamber.  In 
the  two  lower  mixture  ratio  tests  (047-948) ,  the  first  visible 
light  was  detected  at  9  and  10  msec  after  valve  opening  while 
chamber  pressure  was  still  at  the  original  ambient  pressure 
of  10  mm  Eg.  The  flame  filled  moat  of  the  chamber  cross- 
sectional  area  and  required  11  to  13  frames  (about  3  msec) 
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to  reach  the  nozzle  end  of  the  chamber.  During  this  period, 
the  chamber  pressure  increased  to  300  mm  Eg  and  continued 
to  rise  slowly  and  smoothly  thereafter. 

By  contrast,  the  higher  mixture  ratio  tests 
(949-950)  did  not  show  any  visible  light  until  12.5  and  13.5  msec 
by  which  time  the  chamber  pressure  had  increased  to  55  and  70 
mm  Hg  due  to  propellant  vaporization  and  possibly  some  low 
order  reaction.  The  flame  initially  was  visible  only  in  the 
center  of  the  chamber  and  reached  the  nozzle  within  1  msec. 
Chamber  pressure  increased  rapidly  during  the  1  msec  period 
and  the  low  pressure  trace  rapidly  rose  off  scale  during 
this  period.  The  high  pressure  trace  also  rose  abruptly  with 
a  series  of  small  oscillations  of  about  10-15  psi  amplitude. 

In  all  four  tests,  the  photomultiplier  indication  of  ignition 
preceeded  the  abrupt  pressure  rise  by  0.5  to  1.5  msec. 

From  analysis  of  the  motion  pictures,  it  is 
evident  that  ignition  is  initiated  at  the  impingement  point 
of  the  liquid  propellants  and  proceeds  to  fill  the  chamber 
over  a  period  of  1  or  more  milliseconds.  The  time  difference 
between  chamber  pressure  rise  and  photomultiplier  output  is 
due  to  the  finite  time  required  for  ignition  to  travel  the 
length  of  the  chamber.  The  time  delay  will  be  discussed  in 
more  detail  in  the  Compound  A/MHF-5  tests  where  pressure 
transducers  were  locat'd  at  both  ends  of  the  chamber.  As  the 
result  of  the  schlieren  and  direct  photography  tests,  it  was 
concluded  that  the  vertical  deflection  of  the  photomultiplier 
tube  should  be  taken  as  the  true  indication  of  ignition  and 
subsequent  discussions  will  be  based  on  this  time. 


As  discussed  at  the  beginning  of  this  section, 
the  schlieren  tests  indicated  some  variation  in  the  determina¬ 
tion  of  propellant  entry  time.  Both  Compound  A  and  N2H4  entered 
as  a  liquid  between  9  and  11  msec  at  nominal  flowrates.  The 
ignition  times  from  valve  signal  in  Table  XVII  range  from  9  to 
14.5  msec.  Thus,  the  uncertainty  in  propellant  entry  time 
represents  a  large  fraction  of  any  calculated  ignition  delay. 

In  any  event,  ignition  delays  are  certainly  less  than  5  msec 
from  entry  of  the  second  liquid  propellant  and  are  more 
probably  less  than  2  msec  even  at  the  low  ambient  pressures 
used  for  these  tests. 


Effect  of  Thrust  Chamber  Configuration  on 
Ignition  and  Rise  Time 


The  effects  of  thrust  chamber  configuration  on 
ignition  and  rise  time  to  90%  of  steady-state  chamber  are  shown 
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in  Figures  42  and  43,  respectively,  for  the  tests  reported 
in  Table  XVII.  The  times  shown  in  Figure  42  are  from  elec¬ 
trical  signal  to  the  valve  so  that  the  times  required  for 
propellants  to  reach,  the  chamber  (9-11  msec)  are  included. 

From  Figure  42  it  is  evident  that  the  thrust  chamber  configura¬ 
tion  effects  on  ignition  time  are  not  significant.  The  range 
of  ignition  delays  are  quite  similar  for  all  configurations 
except  for  the  50  in.  L*  tests.  Although  the  average  ignition 
times  are  approximately  the  same,  the  shortest  ignition  time 
with  50  in.  L*  chamber  was  11  msec  compared  with  9-10  msec 
for  the  other  L*  configurations,  indicating  the  possibility 
of  a  slight  L*  effect  at  large  L*s.  Also,  close  examination 
of  the  data  indicates  that  the  ignition  times  for  the  200  psia 
design  chamber  pressure  configurations  were  less  than  12.5  msec, 
while  the  ignition  times  for  the  other  design  chamber  pressure 
configurations  were  up  to  14.5  msec."  Although  it  appears  that 
the  average  ignition  times  are  longer  for  the  30  in.  L*  tests 
having  a  contraction  ratio  of  1.5,  similar  tests  with  this 
contraction  ratio  at  other  L*s  do  not  indicate  any  significant 
effect  of  this  parameter. 

Similarly,  there  are  no  signf leant  effects  of 
thrust  chamber  configuration  on  rise  time  as  indicated  in 
Figure  43.  Rise  times  shown  in  the  figure  are  measured  from 
ignition  to  90%  of  steady-state  chamber  pressure.  As  might 
be  expected,  however,  the  rise  times  appear  to  increase  as 
the  L*  of  the  thrust  chamber  is  increased  due  to  the  larger 
volumes  involved,  although  the  effect  is  small  over  the  wide 
range  of  configurations  actually  tested. 

Contrary  to  the  results  of  thrust  chamber  tests 
performed  previously  with  N2O4  antl  hydrazine-type  propellants 
(Ref.  1),  it  is  apparent  that  the  influence  of  thrust  chamber 
configuration  on  ignition  and  rise  time  is  small  with  Compound  A 
and  N2H4.  These  results  are  not  contradictory  however, 
because  it  also  has  been  demonstrated  in  unconfined  tests  that 
ignition  of  Compound  A  and  hydrazine  is  not  influenced 
significantly  by  ambient  pressure.  The  ignition  characteristics 
of  N2C>4/hydrazine-type  propellants  are  a  strong  function  of 
ambient  pressure  and  depend  upon  vaporization  of  the  propellants 
to  increase  the  pressure  in  the  chamber  to  a  level  at  which 
ignition  can  occur.  Although  pressure  also  may  rise  in  the 
chamber  in  tests  with  Compound  A/hydrazine-type  fuels  due  to 
vaporization  prior  to  ignition,  ignition  does  not  depend  on 
pressure  rise  and  there  is  no  direct  correlation  between 
ignition  delay  and  pressure  in  the  chamber  at  ignition. 
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The  short  ignition  delzyr  obtained  with 
Cbnpotind  A  and  tbs  hydrazine  fuels  aluo  resulted  in  anooth 
ignition  transients  in  nost  cases.  No  pressure  spikes  of 
the  eery  short  duration,  high  kaplltuA*  type  associated 
with  M0O4  and  hydrazine  fuels  sere  encountered.  Only  Tests  651, 
892,  915  and  916  had  any  indications,  of  pressure  peeks 
during  the  starting  transient  and  these  sere  lees  than  steady- 
state  c hasher  pressure. 

Three  Compound  A/M4  oscillograph  traces  are 
reproduced  in  Figure  44.  Tests  *16 ,  813,  and  889  sere  node 
at  nosinal  chasber  pressures  ci  200,  79,  and  20  psia» respec¬ 
tively  eith  a  contraction  rat^o  of  9.0  and  an  L*  of  30  in. 

These  traces  are  typical  of  the  smooth  transients  experienced 
with  nost  of  the  Compound  A  tests.  The  upper  trace  in  each 
record  is  the  low  pressure  Eistler  transducer  shoeing  the 
pressure  before  ignition,  fhe  next  trace  is  the  high  pressure 
Mistier  which  is  calibrate 1  for  86,  33,  and  16.5  pei/oa  in  the 
three  records.  The  third  trace  is  the  pho  torn  It  ip  Her  which 
deflects  downward  at  ignition  and  the  lowest  trees  is  the  valve 
current. 


c.  Iflect  of  Mixture  latlo  on  Ignition 

The  effect  of  sixture  ratio  on  ignition  tise 
is  shown  in  Figure  45.  Ignition  tines  are  co spared  for 
COspound  A/N2H4  (Table  XVIII)  and  Compound  A/B3F-5  (Table  XXI) 
at  sixture  ratios  of  1.5,  2.0  and  2.5  with  thrust  chasbers 
having  design  chasber  pressures  of  20,  75  and  200  psia  and 
L*s  of  5  In.,  30  in.  and  50  in.  As  indicated  in  the  figure, 
average  ignition  times  Increased  about  1  asec  as  the  sixture 
rstio  was  increased  to  2.5  frow  the  newinsl  2.0  value  and 
decreased  shout  1.5  asec  as  the  mixture  ratio  was  decreased 
to  1.5.  These  results  are  mitigated  sose what  by  the  effect 
of  varying  propellant  flowrates  upon  entry  tise  into  the 
cbaartK .  As  previously  discussed  is  Section  VCI,  varying 
hydrazine  flowrate  to  obtain  other  Mixture  ratio#  changed 
entry  tise  up  to  about  0.5  sse*  at  20  aw  Mg  pressure.  Entry 
tins  was  shorter  at  the  lower  sixture  rat  lew  and  Imager  at 
the  higher  sixture  ratios  due  to  the  varying  flowrates  through 
the  fixed  seal fold  sad  injector  volumes.  Is  view  of  these 
effects  on  propellant  entry,  however,  it  still  appears  that 
ignition  tines  arm  smaewhat  shorter  with  lower  sixture  ratios, 
•sth  *98*  and  MUF-5  showed  the  sane  trend.  As  is  the  tests 
at  the  sonisal  2.0  sixture  ratio,  thrust  chasber  configuration 
IMS  no  significant  effect  on  Ignition  tine. 
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Test  886 


Pc  -  105  psig 

A  /A.  -  9.00 
c  t 

L*  -  30 


Pc  -  50  psig 
Ac/At  -  7.90 
L*  -  30 


* 


Test  893 

pc  "  5  pai& 
Ac/At  -  7.68 

L*  -  30 


Test  889 

(sweep  -  5  asec/ca) 


Figure  44.  Oscillograph  Records  of  Compound  A/N2H4 
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Figure  45.  Effect  of  Mixture  Ratio  on  Ignition  Tine 
Compound  A/H2H4  and  Compound  A/MHF-5 
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The  apparent  effect  of  mixture  ratio  on 
ignition  time  was  also  observed  in  the  previous  tests  with 
**2®4 /hydrazine  fuels  and  is  discussed  further  in  Reference  1. 

d.  Effect  of  Propellant  Leads  on  Ignition 

The  effect  of  propellant  leads  on  the  ignition 
of  Compound  A/N2H4  is  shown  in  Figure  46.  The  center  point 
of  the  abscissa  represents  simultaneous  application  of  power 
to  both  propellant  valves.  Under  these  conditions  the 
schlieren  characterization  of  the  injector  discussed  in 
Section  VCI  indicates  that  fuel  first  enters  as  a  liquid  at 
10-11  msec  while  Compound  A  first  enters  as  a  vapor  at 
3  msec,  a  mixture  of  vapor  and  liquid  at  7  msec,  and  finally 
as  a  liquid  at  10-11  msec.  For  the  tests  with  simultaneous 
energization  of  both  propellant  valves,  it  is  apparent  that 
ignition  occurs  at  or  just  after  both  propellants  arrive  in  the 
chamber  in  the  liquid  phase.  Since  a  vapor  phase  of  N2H4 
not  detectable  in  the  schlieren  films,  this  time  also  corres¬ 
ponds  to  the  initial  entry  of  NoUg-. 

A  fuel  lead  was  created  by  delaying  the  oxidizer 
valve  with  the  variable  delay  rheostat  in  the  control  system. 
Under  these  conditions,  the  fuel  entry  time  was  not  affected 
as  indicated  by  the  horizontal  band  in  Figure  46  corresponding 
to  a  fuel  entry  time  of  10-11  msec  to  the  left  of  center  of 
the  figure.  If,  for  example,  a  fuel  lead  of  3  msec  were 
desired,  energization  of  the  oxidizer  valve  would  be  delayed 
3  msec  and,  as  indicated  in  the  figure  for  a  3  msec  fuel  lead, 
oxidizer  vapor  first  appears  at  6  msec,  the  mixture  at  10  msec 
and  liquid  at  13-14  msec.  Conversely,  an  oxidizer  lead  was 
obtained  by  delaying  energization  of  the  fuel  valve.  Note  that 
entry  times  of  the  various  oxidizer  phases  under  these  conditions 
are  represented  by  horizontal  lines  or  bands  to  the  right  of 
the  figure. 


Data  points  for  the  three  design  pressure  con¬ 
figurations  shown  in  the  legend  have  been  superimposed  on  the 
propellant  entry  times.  It  is  evident  that  with  an  oxidizer 
lead,  ignition  does  not  occur  until  entry  of  the  liquid  fuel. 
Also,  with  a  fuel  lead,  ignition  again  occurs  when  liauid  fuel 
first  enters  regardless  of  the  physical  state  o'f  the  oxidizer. 
With  a  3  msec  fuel  lead  the  oxidizer  is  a  mixture  of  vapor  and 
liquid  at  ignition  and  with  a  5  msec  fuel  lead  the  oxidizer  is 
in  the  vapor  state.  The  hydrazine  is  the  controlling  propellant 
and  ignition  occurs  upon  its  injection  in  the  liquid  state.  If 
a  vapor  phase  preceeds  liquid  injection  it  cannot  be  detected 
by  the  schlieren  system  and  is  too  small  to  cause  ignition. 
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e.  Triplet  Injector  Teats 

Three  chamber  configurations  were  teste.d  with 
Compound  A/hydrazine  using  the  triplet  injector  shown  in 
Figure  35.  The  results  of  these  tests  are  summarized  in 
Table  XX.  In  this  injector,  the  oxidizer  passages  were  similar 
to  those  in  the  doublet  injector,  but  the  fuel  passages  were 
necessarily  longer  due  to  the  orifice  pattern.  As  could  be 
expected  from  the  results  with  the  doublet  injector,  the 
ignition  times  were  about  4  msec  lodger  with  the  triplet  in¬ 
jector  due  to  the  later  fuel  entry  associated  with  the  longer 
manifold.  The  pressure  at  ignition  was  higher  due  to  the 
longer  oxidizer  lead  but  times  to  90%  Pc  were  comparable  to 
those  with  the  doublet  injector. 

3.  Thrust  Chamber  Tests  -  Compound  A/UDMH,  Compound  A/ 

llHZ-ii 


Thrust  chamber  tests  were  also  made  with  Compound  A/ 
UDMH  and  Compound  A/MHF-5  for  comparison  with  the  Compound  A/ 
K2H4  tests  previously  discussed.  These  tests  were  made  with 
the  same  thrust  chamber  configurations  using , the  doublet 
injector  shown  in  Figure  35.  The  data  for  the  Compound  A/UDSSH 
and  Compound  A/MHF-5  tests  are  tabulated  in  Tables  XXI  and 
XXII,  respectively.  Individual  characteristics  of  these  com¬ 
binations  first  will  be  discussed  and  then  comparisons  with 
Compound  A/N2H4  will  be  made. 

a.  Compound  A/UDME  Tests 

In  tests  855,  856,  857  and  863  with  Compound  A/ 
UDMH,  the  initial  photomultiplier  indication  of  ignition  was 
an  almost  vertical  deflection  at  approximately  7  msec  which 
left  the  oscilloscope  face  in  a  small  fraction  of  a  milli¬ 
second.  The  ignition  times  in  these  four  tests  correspond  to 
the  entry  time  of  the  mixture  of  vapor  and  liquid  fuel  as 
determined  by  the  schlieren  tests  reported  in  Table  XVI.  In 
the  other  five  tests,  one  or  more  weaker  photomultiplier  sig¬ 
nals,  with  a  return  to  or  near  zero  between  pips,  preceeded 
the  stronger  signal.  The  weak  signals  correspond  to  the 
period  of  gaseous  injection  from  3  to  7  msec  as  indicated  by 
the  schlieren  films  and  the  stronger  signal  at  7  msec  corresponds 
to  entry  of  the  fuel  as  a  vapor/liquid  mixture. 

The  occurrence  of  the  weak  reactions  during  fuel 
vapor  injection  appears  to  be  influenced  by. thrust  chamber 
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configuration.  The  thrust  chambers  in  which  the  vapor  phase 
reactions  were  detected  (tests  858  to  861)  were  designed  for  low 
chamber  pressure  and  had  relatively  large  chamber  volumes  and 
diameters.  The  tests  with  thrust  chambers  designed  for  high 
chamber  pressure  which  had  snail  volumes  and  diameters  (for  the 
same  L,*  and  Ac/At)  by  comparison,  however,  did  not  indicate  a 
vapor  phase  reaction  and  ignition  did  not  occur  until  the  UDMH 
was  injected  as  a  vapor/liquid  mixture. 

The  schlieren  test  shown  in  Table  XXI  with  both 
propellants  injected  simultaneously  showed  oxidizer  vapor  appear¬ 
ing  at  3.0  msec,  fuel  vapor  at  3.1  msec,  and  ignition  occurring 
at  3.3  msec.  In  the  same  test,  the  photomultiplier  tube  also 
indicated  ignition  at  3.3  msec  by  slight  deflection  of  the  trace. 
A  stronger  signal  occurred  at  7  msec  when  the  liquid  fuel 
entered. 


The  times  from  signal  to  90%  of  steady  state 
chamber  pressure  in  the  chamber  tests  were  about  2  msec  shorter 
with  UDMH  than  witn  hydrazine  due  to  the  2  msec  shorter  ignition 
times  obtained  with  UDMH.  The  shorter  ignition  times  in  turn 
were  due  to  the  shorter  time  required  for  UDMH  to  enter  the 
chamber  compared  with  N2H4 .  Rise  times  from  ignition  to  90%  Pc, 
however,  were  comparable  for  hydrazine  and  UDMH. 

b.  Compound  A/MHF-5  Tests 

The  tests  ms  je  with  Compound  A/MHF-5  are  summarized 
in  Table  XXII.  The  minimum  volume,  doublet  injector  was  used 
with  three  chamber  configurations  at  mixture  ratios  of  1.5,  2.0 
and  2.5.  Several  tests  (986-  992)  were  made  without  the  chamber 
installed  in  order  to  obtain  ignition  characteristics  with  the 
schlicren  system. 


In  tests  987,  988,  and  989,  ignition  occurred  10.1 
to  11.3  msec  after  valve  signal,  corresponding  to  the  liquid 
entry  time  of  liquid  MHF-5  (10-11  msec) .  The  ignition  times  also 
corresponded  to  the  entry  time  of  liquid  Compound  A  (10-11  msec) 
as  also  determined  in  previous  schlieren  films.  Since  ignition 
times  corresponded  to  the  liquid  entry  times  of  both  propellants, 
Tests  989  to  991  were  made  with  long  propellant  leads  (about  1/4 
sec)  to  determine  if  ignition  would  occur  with  either  propellant 
in  the  vapor  phase.  In  Test  989  with  a  long  fuel  lead  in  order 
to  establish  liquid  fuel  injection  before  entry  of  the  oxidizer, 
ignition  occurred  less  than  3  msec  after  the  "on"  signal  to  the 
Compound  A  valve.  Ignition,  therefore,  coincided  with  initial 
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gaseous  injection  of  the  Compound  A.  Xn  Tests  990  and  991  with 
long  oxidizer  leads,  to  establish  liquid  oxidizer  injection,  the 
ignition  tines,  measured  from  "on"  signal  to  the  fuel  valve,  were 
8.1  and  8.8  nsec,  which  correspond  very  closely  to  the  vapor/ 
liquid  mixture  tines  for  MHF-5.  It  appears  therefore  that  if 
either  the  Compound  A  or  the  MHF-5  is  present  as  a  liquid,  ignition 
will  occur  with  the  other  propellant  in  the  liquid  or  vapor 
phase  (vapor /liquid  mixture  for  MHF-5)  with  ignition  delays  of 
less  than  1  msec. 

The  other  tests  in  Table  XXII  with  the  three  chamber 
configurations  have  ignition  delays  within  2  msec  of  the  delays 
obtained  in  the  unconfined  tests  for  the  same  mixture  ratio. 

As  in  the  hydrazine  tests,  chamber  configuration  does  not  affect 
ignition  time  and  the  delay  in  the  unconfined  tests  are  comparable  | 

to  the  shortest  delays  in  the  chambers.  The  variation  in  ignition  | 

time  with  mixture  ratio  has  already  been  3hown  in  Figure  45.  jj 

In  Tests  965-980,  a  third  Kistler  transducer  was  ! 

installed  at  the  injector  end  of  the  chambers  which  were  about  | 

8  in.  long  in  this  series  of  tests.  In  the  8  tests  with  the  I 

75  psia  design  chamber  pressure  nozzle  (965  -  972) ,  there  was  | 

no  indication  of  ignition  from  the  photomultiplier  tube  until 
the  strong  vertical  deflection  corresponding  to  liquid  entry  of 
the  propellants.  In  all  of  these  tests,  the  pressure  transducer 
at  the  injector  end  of  the  chamber  showed  a  sharp  rise  at  the 
same  time  as  the  photomultiplier  trace  deflected.  The  two  trans-  \ 

ducers  at  the  nozzle  end  of  the  chamber,  however,  did  not  start 
to  rise  until  1/2  msec  later.  The  transducer  outputs  were  inter¬ 
changed  on  the  oscilloscope  in  two  tests  to  be  sure  the  1/2  msec  . 

difference  was  not  a  result  of  a  time  lag  in  the  dual  beams  of 
the  oscilloscope.  The  same  time  lag  was  obtained  in  each  case  f 

indicating  that  the  delay  in  pressure  rise  between  the  injector  jj 

end  and  the  nozzle  end  of  the  chamber  was  real.  Once  ignition  ! 

occurred,  however,  there  was  no  detectable  phase  shift  in  the  ( 

oscillations  of  the  transducer  at  either  end  of  the  chamber.  jj 

I 

Three  of  the  8  tests  with  the  20  psia  design  - 

chamber  pressure  nozzle  (Tests  973,  974,  977)  had  an  indication  ] 

of  the  preignition  reactions  discussed  previously  on  the  photo-  | 

multiplier  trace.  These  reactions,  however,  were  not  detectable 
in  the  pressure  traces.  In  these  three  tests,  the  injector 
pressure  rise  was  delayed  1/2  msec  after  the  stronger  photo¬ 
multiplier  signal  although  the  reason  for  the  delay  is  not 
evident.  The  remaining  tests  were  similar  to  the  75  psia  chamber 
pressure  tests  in  which  injector  pressure  rise  corresponded  to 
the  strong  deflection  of  the  photomultiplier  trace,  and  up  to 
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1/2  msec  delay  in  the  pressure  rise  at  the  nozzle  end  was 
observed.  These  data  corroborate  the  photographic  evidence 
previously  noted  that  ignition  of  Compound  A/N2H4  occurs  at 
the  injector  and  takes  a  finite  time  to  progress  to  the 
nozzle  end  of  the  chamber. 

c.  Comparison  of  Compound  A/Hydrazine  Fuels  Test 

The  ignition  times  and  the  times  to  90%  of  steady- 
state  chamber  pressure  are  compared  for  Compound  A  and  the  three 
hydrazine-type  fuels  in  Figure  47  and  48.  The  data  are  for 
thrust  chamber  configurations  having  an  L*  of  30  in.  and  a 
contraction  ratio  of  3.5  for  each  of  the  three  design  chamber 
pressures.  The  tests  were  made  at  a  mixture  ratio  of  2.0. 

As  indicated  above,  inspection  of  Figure  47  shows 
that  UDMH  has  the  shortest  Ignition  time  of  the  three  fuels. 
Ignition  occurs  at  about  7  msec  corresponding  to  the  entry  of 
a  liquid/vapor  mixture  of  UDIIH.  With  N2H4  and  MHF-5,  ignition 
does  not  occur  until  after  10  msec  at  which  time  both  propellants 
are  being  injected  in  the  liquid  phase.  As  indicated  in  the 
figure,  design  chamber  pressure  has  no  significant  effect  on 
the  ignition  times  of  the  respective  combinations.  Although 
not  shown  in  this  figure,  previous  discussions  have  also  shown  that 
the  other  thrust  chamber  configurations  (L*,  Ac/At)  have  little 
effect  on  ignition  time. 

As  shown  in  Figure  48,  90%  of  steady-state 
chamber  pressure  is  also  reached  faster  with  UDMH  than  the 
other  fuels  by  virtue  of  the  earlier  ignition  obtained  with 
UDMH.  Rise  times  from  ignition  to  90%  of  steady-state  chamber 
pressure,  however,  are  similar  for  all  three  fuels.  As  might 
be  expected,  the  rise  times  associated  with  the  low  design 
chamber  pressure  conf igurations  are  slightly  longer  than  with 
the  high  design  chamber  pressure  configurations  because  of 
the  larger  chamber  volumes  involved. 

The  ignition  transients  with  both  UDMH  and  MHF-5 
were  smooth  and  there  was  no  indication  of  pressure  peaks  or 
spikes  in  any  of  the  tests. 
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D.  Xxpcrlme  ntal  Besults  -  Fluorine/Hydrogen  T»»t« 

A  complete  Mrivs  of  ignition  tests  ess  also  Bade  with 
gaseous  fluorine  and  hydrogen  using  the  thrust  c hasher  con¬ 
figurations  shown  in  Table  XV.  The  chambers  and  nosales  were 
identical  to  those  used  for  the  Co«pound  A  tests,  but  the 
concentric  injector  shown  in  Figure  36  was  used  with  the  larger 
propellant  waives  previously  discussed  in  order  to  obtain  the 
desired  flowrates  with  gaseous  propellants.  The  fluorine  was 
injected  axially  through  a  single  orifice  and  the  hydrogen  was 
injected  through  a  concentric  annulus  so  that  it  inpinged  on 
the  fluorine  stress  at  a  30°  angle. 

Data  fron  the  tventy-seven  chamber  configurations  tested 
are  suaaarised  in  Tabl^  XXIII.  The  tabulated  parameters  are 
defined  as  they  were  for  the  Coapound  A  tests.  Ignition  tine 
shown  in  the  table  is  based  on  the  photonultiplier  output.  Ho 
preignition  reactions  were  noted  with  this  propellant  combina¬ 
tion.  The  shorter  i^nitiou  and  rise  tines  experienced  with  this 
combination  permitted  a  faster  oscilloscope  sweep  which  resulted 
in  ignition  time  resolution  of  0.1  msec  compared  with  0.5  nsec 
for  the  Compound  A/HgHg  testa.  The  "valve  current"  column  in 
the  table  is  the  tine  from  valve  "on"  signal  to  the  first 
minimum  in  the  valve  current  trace.  Six  tests  were  made  vith 
the  fuel  valve  delayed  2-3  msec  as  indicated  in  the  table. 

Early  in  the  test  program  it  was  found  that  a  2-3  nsec 
spread  in  ignition  tine  was  being  experienced  which  could  not 
be  correlated  with  thrust  chamber  configuration.  The  ignition 
time,  which  is  defined  aa  the  time  Iron  valve  "on"  signal  to 
ignition,  varied  from  2.8  to  6.0  for  tests  vith  simultaneous 
propellant  valve  energization.  A  spread  of  up  to  2.3  nsec  was 
noted  for  successive  tests  vith  the  sane  chamber  configuration. 
In  order  to  determine  the  cause  of  these  small  but  unexplained 
variations  in  tbo  results,  the  various  components  of  the  test 
systen  were  examined.  It  vaa  found  that  the  DC  power  supply 
cced  for  these  toats  had  a  120  cps  ripple  when  drawing  the  12 
anps  required  to  operate  the  larger  propellant  valvea.  This 
ripple  was  visible  on  the  valve  current  trace.  The  effect  of 
the  r 'pp?e  was  somewhat  accentuated  by  the  fact  that  the  tank 
pressures  required  to  obtain  the  desired  gaseous  propellant 
flowrates  were  approaching  the  limiting  pressure  at  which  the 
large,  direct-acting  eolesoid  valvea  would  open  consiptently . 

The  variation  is  power  to  the  ripple  at  the  critical  openly 
point  could  effect  actual  opening  time  by  several  msec,  it  is 
evident,  therefore,  that  the  apparent  variation  is  ignition  tines 
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encountered  is  due  to  the  ripple  rather  than  the  result  of  any 
propelHant  or  chamber  characteristics.  In  Figure  49 ,  the 
ignition  time  is  plotted  as  a  function  of  the  minimum  point  on 
the  valve  current  trace.  It  is  apparent  that  all  ignition  times 
can  be  Included  in  a  1  msec  wide  band  irrespective  of  chamber 
configuration  and  that  with  a  more  stable  power  supply,  all 
chamber  configurations  would  have  had  equal  ignition  times  within 
+1/2  msec. 

Schlieren  films  also  were  taken  of  the  injection  and  ignition 
of  the  unconfined  gaseous  propellants.  Both  the  fluorine  and 
hydrogen  gases  were  difficult  to  detect  in  the  schlieren  pictures 
at  ambient  pressures  below  10  mm  Hg.  Because  of  the  low  density 
of  the  gases  it  is  possible  that  some  vapors  may  have  entered 
for  a  short  period  before  the  flowrate  was  sufficient  to  detect 
propellant  entry  at  the  low  ambient  pressure.  The  propellant 
entry  times  in  the  four  Fo  tests,  three  H2  tests,  and  three 
F2/H2  schlieren  ignition  tests  are  included  in  Figure  49  and  are 
all  within  1/4  msec  of  the  center  of  the  1  msec  spread.  All 
three  ignition  times  of  these  unconfined  tests  are  within  the 
1  wove  bauu  also.  In  tne  unconfined  tests  with  both  propellants 
being  injected,  the  leading  propellant  obscured  entry  of  the 
second  propellant  because  of  the  concentric  design.  However, 
ignition  delays  from  the  leading  propellant  were  only  0.70  msec 
at  100  mm  Hg  ambient  pressure,  0.29  at  50  mm,  and  0.24  at  7  mm. 

The  apparently  longer  delay  at  100  mm  may  have  been  due  to  more 
accurate  detection  of  the  first  gas  at  the  higher  pressure  rather 
than  a  truly  longer  delay.  Ignition  was  clearly  visible  near 
the  injector  in  the  films  at  50  and  100  mm  Hg.  At  7  mm  Hg, 
however,  the  origin  of  ignition  could  not  be  determined  and  it 
may  have  occurred  out  of  the  field  of  view  of  the  schlieren 
system  which  is  limited  to  about  2  inches  from  the  impingement 
point . 

Several  thrust  chamber  tests  were  made  with  the  third 
Kistler  transducer  at  the  injector  end  of  the  chamber.  As  in 
the  Compound  A  tests,  the  initial  rise  in  chamber  pressure  at 
the  injector  end  coincided  with  the  initial  visible  light.  No 
delay  in  chamber  pressure  rise  at  the  nozzle  end  was  experienced 
with  the  shorter  (1  -  4  in.)  chambers.  However,  delays  in 
pressure  rise  at  the  nozzle  end  with  longer  chambers  ranged  from 
0.4  to  0.6  msec  for  chambers  about  12  in.  long  to  a  maximum  of 
0.8  msec  for  the  22  in.  long  chamber. 

The  rise  times  from  ignition  to  90%  P~  were  less  than 
7.5.  msec  except  for  one  configuration  (Tests  1076-1078)  which 
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does  not  seem  to  fit  the  other  data.  Rise  times  for  these 
three  tests  were  unusually'  long  for  no  apparent  reason.  Except 
for  these  tests,  rise  times  were  proportional  to  chamber  volume 
or  L*  as  shown  in  Figure  50.  With  the  30  and  50  in.L*  con¬ 
figurations,  the  rise  time  data  spread  was  +1  msec  and  there 
was  no  significant  effect  due  to  either  contraction  ratio  or 
design  chamber  pressure.  Data  spread  was  slightly  larger  at 
the  smaller  L*'s  and  the  lowest  design  chamber  pressure  (20  psia) 
appeared  to  have  shorter  rise  times.  However,  combustion 
efficiency  with  the  smaller  L*,  20  psia  chambers  was  low  and 
the  time  required  to  reach  the  low  steady-state  pressure  obtained 
in  these  tests  is  undoubtedly  less  than  that  required  to  reach 
higher  pressures  consistent  with  good  performance. 

Six  tests  (Nos.  1008,  1011,  1002,  1022,  998,  and  1017)  are 
included  in  Table  XXII 1  in  which  the  fuel  valve  was  delayed  for 
2-3  msec.  Ignition  times  were  delayed  corresponding  to  the 
delay  in  fuel  entry.  Rise  times  were  comparable  to  the  rise 
times  in  the  tests  with  simultaneous  valve  energization  for  the 
10  in.  L*  chambers  and  1-2  msec  shorter  for  the  50  in.  L*  chambers. 
Some  pressure  rise  due  to  the  fluorine  could  be  detected  in  the 
higher  design  pressure,  shorter  L*  chamber  configurations  when 
the  fuel  valve  was  delayed.  In  all  the  other  tests  ignition 
proceeded  any  detectable  rise  in  chamber  pressure. 

There  was  no  indication  of  any  pressure  spikes  during  the 
rise  time  with  the  F2/H2  combination  in  any  of  these  tests. 

Ignition  delays  were  short  as  indicated  and  the  pressure  tran¬ 
sients  were  smooth. 

The  oscillograph  traces  from  Tests  1086,  1084.  and  1081  are 
shown  in  Figure  51.  The  trace  positions  are  the  ror  tne 

Compound  A  tests  in  Figure  44,  but  the  sweep  is  increased  from 
5  to  2  msec/cm.  As  in  the  case  of  the  Compound  A/hydrazine- 
propellant  tests  the  short  delays  and  smooth  pressure  transients 
should  be  noted. 


-166- 


rise  Time,  msec 


A7BPL-7B-65-257 


Characteristic  Length,  h*,  in. 


Figure  50,  Effect  of  Chamber  L*  on  Eise  Time  -  F2/% 

Chamber  Tests 
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Pc  -  132  psig 
Ac/At  -  3.22 
-  30  in. 


Test  1086 


■  Pc  - 

psig 

■  VA«  ■ 

3.50 

ill  I>  -  30 

in. 

Test  1084 


<• 


Pc  -  10  psig 
Ac/At  **  3.54 
L*  -  30  in. 

Test  1081 

(«weep  -  2  nsec/cm)  j 

Figure  51.  Oscillograph  Records  of  F2/H2  Tests  *' 
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OOKOiliaiONS 

(B)  Significant  progress  has  been  made  in  obtaining  a 
fundamental  understanding  of  hypergoiic  ignition  under  vacua* 
conditions.  Although  further  refinement  Incorporating  additional 
ef$e$$s  such  ss  propellant  "flashing"  in  injector  manifolds 
is  ^|H|ifed  to  optimize  the  model,  a  working  mathematical  model 
q^icll  p^tditts  the  ignition  delay  of  /hydrazine- type  propel- 
Xante  with  reasonable  accuracy  has  been  developed,  further, 
the  Ignition  model  forms  the  basis  for,  and  can  be  extended  to, 

9  mathematical  model  which  considers  the  pressure  spiking 

Phenomenon. 

(U)  based  on  the  results  of  the*  experimental  and  analytical 
investigations  performed  under  Part  II  of  Contract  AF04(611>- 
9946  the  fallowing  conclusions  are  drawn  with  regard  to  the 
various  tasks  and  phases  of  the  program. 

A,  T^ab  I  -  Measurement  of  Reaction  Kate 

(U)  1.  At  sub-ignition  pressures,  a  pre-ignition  reaction 

product  (adduct)  is  formed  by  each  of  the  NG^/hydra^ine  fuel 
combinations  tested.  Chemical  analyses  of  the  clear,  yellow, 
viscous  liquid  adduct  formed  in  appreciable  quantities  by 
MME  and  NOj  vapors  indicate  the  adduct  is  an  associative  product 
of  the  rpactants,  has  the  characteristics  of  a  monopropellant , 
and  contains  considerable  energy.  Its  relation  to  pressure 
Spiking  during  engine  start  transients  is  not  as  determined. 

(©)  2,  the  initial  rate  of  pressure  rise  ia  a  thrust  chamber 

is  retarded  in  the  case  of  very  volatile  propellants  such  as 

apd  Compound  A  by  substantial  propellant  "flashing"  within 
the  injector  volume  upon  propellant  valve  opening.  A  significant 
time  elapses  before  the  mass  flowrate  into  the  thrust  chamber 
reaches  its  full,  nominal  value. 

(U)  3.  Although  sufficiently  accurate  chamber  pressurization 

histories  for  K2O4  cannot  be  calculated  at  present  due  to  its 
trgnmlent  flowrate  behavior  upon  propellant-valve  opening,  the 
■' Chemical  kinetics  part  of  the  hyper go lie  ignition  model  can  be 
cpmba^md  with  experiment.  Theoretically  calculated  ignition 
dfjLjky  tipes  in  ah  engine,  based  on  estimated  vapor-phase  com¬ 
position  and  temperature,  agree  reasonably  with  experimentally 
determined  engine  ignition  delays. 
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(U)  4.  In  agreement  with  experimental  results,  the  theoretically 

derived  hypergolic  ignition  w odel  indicates  that  the  optimum 
vapor-phase  composition  for  shortest  ignition  delays  is  the 
equiMOlecular  mixture  which,  for  ITAIH,  MMB  and  50-50  with  N3O4, 
are  markedly  more  fuel  rich  than  typical  operational  mixture 
ratios*  Advantage  of  this  fact  to  obtain  shorter  ignition 
delays  may  be  taken  through  suitable  injector  designs  Coupled 
with  propellant  leads  and  lags, 

(U)  5.  Dominant  ignition  reactions  of  M^Ga/hydrazine-type  fuels 

lew  pressures  are  thermal,  gas  phase  reactions  which  are 
blmolaoular  and  have  low  activation  energies  (less  than  10  kcal/ 
melf)  and  low  pre-exponential  factors. 

(U)  6.  Propellant  vaporization  in  a  thrust  chamber  prior  to 

ignition  is  non-adisbatic  due  to  heat  addition  to  the  vapor-drop 
System  from  the  thrust  chamber  walls.  / 

B.  Task  II  *  Evaluation  of  Additives  to  Seduce  Activation 


(U)  1.  The  low  pressure,  premixed  vapor  stream  apparatus  is 

well  suited  to  screeuing  of  additives.  The  measured  variable 
ignition  pressure  limit  is  a  sensitive,  reliable  measure  of 
ignition  delay.  Additives  may  lie  evaluated  relatively  easily 
with  good  accuracy  and  lees  expensively  than  in  engine  firings 
at  low  ambient  pressure . 

(U)  8.  Of  the  five  organic  fuel  additives  tested  in  the 

apparatus,  furfuryl  alcohol  was  found  to  have  a  significant 
beneficial  effect  on  the  ignition  characteristics  of  N2O4/IMH. 

(0)  3.  The  one  principal  oxidizer  additive  tested,  Cospound  R- 

PC(NFa)St  not  only  worsened  ignition  characteristics  but  also 
reduced  reproducibility  considerably.  Only  after  thoroughly 
washing  the  apparatus  were  the  normal,  repeatable  results  obtained 
again  with  the  neat  propellants. 

C.  iAak  HI  -  Thrust  Chamber  Design  Parameter  Study 


(C)  a.  Ignition  delays  in  the  thrust  chamber  tests  with 
thp  Compound  A/hydyasinp-type  fuels  were  very  short.  Ignition 
decays  with  ODS®  are  somewhat  shorter  than  those  with  N2H4  or 
KDErb,  although  rise  times  (ignition  to  90%  Pc)  are  comparable. 
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((S  b.  Because  of  the  short  ignition  delay,  tbs  pressure 
trkrSlbht*  Mrs  smooth  with  no  indication  of  the  very  short 
duration,  high  amplitude,  xandom  pressure  spikes  characteristic 
of  the  XgQg/hydrasise-type  fuels. 

<C)  e.  Ignition  occurs  at  tht  injector  and  proceeds 
rapidly  and  smoothly  through  the  chaster. 

(0)  d.  Ignition  in  a  thrust  chaster  is  independent  of 
pressure  rise  due  to  propellant  vaporisation.  Ignition  delays 
is  the  chaster  tests  with  each  propellant  were  comparable  with 
the  delays  measured  in  the  unconfine *3  echlieren  tests.  By 
contrast.,  previous  thrust  chaster  tests  with  WaOg/hydrar.in  - 
type  fuels  (itef.  1)  indicated  an  extrese  pressure  dependency. 

(C>  i.  Ignition  characteristics  are  not  significantly 
influenced  by  thrust  chaster  configuration  parameters  because 
of  the  lack  of  pressure  dependency. 

(C)  t.  Although  the  influence  of  thrust  chaster  geometry 
on  rise  tii*e  to  P0%  of  tbs  steady-state  chamber  pressure  was 
ssall,  the  configurations  with  the  pnaller  volumes  generally 
produced  slightly  shorter  rise  tines. 

(v,)  g.  Mixture  ratio  tests  indicated  somewhat  shorter 
ignition  delays  with  fuel-rich  mlxtur  ratios  than  with  oxidizer- 
rich  six tore  radios. 

(C)  k.  Propellant  lend  tests  with  Compound  A/NgR* 
indicate  that  Ignition  occurs  rapidly  upon  entry  of  hydrazine 
regardless  of  the  state  (vapor,  mixture,  or  liquid)  of  Compound 
A. 


2.  Fg/Ha 

(9)  a.  Ignition  oelayc  were  very  short  (within  i  nsec) 
in  the  thrust  chastex  and  uncounted  schl-eren  tests. 

(*f)  b.  As  with  Compound  A/hydrazlne-t?pe  fuels,  ignition 
does  not  depend  upon  pressure  rise  In  the  chaster.  In  most 
tests  ignition  occurred  before  any  dev  actable  rise  in  pressure 
due  to  pro^ella-t  eatry. 

(0)  c.  Ignition  delay  was  not  influenced  by  thrust 
chaster  design  parameters. 
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(Q)  d.  Ignition  occurred  at  the  injector  and  pressure 
treasiaints  sere  smooth  rith  no  indication  of  pressure  spikes 
la  any  of  the  tests. 


(H)  e.  Rise  time  tc  90%  of  steady-state  chamber  pressure 
was  eery  rapid.  Else  time  sas  someshat  shorter  with  smaller 
L*'s  than  *ith  larger  L*’s  as  might  be  erpected. 
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VII 

RECOMMENDATIONS 


Several  factors  pertinent  to  hypergolic  ignition  in  reaction 
control  systems  that  were  either  outside  the  3Cope  of  the  present 
program  or  require  further  study  beyond  that  which  time  permitted 
in  this  six-aontfa  effort  are  enumerated  below: 

1.  The  pre-ignition  reaction  product  (adduct)  appears,  at 
present,  to  be  a  likely  candidate  for  the  cause  of  the  pressure 
spikes  experienced  during  engine  start  transients.  It  should  be 
determined  whether,  in  fact,  the  adduct  is  the  cause.  Properties 
of  the  adducts  formed  by  the  hydrazine  fuels  with  N2O4  are  re¬ 
quired,  including  their  rates  of  formation  and  decomposition. 

It  should  also  be  determined  whether  the  adducts  are  necessary 
intermediates  in  the  reaction  mechanisms  for  ignition  or  Whether 
they  are  side  products  which  tend  to  impede  ignition. 

If  the  adducts  do  prove  to  be  the  major  cause  of  pressure 
spikes,  methods  should  be  investigated  to  vl)  avoid  its  formation 
if  it  is  not  a  necessary  intermediate,  (2)  minimize  its .accumulation 
on  thrust  chamber  walls,  during  ignition  delay  times  (and  during 
tail-off  too)  and/or  (3)  retard  its  decomposition  rate  in  order 
to  lengthen  the  time  over  which  the  accumulated  energy  is  released. 

2.  To  be  made  general,  the  analysis  for  chamber  pres¬ 
surization  due  to  prcpellant  vaporization  requires  inclusion  of 
expressions  tc  account  for  the  initial  transient  flowrate 
behavior  of  very  volatile  propellants  upon  propellant  valve  open¬ 
ing.  Also  required  is  a  general  expression  for  the  heat  addition 
to  the  vapor-drop  system  from  the  thrust  chamber  walls. 

3.  Then,  the  complete  mathematical  model  for  hypergolic 
ignition  in  reaction  control  systems  should  be  further  verified 
by  comparisons  with  experimental  ignition  delays  in  engines  of 
various  chamber  configurations,  thrust  levels,  etc. 

4.  A  mathematical  model  of  pressure  spiking  during  start 
transients  should  be  developed  to  permit  evaluation  of  chamber 
geometry  and  ignition  hydraulics  from  a  pressure  spiking  view¬ 
point.  The  hypergolic  ignition  model  developed  in  the  present 
program  is  a  major  part  of  the  pressure  spiking  model  and  is 
readily  extended  to  include  the  pressure  transient. 
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5 ,  Since  a  reliable,  low  cost  method  is  now  available, 
k tensive  screening  of  chemical  additives  to  shorten  ignition 
slays  is  recommended .  The  sore  costly  tasting  in  engines  at 
low  ambient  pressures  can  then  be  limited  to  only  the  most 
promising  additives. 

6a  As  indicated  above,  thrust  chamber  tests  to  verify  th® 
results  of  the  eirperinents  and  to  determine  the  significance  of 
the  results  under  actual  operating  conditions  are  required.  It 
is  recommended,  therefore,  that  a  concurrent  fundamental  investi¬ 
gation  end  thrust  chamber  verification  program  be  undertaken  to 
accomplish  the  above  goals. 
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is.  abstract  (U)  By  theoretical  and  experimental  means,  a  mathematical 
model  of  hypergolic  ignition  in  reaction  control  systems  has  been 
developed  based  on  physical  kinetics  of  droplet  evaporation  and  chem¬ 
ical  kinetics  of  ignition  reactions.  Reasonable  agreement  between 
theoretical  and  experimental  ignition  delays  are  obtained.  The 
dominant  ignition  reactions  of  N2O4  with  the  hydrazine  fuels  at 
reduced  pressures  are  found  to  be  thermal,  gas  phase  reactions  which 
are  bimolecular  and  have  low  activation  energies  and  pre-exponential 
factors.  A  pre- Ignition  reaction  product  was  found  which  is  a  clear, 
yellow,  viscous  liquid  with  a  very  low  vapor  pressure.  The  "adduct" 
has  the  characteristics  of  a  monopropellant  and  contains  considerable 
energy.  •  Its  relation  to  pressure  spiking  during  engine  start  tran¬ 
sients  is  as  yejt  undetermined.  Of  six  additives  tested,  only 
furfuryl  alcohol  had  a  significant  beneficial  effect  on  the  ignition 
characteristics  of  N2°4/MMH. 

The  influence  of  thrust  chamber  design  parameters  on  the  igni¬ 
tion  delay  and  pressure  transients  of  Compound  A  with  N2H4,  UDMH  and 
MHF-5  and  also  gaseous  F2/H2  was  Investigated. 
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